
Documentation, Analysis and Modelling of 
Eucalyptus based system for development of 
consolidated models at Zonal / National level 

Thesis 




Sobmitted for the award of 


Doctor of Philosophy 


AGROFORESTRY 

(Tree Mensuration/Modelling) 


Institute of Agricultural Sciences 

Bundelkhand University, Jhansi 


Under the Supervision and Guidance of 


Senior Scientist 


+33 -H 


National Research Centre for Agroforestry 

(NRCAF) 

Gwalior Road, Jhansi 


Declaration 


I hereby certify that the thesis entitled “Documentation, Analysis and Modelling of 
Eucalyptus based system for development of consolidated models at Zonal/National 
level” being submitted for partial fulfillment of Degree of Doctor of Philosophy 
(Agroforestry), Bundelkhand University, Jhansi (UP) is an original piece of work done by 
me under the supervision and guidance of Dr. Ajit, Senior Scientist, NRCAF, Jhansi and 
does not incorporate, without acknowledgment, any material previously submitted for a 
degree or diploma in any university, and that to the best of my knowledge and belief it 
does not contain any material previously published or written by another person except 
where due reference is made in the text. 


Dated: 



(NighopJabeen) 




ICAR 


Gram : KRISHIVANIKI 
Itkphmt CO): §517- 2730214 
-2730213 
- 2730470 
Fax : 91-IB17-2730364 
E.miil ; iirwIt^liiibLiiiciii 

^ 3rg?T€noT 

National Research Centre For Agro-Forestry 

PahuJ Dam, Jhansi Gwalior Road, >75^ g i f ^ Tmj ps gRqpf 

JHANSI-284003 (U.P.) India sM-284003 (g.s.) iTHcf ' 



NRCAF 


Dated: 0<^-|0-2oo2> 

Dr. Ajit 
Senior Scientist 

Email: umaajitgupta@yahoo.co.in 
ajit@nrcaf.emet.in 


CERTIFICATE 

It is certified that this thesis entitled “Documentation, Analysis and Modelling of 
Eucalyptus based system for development of consolidated models at Zonal/National 
level” is an original piece of work done by Ms. Nighat Jabeen under my supervision and 
guidance for the degree of Doctor of Philosophy in Agroforestry, Bundelkhand 
University, Jhansi (India). 

I further certify that: 

> It embodies the original work of candidate herself. 

> It is upto the required standard both in respect of its contents and literary 
presentation for being referred to the examiners. 

y The candidate has worked under me for the required period at National Research 
Centre for Agroforestry, Jhansi. 

> The candidate has put in the required attendance (more than 200 days) in the 
department during the period. 

Senior Scientist 
NRGAF, Jhansi 




< - ct i {^*iy ^inf, sff# - 284 003 (’3.1.) 

National Research Centre for Agroforestry 

Neai' PatmJ Darn, Gwalior Road, Jhansi - 284 003 (U.P.J Indtei 
Phoije& : 0510 - 2730214 {Oi. 2730366 fR}. Fax : 0510 - 2730364. 

E inaO ; kiislitv'a.atki@iriall.r»rcaf.eraei,.tii . Gram : KRISHIVANIKI 

SI. ma*. HTPli 

P»^!>rc(i 

Dr. S.a. OKyani 

Director 






Dated:09/ 10/2008 


Certificate 


This is to certify that the thesis entitled “Documentation, Analysis and Modelling of 
Eucalyptus based system for development of consolidated models at Zonal/National 
level” submitted by Ms. Nighat Jabeen to the Department of Agroforestry, Institute of 
Agricultural Sciences, Bundelkand University, Jhansi for partial fulfillment for the award 
of the degree of Doctor of Philosophy in Agroforestry is a record of bonajide research 
work carried out by her at this centre under the supervision and guidance of Dr. AJit 
(Senior Scientist). Ms. Nighat Jabeen has worked on this problem for a period of more 
than two years and the thesis is in my opinion worthy of consideration for the award of 
the degree of Doctor of Philosophy in Agroforestry in accordance with the regulation of 
this Centre. The results embodied in this thesis have not been submitted to any other 
university or institution for the award of any degree or diploma. 



Dr. V.K.Gupta 
Principal Scientist and 
Chairman of Higher Studies, 
NRCAF, Jhansi 



NRCAF, Jhansi. 

kJKMKi^ln; 
iinUfwjl llBaaftfe fiessUe fw 
.MSMMOOa 


Bundelkhand University, Jhansi 



Department of Agroforestry 




CERTIFICATE 

This is to certify that thesis entitled “Documentation, Analysis and Modelling of 
Eucalyptus based system for development of consolidated models at Zonal/National 
level” has been carried out by Ms. Nighat Jabeen at National Research Centre for 
Agroforestry, Jhansi under the Supervision of Dr. Ajit, Senior Scientist and is submitted 
in partial fulfillment for the requirements of the Degree of Doctor of Philosophy in 
Agroforestry. 



Head of Department 
Institute of Agriculture Sciences 
Bundelkhand University, Jhansi 


0)r. CoS. Dhanai) 

Head 

Division of Agroforestry 
Bundelkhand University, Jhansi 



Acknowledgement 

There are lots of people I would like to thank for a multiplicity of reasons. 

Firstly, I would like to thank my Supervisor, Dn jIJit (Senior Scientist, NRCAF, Jkamif I could not have imagined hmmg 
a better advisor and mentor for my PkD, and without his common-sense, knowledge, perceptiveness and cracking-ofthe- 
whip, I would never have finished. Not only a great mentor, but he has also been a comers tone in my professional 
development Words fall short to describe what I learned from him. My deepest personal regards are due for him forever. 

I am also indebted to Dn S. K Dhymi, Director, NRCAF, Jhansi, for providing me all the necessary facilities to complete 
this study in time. Thanks are due to Dr* KK.Gupta (Principal Scientist) for permitting me to work at NRCAF, Jhansi, 

It is my profound privilege to verbalize my feelings of gratitude to the members of my PkD, advisory committee namely 
Dr* O* P* Chaturvedi (Principal Scientist) and Dr* AL K Handa (Senior Scientist), who monitored my work and took effort 
in reading and providing me with valuable comments on this thesis. Thanks are due to Dr* AnU Rm (Senior Scientist, 
lASRI, New Delhi) for permitting me to host Tree Growth Modelling website at lASRl server. New Delhi in time, 

I also express my gratitude and indebtedness towards Dr* C S* Dhanai and Dr* G* S* Pawar Head, Department of 
Agricultural Sciences, Bundelkhand University, Jhansi for their support to carry out this research work 

I also wish to thank the scientific galaxy and administrative staff of NRCAF, particularly Dr* Rajeev Tmari, Mr, Rajenika 
Singh and Mr* S* G* Dhamdere for their prolific cooperation, inspiration, generous support, encouragement and countless 
ideas throughout my research work, when ever and where ever I needed, 

I would like to express my heartiest thanks to Sadiya for assisting in thesis typing and her family members for their 
affection and for never letting me feel that I am away from my parents. 

I feel a deep sense of gratitude for my late father who formed part of my vision and taught me the good things that really 
matter in life. The happy memories of my father still provides a persistent inspiration for my journey in this life. 
Unfortunately, he is not here to share the moments with me, but I am sure he keeps an eye on me from heaven. 

1 express my sincere thanks and immense gratitude to my maternal uncles Mr* Ab. Rashid Bulla, Haji Khurshid Ahmad 
Bulla and Mr* Zahoor Ahmad Bulla for their continued encouragement, affection and financial support during my study 
period. 

Great reverence and many many thanks to Mrs* Uma (Technical Officer, NRCAF) who appeared to be my guardian here,! 
truly appreciate her expert guidance during awestruck times. 

My heartfelt thanks to Dr* Irshad (Postdoctoral Fellow, University of Amsterdam) for his quick responses to my emails and 
ready to help attitude. Although we met only for a short period of time, but his enthusiasm for research always inspires me. 

Here, I would like to take an opportunity to acknowledge my friends, Fozia, Dr* Ruchi, Dr* Ritu, Munazah, Snober and 
Mussarat for inspiring and motivating me throughout my research work Warm thanks are extended to my roommates 
Swad and Chanchal (Senior Research Fellows) for good suggestions for this thesis and funny time together, which re- 
energized me at the times I needed most. 

With core of my heart, I am thankful to my dear mother, for her unconditional love and prayers; sister^s Neelofar, Rehana 
and Nighat whose continuous encouragement and support have been a source of inspiration in completion of this tough 
task 

Thanks are due to Arvind and Ashok for assisting in lab work andPremnarayanfor assisting in fieldwork 

Last but not the least; I am thanfful to the Almighty, who has provided me sound health and strength to complete this work 
Whatlamandwhatiwouldbe, I owe to the Almighty for leading me the path of success. 


Perhaps, I forgot someone... so, just in case: thankyou to whom it concerns! 


(Nighat Jabeen) 


'Fr^nc* 


Preface 

The present work is the outcome of the research carried out by me in the field of 
“mathematical modeling” of tree growth in agroforestry ” registered at the Institute of 
Agricultural Sciences. Bundelkhand University Jhansi, under the guidance and supervision 
of Dr. Ajit, Senior Scientist at National Research Centre for Agroforestry (NRCAF), Jhansi, 
as a research scholar. 

The work presented in the thesis entitled “Documentation, Analysis and Modelling of 
Eucalyptus based system for development of consolidated models at Zonal/National level” 
is based on my following research papers [copy of acceptance/publications attached in 
AnnexureJ. 

1. “Generalized models for prediction of aboveground biomass for Eucalyptus species in 
India” 

(Abstract accepted in Proceedings of National Symposium on “Agroforestry Knowledge for 
Sustainability, Climate Moderation and Challenges being conducted during /5**-i/* 
December, 2008 being organized by National Research Centre for Agroforestry, Jhansi. 
India. 

2. “Designing and populating Tree Growth Modelling Website in Indian Perspective” 

(Abstract accepted in Proceedings of National Symposium on “Agroforestry Knowledge for 
Sustainability, Climate Moderation and Challenges being conducted during I5‘^-I7^ 
December, 2008 being organized by National Research Centre for Agroforestry, Jhansi, 
India. 

3. “Comparison of initial growth and yield of Eucalyptus tereticomis clones under agroforestry 

system vis-a-vis block plantations in semi-arid conditions of India”. 
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(Paper was presented in the National Symposium on ‘‘Agroforestry for livelihood security. 
Environmental Protection and Biofuel Production ” held at NRCAF, Jhansi during 16-18 
Dec., 2006. In the ABSTRACT of National Symposium, 125-126 pp. 

4. “Statistical prediction of height through height-diameter models for Eucalyptus tereticomis 
in central India ”. 

(Abstract published in Proceedings of National Symposium on “Intensive Forest 
Farming: The State of the Art”, held during Feb 12-14, 2008 organized by Punjab 
Agricultural University (PAU), Ludhiana. 146-147 pp). 

5. “Height estimation model for Eucalyptus tereticomis grown under semi arid conditions of 

India”. 

(Indian Journal of Agroforestry (2008). 10 (2): 34-39). 

6. “Growth model of Eucalyptus tereticomis and its validation through SPSS: a case study 
from semi-arid conditions in India ”. 

(In Press: SPSS Analyst (2008) 


A broad sketch of the thesis presentation is as follows. 

Chapter-1 entitled ‘Introduction’ primarily defines the term ‘agroforestry’, its various 
definitions, types (viz agrisilviculture, compact block plantations, agrihorticulture etc.), and^ 
its taxonomic classification. Moreover, the significance of Eucalyptus species has also been 
highlighted (a tree of great economic importance, as it is the basis of matchbox, pulp and 
paper production industry) along with its utilities (Eucalyptus trees are valued for their oil, 
gum and timber). Eucalyptus oil is used in medicine to treat many sicknesses such as 
infections, colds, flu, sore throats, bronchitis, and even some skin infections. The timber is 



claimed to be the strongest, densest and most durable in the world, used for the building of 
ships, railway ties, piers and telegraph poles. Eucalyptus oil is readily distilled from the 
leaves and can be used for cleaning, deodorizing, and in very small quantities in food 
supplements, especially cough drops and decongestants. Eucalyptus has been widely 
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accepted by the farmers not only for pure plantations but also in agroforestry systems since 
Eucalyptus is one of the fastest graying species of the world with high productivity of 20-40 
m^/ha/yr. In view of the above, the necessity of Eucalyptus growth and yield modeling in the 
country has also been elaborated in this chapter. 

Chapter-2 entitled ‘Materials and Methods ’ describes the complete details of the study site 
(rainfall and soil characteristics), experimental description (systems, clones, spacing etc.), 
growth observations (time, frequency and number of observations), harvesting details 
(above and below ground biomass, stem volume and index estimation representing the 
productivity). This chapter also includes the mathematical and statistical analysis details 
(including descriptive statistics viz mean, standard deviation, coefficient of variation, 
skewness, kurtosis etc, probability plots, scatter diagrams, fitting of the functions 
(parameter estimates with confidence interval, asymptotic standard error, value). Linear 
and non-linear models of regression analysis were fitted using various functions on the 
observed data. Residual diagnostics (mean residual, mean absolute residual) and various 
residual graphs pertaining to model validation procedures (viz the plot of residuals against 
their expected values for testing the normality of residuals, the auto correction plot of 
residuals for testing the independence of residual, the plot of residual vs. independent 
variate to confirm that the residual are not continuously being over/underestimated and the 
plot of residual against estimate to confirm that the residual have constant variance) were 
plotted to confirm that the regression pre-requisites are fulfilled and the models developed 
are statistically sound. 

Chapter-3 entitled ‘Collection and compilation of existing models on Eucalyptus spp. in 

India and their critical analysis’ deals with the objective of performing a critical analysis of 

available models specifically on the basis of type of mathematical functions employed. The 

information have been compiled on the basis of various available parameters like species, 

state, measurable parameters like DBH (diameter at breast height 1.37m) and heightrange, 

* 

spacing, age, type of equation, number of trees used for developing model etc. so as (i) to 
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infer what extent of work on eucalyptus modeling has been reported in India (ii) to examine 
model accuracy and limitations of various equations used in eucalyptus modeling like 
linear, allometric, sigmoid etc. In addition to the critical analysis/interpretation of the 
existing work, further possible studies and implementations have also been proposed. A 
total of II equations for biomass and 26 equations for volume estimation could be traced 
from the published literature. The common phenomena associated with most of the linear 
growth models, is the negative estimation of size that is the dependent variate (volume, 
biomass) performed well with the observed range of independent variate viz DBH or ifH, 
whereas the predictive values of size (volume, biomass) comes out to be negative for the 
lower range of independent variate as we start moving below the minimum observed value. 
The problem of extreme overestimation is associated with the allometric functions 
particularly just outside the observed higher range of independent variate, whereas the 
intricacy of constant estimation of size is observed with the S-shaped Junctions like 
Gompertz, Richards, Slodoba etc for the higher range of the observed independent variate 
like DBH or D^H. A review of stem volume and biomass equation for eucalyptus in India is 
presented, with a brief introduction of different mathematical equations employed along 
with a critical discussion on their prediction capabilities and limitations. 


Chapter-4 entitled ‘Modelling and comparison of growth attributes (height and dbh) of 
Eucalyptus tereticomis under different systems, spacing and clones’ deals with the growth 
response of different parameters like tree height, diameter at breast height (DBH) of 
Eucalyptus tereticomis planted in different systems with different spacing [agrisilviculture 
(5x4m, 10x2m, I0x5m, 8x4m and 5x5m), compact block (3x3m and 2.5x2.5m), and 
boundary plantation (2.5m)] up to the age of 52 months evaluated in Bundelkhand region of 
Uttar Pradesh. At the age of 52 months, the Mean Annual Increment (MAI) of growth 
attributes under different systems were of height (3.27m) and DBH (3.41cm) in 
agrisilviculture, height (3.14m) and DBH (3.88cm) in compact block and height (3.62m), 
DBH (4.8cm) in boundary plantation respectively. The growth performance was statistically 
at par in all the three systems during the initial period (upto first three years of age) 
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however significant and distinctly higher growth was recorded under the two agroforestry 
systems (AS, BP) as compared to monocropping system (CB) from the third year onwards. 
This chapter also interprets the effect of spacings, clones on tree growth attributes. The 
results suggests a positive effect of intercropping and wider spacing on height and diameter 
growth in eucalyptus based agrisilviculture system as compared to sole tree plantations in 
compact block plantation. This may be attributed to the fact that due to the crop 
management practices (fertilizers and irrigation given to the crop) in agrisilviculture 
system, the root system of tree component got developed at early stage as compared to 
compact block where there were no such management of crop and accordingly it took more 
time in compact block for proper development of root system and thus the growth picked up 
from second year onwards. Observations were recorded on tree height and DBH (diameter 
at breast height-1.37 m) to develop stand growth model. Model fitting was done using non- 
linear least squares algorithm. Schumacher, Richards and Allometric junction have been 
attempted for fitting height-dbh models on the observed data, since these three functions are 
most preferred one’s to fit height —diameter curves in majority of the published models. All 
the three functions shows comparable value of however, by going through future 
prediction capabilities of tree height with respect to extrapolated DBH values, it is clear 
from the fitted graphs that only allometric function gives reasonably accepted predicted 
results, while other models leads to the problem of constant estimation of predicted height 
outside the observed range of independent variate. Residual diagnostics, through pertinent 
graphs, confirmed the accuracy of model estimation. A comparison of models for different 
systems (compact block, agrisilviculture and boundary plantation) is also discussed in this 
chapter. The developed functions are meant to form part of a system for long term 
forecasting of height-diameter growth curves. 


Chapter-5 entitled ‘Modelling and comparison of yield attributes (biomass and volume) of 
Eucalyptus tereticomis under two contrasting systems [agrisilviculture vs. compact block] ’ 
deals with the harvesting of Eucalyptus trees for development of models for growth and 
yield. A total of one hundred and twenty eight trees of Eucalyptus tereticomis (sixty four 
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from compact block and sixty four from agrisihiculture plantation respectively) "were 
harvested to develop regression equation for predicting standing tree’s biomass (above 
ground, below ground, total biomass in kg/tree and biomass index in tons/ha/yr), volume 
(m^/tree) and volume index (m^/ha/year) with respect to tree diameter at breast height 
(DBH). A comparison has also been made, by developing separate prediction equations 
under the two systems, namely agrisilviculture and compact block plantation. DBH of 
harvested trees ranged from 0.95 to 13.69 cm and 0.96-19.43 cm, volume from 0.00034 to 
0.1 16 m^ /tree and 0.00033 to 0.2219 m^/tree, index from 1.08 to 53.08 m^/ha/yr and 0.334 to 
27.82m^/ha/yr, total biomass from 0.319 to 106.81 kg/tree and 0.44 to 184.34 kg/tree, above 
ground biomass from 0.30 to 70.05 kg/tree and from 0.266 to 152.82kg/tree, belowground 
biomass from 0. 1 09 to 36. 76 kg/tree and from 0.113 to 46.52 kg/tree and height from 2.41 to 
18.30 and 2.34 to 18.35 m under compact block and agrisilviculture plantation respectively. 
Component wise separate models were also developed for above ground biomass of the tree 
viz, leaf, branch and bole biomass. Tree wise volume and biomass values were observed to 
be on higher side in agrisilviculture plantation as compared to compact block, and this may 
be attributed to the fact that the nutrients and irrigation given to the crop in agrisilviculture 
system indirectly also helped the tree component and thus resulting in better growth 
attributes. Non-linear functions were fitted to the observed data and thoroughly evaluated 
through residual diagnostic which revealed that the variant of the allometric equations. 
Total biomass = exp (6.21+77.56/(DBH)^-^^ (R^=0.92) and 
Volume= exp (-0.70+73.75/(DBH)^'^^ (^=0.96) under compact block and 
Total biomass= exp (6.1 5+78.97/(DBHy-^> (R^ =0.96) and 

Volume= exp (-0.50+87.1 2/(DBH) (R?=0.95) and under Agrisilviculture were 

statistically valid for the observed range of DBH (l-19cm), i. e. for early stage of tree 
growth under compact block and agrisilviculture plantations respectively for semi-arid 
conditions. These equations, though complex in nature and not as simple as the allometric 
ones, leads to logical and reasonable predictions. These equations may be used for 
estimation of plant biomass and volume from DBH measurements in permanent study sites 
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where destructive sampling is not possible, in addition to their utility for the policymakers, 
planners and research managers. 



Chapter-6 entitled ‘Development of consolidated yield models (volume and biomass) of 
Eucalyptus spp. at national level by employing primary, secondary and simulated data ’ 
deals with the development and validation of generalized equations at the country level. In 
fact, most published biomass/volume equations were developed using trees sampled from 
specific study sites or from sites that represent small regions only. As a result use of existing 
volume/biomass equation with forest inventory data at large spatial scale is sometime 
unreliable because the equations of previous studies may be site specific often disorganized 
and at times in consistant. Furthermore, unless an equation was developed exclusively for 
the species and study region of interest under conditions typical for the study site, it is 
impossible to know which equations to choose for a particular species and site. 
Accordingly, with the objective of developing site independent equations, we compiled all 
the published biomass/volume equation in India at one platform. These site specific 
equations were used to generate five equally space points (dbh points) and then the volume 
were computed. These points covered the range of dbh values specific to each regression 
and were used to generate a new simulated data set. This model ensures that each site- 
specific allometric equation was weighted equally in the formulation of the generalized 
equation. The simulated data set forAGB covered psedo harvested data from 4 states (Uttar 
Pradesh, Haryana, Karnataka, Tamil Nadu) and the volume-simulated data set encomposed 
harvested observations from 6 states (Uttar pradesh, Tamil nadu, Uttranchal, Rajasthan 
and Punjab). The simulated generalized equations are as follows: 

Y~0.00093*(jQ^'^^ for volume-dbh relationship 
Y^ 0.74*(X)^-^^ for AGB-dbh relationship 

These equations have been validated using geophysically and statistically independent data 
sets. The validity of the mathematical and physical assumptions used in developing a model 
and estimating the coefficients is less open to question if the model gives accurate 
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predictions of new data. In effect, the collection of new data provides an overall check on 
the entire model construction process. 


Chapter-7 entitled ‘Development of Tree Growth Modelling Website in Indian Perspective 
and its launching on World Wide Web ’ has been designed with the basic objective of 
providing a platform where the user can have the basic concepts used in development of 
tree growth models at one end and the comprehensive collection of already developed 
model for a specific tree species at other end so as to use them directly for prediction 
purposes. At the initial phase of development of this site, we have concentrated on 
Eucalyptus species modelling. This website has been hosted at the server of lASRI (Indian 
Agricultural Statistics Research Institute), New Delhi and can be accessed at the URL: 
http://mirror.iasri.res.in/net/tsmi/index.htm. This site is a unique collection of published 
biomass and volume equations, on Eucalyptus off course in Indian context only. The general 
portion of this site describes definition of modelling by various authors and principles of 
modeling; types of modelling describing about various types of models like descriptive 
models, prescriptive models, predictive models and simulation models; model development 
process depicting various steps of model development process like model estimation and 
model validation; tree growth modelling techniques describes essential steps of modelling 
tree growth like sampling of trees, measurement of height, diameter measurement, volume 
estimation, height measurement, estimation of biomass components etc.; models relevant to 
agroforestry provides the working outlines of the models relevant to agroforestry worldwide 
and finally have added tree growth modelling publication displaying a compilation of 
publication on tree growth modeling along with full length papers available in pdf format 
for downloading. This website has been designed to help for organizing, managing and 
sharing research information on Eucalyptus modelling in India. HTML, DHTML have been 
used for creating the dynamic pages of the database at the front end and the source codes 
written in HTML for designing the home page of Tree Growth Modeling Web site. The 
unified database has been equipped with interactive web pages to browse information on 
Eucalyptus introduction, modelling introduction, types of modelling, model development 

Ph.]). Thesis on ‘^Documentation, Analysis and Modellii^ of Eucalyptus... at National Level” by NIGHAT JABEEN (2008) 



process, detailed information on biomass and volume models in India related to Eucalyptus 
species. 

As the concluding remarks, in this thesis attempts have been made to transform the complex 
biological tree growth process into meaningful mathematical Junctions whose parameters 
have biological interpretation along with mathematical Justification. These developed 
models are proposed to be used for interpolative as well as extrapolative predictions for the 
biomass, volume and index of standing tree (with out destructively harvesting the tree) by 
simply measuring a single and simply measurable variable namely the diameter at breast 
height. The models developed in this study will immediately provide us the values of 
biomass, volume and index of standing Eucalyptus trees by simply measuring the dbh values 
only. These models have been thoroughly validated for their prediction capabilities using 
mathematical and statistical tools. 
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Objective~l: To aMeet, compUe and collate all existing modds on Eucafyptus spp m India at one 

platform in tite form of a database. 

Achievements: A total of eleven biomass and twenty-six volume equations of 
eucalyptus spp.was compiled all over Irulia. The information have been 
compiled on the basis of various available parameters like species, state, 
measurable parameters like DBH (diameter at breast height I37m} and 
height range, spacing, age, type of equation, number of trees used for 
developing model etc. so as to infer what extent of work on eucalyptus 
modeling has been reported in India. In addition to the critical 
analysis/interpretation of the existing work, further possible studies and 
implementations have also been proposed. The rekdts are detailed in 
Chapter 3. 

Obiective-2: To perform a critical analysis of existing models. 

Achievements: Critical analysis of existing models of eucalyptus spp was ‘*1 1 “[ [7” 

done on the basis of type of equation used viz linear, sigmoid, allometric, Ibdanet 

parabolic etc, their problems were identified and the solution to these § 2 

Imitations have been proposed A review of stem volume and biomass | 

equation for eucalyptus in India is presented, with a brief introduction of ^ — 

different mathematical equations employed along with a critical discussion f 

on their prediction capabilities and limitations. The restdts are presented in -2 - 7 ^ ’ — 

Chapter 3. / 
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Objective-3: To compare the growth and yield in different systems/re^ns. 
Achievements: Growth attributes (height-dbh): Height-dbh growth was 
evaluated in three agroforestry systems viz agrisilviculture, compact block 
and boundary plantation and height-dbh models were developed and 
thoroughly validated. A comparison of models for different systems (compact 
block, agrisilviculture and boundary plantation) is also discussed in this 
chapter. The developed functions are meant to form part of a system for long 
term forecasting of height-diameter growth curves. The results are outlined 
in Chapter 4. 


Achievements: Yield attributes (biomass and volume): yield attributes were 
evaluated in two contrasting systems (agrisilviculture and compact block). 
Biomass and volume equations were also developed and thoroughly 
validated. A total of one hundred and twenty eight trees of Eucalyptus 
tereticomis (sixty four from compact block and sixty four from 
agrisilviculture plantation respectively) were harvested to develop ' 
regression equation for predicting standing tree’s biomass (above ground, 
below ground, total biomass in kg/tree and biomass index in tons/ha/yr), 
volume (m^/tree) and index (m^/ha/year) with respect to tree diameter at 
breast height The results are revealed in Chapter 5. 
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Objective-^: T 0 iewehp c^m^lidated m 0 dels at bmader lewi -^acMt^e^mal #r state/m'Oimm'M lewd} If 

employing primary, secondary and simulated data. 

AcMevemenis: Volume simtdation: Consolidated model of volume of ziir — , , — 

eucalyptus spp was developed at national level covering five stotes .* 

(Uttar pradesh, Tamil nadu, Uttranchal, Rajasthan and Punjab) and * 
thoroughly validated with independent datasets. The simulated - I / 

generalized vohme is: Y=0.00093*(}Q^'^ ^ •»*“’ 

The results have been portrayed in Chapter 6. ^ ^ j * 


Achievements: Aboveground biomass simulation:The simulated data set 
for AGB covered psedo harvested data from 4 states (Uttar Pradesh, 
Haryana, Karnataka, Tamil NadujThe resulting above ground biomass 
equation is Y- 0.74*(fCf'^^ . This equation has been validated using 
geophysically and statistically independent data sets. The validity of the 
mathematical and physical assumptions used in developing a model and 
estimating the coefficients is less open to question if the model gives j 
accurate predictions of new data. In effect, the collection of new data 
provides an overall check on. the entire model construction process. The 
results have been presented in Chapter 6. 


Objective-5: To develop and launch a “tree growth modelling" website in Indian pei 
Achievements: The “Tree Growth Modelling website” been designed 
with the basic objective of providing a platform where the user can have 
the basic concepts used in development of tree growth models at one end 
and the comprehensive collection of already developed model for a 
specific tree species at other end so as to use them directly for prediction 
purposes. At the initial phase of development of this site, we have 
concentrated on Eucalyptus species modelling. This website has been 
hosted at the server of lASRI, New Delhi and can be accessed at the 
URL: http://mirror.iasri.res.in/net/tem/index.htm. This site is a unique 
collection of published biomass and volume equations, on Eucalyptus, off 
course, in Indian context only. This website is compilation of biomass 
and volume equations along with their abstracts; modelling 
introductions, types, tree growth modelling process etc. The results are 
depicted in Chapter 7. 
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Introduction 


1.1 ^g^’olorestry and Tree Growth Modedng; An overview 

Agro^ , „,e of 4e .os. ™de.y p.o«oed s.s.e. io .he world, in which 

woody pere».als (tees, shmbs etc.) are grown in associahon with herbaoeops plants 
^rops, p^) and/or lives.och in spahal and tempond anangenrent m ag^foreshy 
fliere a. bofl. ecological and economic interactions between fte woody and non-woody 
^nrponents of tte system. Ihns. tire term •agroforestiy' encompasses a divmse se. of 

I ariu ^ 

m almost all temestital ecosystems and provide a range of product and services to rmal 
mtd ^an people. Witi. clearance of namral foreste for agricultnm and oti.er types of 
development flte economic and ecological benefite of tiees could be susteined by 
mtegranng i, inte agriculhtraUy productive landscapes tinongh agroforestiy. 

Agroforestiy combines agricultitre and fomstiy technolo^es to create more integrated 

til! T’ 

utipormnt agroforestiy practices include: agrisilviculmre, agrihorticultitre. sivipasture 
a ey croppmg ete. According to the Association for Temperate Agroforestiy / 
grofores^ practices am mtennw combinations of tmes with crops and/or 
avestock titet ntvolve mrensfve matmgement of 4e i„,eracHons betivemt tire 
component as an integrated agro ecosystem”. To be called apoforestiy, a land use 
prac ce must satisfy all of fltese criteria: Intentional: combination of trees, crops 
an^or an^ are intentionaUy designed and managed as a whole unit, ratirer titan as 
m ivi ua e ements that may occur in close proximity but are controlled separately 
m ensrve: agroforestiy practices am intensively managed to maintein their productive 
and protective functions; these practices often involve annual opemtions such as 
cultivation and fertilization; Interactive: agroforesTry management seeks to actively 
mampulate tire brological and physical interactions between tme, crop and animal 
componenrs. The goal is .o enhance the production of more titan one harvestable 
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component at a time, while also providing conservation benefits; Integrated: Tlie tree, 
crop and/or livestock components are structurally and functionally combined into a 
single, integrated management unit. Integration may be horizontal or vertical, and above 
or below ground. Such interaction utilizes more of the productive capacity of the land 
and helps in balancing economic production with resource conservation. 


Productive area available in the county is around 300 million ha of which 76.5 m illion 
ha is recorded forest area. Actual forest cover is 63.34 million ha of which 26.13 
million ha is degraded. About 20 million ha is covered under private tree planting 
(agroforestry, farm foresty, social foresty and other plantations). In order to meet the 
increasing demand of our fast growing population, we would require to boost the 
production of food grain and fuel wood for human and green and dry fodder for 
livestock to the tune of 250,350 and 2085 million ton, respectively, besides 75 million 
m^ of timber. Thus, to meet the gap between demand and supply of fodder, fuel and 
timber, the emphasis is to put on the integration of tree with agricultural crop, which is 
precisely agroforesty. Approximately 1.2 billion people, making 20% of the worlds 
population depend directly on agroforesty products and services in rural and urban 
areas of developing countries (Leakey and Sanchez, 1997) [109]. 

Almost, the entire 76 million ha recorded forest area is owned and managed by the State 
Governments in India. India occupying just 2.5 per cent of the land area of the planet 
earth has to support nearly 15 per cent of the world’s human population and equally 
large, but mostly unproductive livestock. Therefore, the forests are under intense biotic 
pressures leading to degradation of forest resources. In India, 32 million ha of forest 
area has less than 40 per cent crown cover density. Forests in India have vey low 65 
m^/ha growing stock compared to the world average of 110 m^/ha. Likewise, mean 
aimual increment is vey low at 0.5 m^ha'^yr'^ compared to the world average of 2.1 
mVy‘ (Lai, 2008) [104]. 
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Substantial improvement in productivity of forest resources and promotion of large 
scale farm-forestry plantations particularly of fast growing tree species are most 
essential to meet the basic needs for fuel-wood, timber and wood products and 
conservation of the natural resources on sustainable basis. The present forest resources 
with growing stock about 6, 414 million cubic mt, average volume of 61 .72 cubic meter 
/ha and total annual increment of 87.62 million cubic meter cannot support the much 
needed demands of fuel wood, timber and raw material to industries and domestic me 
(Chauhan et ah, 2008)[39]. Therefore, more attention is required towards production 
forestry, which involves raising plantations of fast growing species under short rotation 
and intensive management. Short rotation forestry is a silvicultural practice by which 
high density sustainable plantations of fast growing tree species with annual woody 
production of at least 10 tonnes or 25 m^ per hectare within a rotation period of less 
than 30 years are established (Verma, 2008)[219]. In India, industrial plantations of 
short rotation species like Eucalyptus, Poplars, Casuarinas have been raised for various 
purpose. The area under agroforestry with clonal Eucalyptus plantations is expanding 
fast in many states as farmers are expanding production of timber and pulpwood to 
meet tiie growing demand. Clonal Eucalyptus based agroforestry plantations are 
prevalent in many states like Andhra Pradesh, Orissa, Maharashtra, Punjab, Haryana, 
Uttrakhand and Uttar Pradesh. Many intercrops like oats, barseem, cowpeas, ginger, 
wheat, blackgram, turmeric and millets have been successfully raised during the Rabi 
and Kharif season. Large-scale eucalyptus plantations have been raised on forest and 
farm lands, community lands, road/rail/canal strips. Intensively managed plantations 
based on genetically improved clonal planting stock and supported with soxmd 
silvicultural practices can achieve productivity levels of 20-40 m^a'^yr'^ even on the 
marginal lands. Significant improvements in quality of produce and reductions in per 
unit production costs have also been possible with the use of true to type, uniform and 
genetically improved clonal planting stock. 
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clones Uke C3, C6. C7, CIO and C27 developed a. Indian Tobacco 
onK,mbon (ITC) Bhadrachalam. Andhna Pradesh formed flre basis of initial clonal 
P^hons sh.ce 1992,Average pmducdvit, of comnrerCl clones is around 20 toS 

------ 

athacnve .a ■ forestry an economicaUy 

2007)1103] ”” '™“ti'">dl crops. (Lai. 


mtroduced m Indra as early as the year 1790, at nandi hills, Mysore. In 1 843 it was 

— yh.h„ducedinnil^hiOsh.s„uthlndia.Howeverre^arpl^^^^^^^ 

E,T r. (Ta^anian snow gum), 

^ ^ch^p/eunu (Mountain gum, 71c, /olio (Red flowering gum, gC 

f “ blue gmn, g„„n,,. (cider gum, yoWnll (VeL ^,1 

^^/on (Whrte rmnbaflr, £ purvr/o/iu, £ pauciflora subsp. n^HophUa (slow 

gem , £ pemnrunu (Spinning gum). £. ^ero^on (Red ironbark, £ umigera (Um 

^), fte ^di provenance of B.cafypp. reren'comlr. popularly known as mysore 

^ most pronnsmg. Synonyms include. Wypms remneom. var, pruiniflora 

^ y Carnemn; Eucalypp^ Naudin; Eucalypp^ populifoUa Desf. ■ 

suMaU. Schauer ; Euca,ypp. „.he/lum (Gaerm.) Domin nom. illeg : 

uc^s ™wum van pm^yiara Blakely. Them have also been nnmeroul 

^spaces and varieties published, but Ore only one ftat remain current is 

.. mbsp. me., ana and the antonym K reredcomlr subsp, temneom. 

E* specres is of great economic irnpormnee. as it form flre basis of mamhbox' 

P md paper producuon indushies and valued for oil, gum and timber. Eucalyptus oil 

used n. nredrcine to treat many sreknesses mch as infections, colds, flu, sore throats. 

bmnchrtrs. pneumonm. aching, stiffiress. neuralgia and even some skin infections. The 

timber rs clanned to be the stiong. dense and most durable in tire world and used for the 
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building of ships, railway ties, piers anH i. , 

H- ^ ^ telegraph poles. Eucalyptus oil is readilv 

distiUed ftom the leaves and used for cleaning deodoriain, = a ■ 
in fr.r.,r I te small quantities 

m food supplement, especially cough drops and decongesmnt £ JT 

plantationsaccountedfor4150001iaoutofthetotol a ■ 

fxrnwo tr, u • u* d? reported in India by 1974 it 

grows to a height of from 20 to 50 mePes nn^ ^kk ^ a 

^ ^ ^Pte 2 metres. In eeneral > 

lerencomu: ts toly fine of pests and except teimites attach ’ 

insecticide not used while planting. In India the m t • ^ 

canlter caused hy Cn„ic,„m sn Jtcnlo ^ 

species for raising plantations in denuded and bririrrs ^Talt T'" 
value natuml cmps. This species was fust mised on 1 
stste in 1952 and experiments were subseouentlv 

plantations were undertaken in Punjab and H ®atea Extensive 

demands of fuel wood, smaU timber and pulpwo^,”"^'' >«raforestry to meet the 

Botanical Classification 


Kingdom: Plantae 
Division: Magnoliophyta 
Class: Magnoliopsida 

Order: Myrtales 

Family: Myrtaceae 
Genus: Eucalyptus 
Species: E. tereticomis 
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1.2 Present state of knowledge 

Eucalyptus has been widely accepted by the fanners not only in pure plantations but 
also along with crops in agroforestry because of its short rotation and hi^ economic 
return. Eucalyptus has also been recommended as one of the six tree species by the 
Planning Commission through its report on “Task Force on Greening India” (Planniiig 
Commission, 2001)[143]. 

The yields that have been obtained in existing short rotation fast growing plantations 
(often with exotic species) have been many times greater than tiiose in natural forests. 
There are many documented reports of tire growth and biomass of tropical plantation 
species in relation to adjacent native forest but there are relatively fewer reports of 
biomass response to experimental manipulation of silvicultural inputs. Various aspects 
of Eucalyptus have been studied and published in recent years viz Optimum spacing 
(Rawat, 1990[157], Naugraiya and Puri, 2001[131]); Antibacterial activity and gastro 
protective properties (Nair et ai, 2008[129]; Sin^ and Sharma 2005[181]; CeimeUi et 
al., 2008[33]; Adeniyi et ah, 2006[1]); Diversity of termite (Varma and Saran, 
2007[215]); Sustaining productivity over multiple rotations (Sankaran et al, 2007[170]); 
Height-growth and site index equations (Tewari et al, 2007[206]); Utilisation in 
makrng base paper of playing cards (Dutt et al, 2007[50]); Nutrient dynamics through 
htterfaU (Rana etal, 2007[153]; Turner and Lambert, 2008[214]; Guo etal., 2006[66]); 
Carbon sequestration potential (Sudha et al, 2007[193]; Stape et al, 2008[190]); 
Bioremediation (Maiti, 2007[1 13]); Transpirational water loss (Kallarackal and Somen, 
2008) [89]; Fungicidal and pesticidal properties (Tripathi et al, 2008[21 3]; Richard, 
2006[164]); Electricity generation (Amatayakul and Azar, 2003 [9]) and many more. 

In recent years there has been an increase in the number on studies of forest biomass 
(Munoz er al, 2008[126]; Mani, and Parthasarathy, 2007[115]; Harmand et al, 
2004[70] ; GiU et al, 2001 [62]; Rana et al, 2002[152] ; Mishra and Nayak, 2000[122]). 
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Thakur and Singh, 2005 [209], studied growth (diameter, height and survival 
percentage) and biomass estimates of Eucalyptus tereticomis plantations of different 
provenances and reported that highest aboveground biomass was observed in Raipur 
provenance (21.69 t/ha). Hi^er proportion of aboveground biomass was allocated to 
bole (60.4-63.3%) followed by branch (12.01-14.46%) and foliage (4.14-4.81%). 

Turner and Lambert, (2008)[214] quantified above ground tree biomass of Euclayptus 
grandis and E. pilularis at north coast New South Wales and reported that maximum 
aboveground tree biomass of E. grandis and E. pilularis was 14.62t/ha/yr at 27 yrs and 
8.18 t/ha/yr at 33 yrs, respectively. The estimated maximum current annual 
accumulation of tree biomass was 16.4 t/ha/yr at 5 yrs and 15.7t/ha/yr at 7 yrs of age for 
E. grandis and E. pilularis respectively. 

Datta and Singh, 200.7[43] reported a timber volume of 19.68 m^/ha/yr for Eucalyptus 
hybrid at a spacing of 3mx3m at 1 6 yrs of age with MAI of 1 .77x1 0'W/ha/yr. 

Tandon et al, (1993)[198], evolved biomass equations for 4, 6, 8 and 10 year old 
plantations of Eucalyptus hybrid from Haryana state. Biomass estimates were made 
using the stratified tree technique: a total of 17 sample trees were harvested, and 
measurements were made on height, dbh along with aboveground biomass (stems, bark, 
leaves and twigs). Various models were used for the biomass estimates, and those using 
diameter at breast height (DBH) as independent variable were selected. Biomass = 
0.1353357*(DBH)^‘‘^^^*‘^ has been reported to be the best function. Regressions were 
derived for total biomass, and for each component separately. Density and total 
aboveground standing biomass of the plantations were 1805, 1116, 844 and 955 
stems/ha and 20.1, 34.6, 88.7 and 137.7 t/ha, for 4, 6, 8 and 10 year old plantations 
respectively, mean DBH range for plantations of different ages were 3.3-20.29 cm. 
Validation of model was not attempted. Biomass components also increased with age, 
however, while stem biomass increased with increasing diameter and age from 54 to 
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78%, the per cent allocation for biomass of the other components decreased 
proportionately. 

Gill and Ajit (2005)[61], compared the growth characters of the twelve multipurpose 
tree species including Eucalyptus tereticomis under semiarid conditions and reported 
that it maintained its superiority over remaining trees for growth characters like plant 
height and diameter. Moreover, the growth attributes were better in agrisilviculture 
system as compared to compact block. 

Patil et ai, (2000)[141], studied the growth and yield of eucalyptus at the age of 17 
years and reported that the height was 20.4m and timber volume was 65.6 m^ / ha that 
led to the income of Rs 67,883/ha/yr. 

Height-diameter curves have been developed frequently in forestry studies by several 
authors (Sharma and Parton, 2007[175]; Lootens et al, 2007[111]; Castedo et al, 
2006[30]; Reed et al, 2003[161]; Mehtatalo, 2005[121]; Inoue and Yoshida, 2004[78]; 
Sochacki et al, 2007[188]). Bachpai et al, (2005) [15] developed multiple regression 
models for prediction of height using dbh and age for Bambusa tulda Roxb. plus 
clumps in regions of Assam and Meghalaya. Tewari and Gadov, (2003)[203] developed 
a hei^t-diameter relationship for Prosopis cineraria for the hot arid areas of India. 
Tewari et al, (2002) [207], have attempted height-age and diameter-age models for 
irrigated plantations of Rajasthan and Soarces and Tome, (2002)[185], have developed 
height-diameter equation for Eucalyptus in first rotation. 

Ajit et al, (2006)[6], fitted non linear functions for Above Ground Biomass (AGB) to 
the observed data and thoroughly evaluated the model using through residual diagnostic 
which revealed that the allometric equations: AGB=0.264* (DBH)^’^^^ and AGB 
=0.421* (DBH)^'®®^ were statistically valid for a wide range of DBH (1-40) under 
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energy and boundary plantations, respectively in semi arid conditions for estimation of 
above ground tree biomass. 

Harvesting of fast growing tree species at young age is now a days increasing due to 
pressing demands of wood in paper and wood industry for manufacturing plywood, 
fiberboards, match sticks, packing cases etc. This has led to change in wood utilization 
pattern. As a result, growth and yield models at early phase of tree development are 
viewed as an important research and management tool to estimate the values of volume, 
Diomass and height of standing growing stock without destmctive sampling. Models, 
hough very few as discussed above, have been developed for volume / biomass 
irediction, but these are region / location specific and moreover utilize numerous input 
ariables Thus most of the models are complex in nature and consequently cannot be 
hrectly used. 

ree growth models have ability to predict the future yields and provide chance to 
cplore silvicultural options. They provide an efficient way to prepare resource 
•recasts to explore management options and silvicultural alternatives. Models are 
squently used to visualize cause-effect of relationships to explain and anticipate the 
havior of systems. Models may be descriptive or predictive. A mathematical model is 
;e a descriptive model, uses mathematical relationships that are more concise and less 
ibiguous than natural language. A tree-stand growth model is an abstraction of the 
utral dynamics and generally refers to interlinking of series of mathematical 
rations related to component growth as per an assumed flow chart of the system. 

model is considered as simplified representation of reality. Many complex 
Tactions and results are depicted with simplicity to reach a decision. It is expected to 
a good representation of those factors that influence production and management, 
dellmg provides tools to predict the consequences of an action that would be 
irwise expensive, difficult or destmctive but is essential for launching a programme. 
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diapter-l: Introductioii 


It also helps in understanding the needs of inputs to optimize the desired output and 
system sustainability. Models can be broadly classified into two categories namely 
linear and non-linear. Linear ones, though simple to fit and easy to use, do not in 
general properly mimic the real tree growth. On the other hand, non-linear ones more 
closely approximate the real biological tree growth process. Accordingly, the non-linear 
functions have been addressed, attempted and analyzed for this work. 

Modelling has been used in natural sciences since centuries in one way or the other. 
However, system modelling as a discipline in itself established in the current century 
only. Models have been built in the physical, biological and social sciences. The 
building blocks for these models have been taken from calculus, algebra, statistics, 
geometry, probability and nearly every other field within mathematics. Various 
researchers have defined modelling in their own way, some prominent and worth 
mentioning are: 

Mclone and Andersons, 1976[120] described modelling as the representation of our so 
called ‘real-world’ in mathematical terms so that we may gain a more precise 
understanding of its significant properties and which may hopefully allow some form of 
prediction of future events. 

Hall, 1963 [68] said “The goal of modelling is to imderstand reality mathematically”. 

Meyer, 1985[119] defined model to be an object or concept that is used to represent 
something else. It is in reality scaled down and converted to a form that we can 
comprehend. He further adds that a mathematical model is a model expressing 
mathematical concepts, having constants, variables, functions, equations, inequalities, 
etc. 
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Agroforestry modelling is a complex process that essentially involves at least two 
elements viz tree and crop, where tree is the long term and crop is the short term 
component. Agroforestry models, usually, predict tree growth and yield and the 
reduction/enhancement in crop yield as influenced by the tree component in temporal 
and spatial sequence. Models providing the most accurate prediction of hei^t include 
basal area, trees/ha, dominant stand height, and diameter at breast hei^t (Sharma and 
Parton, 2007)[175], whereas the models for prediction of biomass components or 
volume utilizes dbh(D) or height (H) or a combination say D^H as independent variate. 

A model is a way to represent, in mathematical terms, the relationship tihat exists in the 
real world, we can use models to describe what might happen in future based on certain 
conditions and assumptions. Intuitively we use models every day although we probably 
do not think about it consciously. For example, you know how much time you must 
allow to get home from your work place. If you choose a mitii bus it usually requires 
less time than a regular bus. If you wanted you could have described that relationship 
with a mathematical model. Models are handy tools because they parallel the way we 
human beings think. Once you describe relationships that exists in the real world you 
can then extend the model you developed to predict what will happen if those 
relationships does not remain the same. Thus, modelling is done to aid the 
conceptualization and measurement of complex systems and some times to predict the 
consequences of an action that would be expensive, difficult or destmctive otherwise to 
do with the real system. Thus, you can easily do “what-if ’ modeling. Just change your 
model assumptions and conditions, repeat the calculations, and observe the predicted 
outcomes. 

Although review of literature, described in the preceding sections reveals that although 
model have been attempted for eucalyptus based systems in India, but they are few in 
number and of limited predictive use since (1) they generally utilize data pertaining to a 
single harvesting only (2) they are strictly location specific (3) validation of the models 
have not been attempted in most of the published papers. Accordingly to overcome 
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these limitations and to fulfill these research gaps, this study was initiated with the 
following objectives 

Objectives 

The overall objective of this study was to develop growth models that will hold good at 
zonal / national level for prediction of growth and biomass of Eucalyptus tereticomis 
plantations. Specific objectives were: 

1. To collect, compile and collate all existing models on Eucalyptus spp in India at 
one platform in the form of a database. 

2. To perform a critical analysis of existing models. 

3 . To compare the growtih and yield of E. tereticomis in different systems/regions. 

4. To develop consolidated models of Eucalyptus at broader level 0ocal/regional or 
state/national level) by employing primary, secondary and simulated data. 

5. To develop and launch a “tree growth modelling” website in Indian perspective. 
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Chapter-2 

Materials and Methods 


MateMs aaidMefiiMis 


Materials and Methods 

2.1 Study site 

The experiment was conducted at the Central Research Farm of National Research 
Centre for Agroforestry, Jhansi which is located at an elevation of 300 m above sea 
level and is situated between 24°ir N latitude and 78°17’ E longitude having tropical 
semi arid climate with mean annual rainfall of 900 mm. More than 75% of the rainfall 
is received during monsoon season (last week of June to first week of September). 
Mean maximum temperature ranges from 47.4° C (June) to 23.5° C (December) and 
mean minimum temperature from 27.2° C (June) to 4.1° C (December). The soil of the 
experimental area is sandy loam with low organic carbon, nitrogen, phosphorous and 
medium to high in potassium. 

2.2 Experiment description 

The study was initiated in August, 2003 with the objective of studying the structure, 
biomass, productivity and modelling the growth of Eucalyptus based systems. The 
experiment was laid with three systems of Eucalyptus tereticomis namely 
agrisilviculture (AS), block plantation (CB) and boundary plantation (BP). Tree 
spacmgs considered in agrisilviculture are 5x4 m, 10x2 m, 10x5 m, 8x4 m and 5x5 m; 
compact block are 3x3 m and 2.5x2.5 m and boundary plantation is 2.5 m. Wheat and 
Black gram were grown in the interspaces during Rabi and Kharif seasons respectively 
under boundary plantations and agrisilviculture systems. Four Eucalyptus clones 
namely C-3, C-6, C-7 and C-10 were obtained from ITC, Bhadrachalam, Andhra 
Pradesh and planted in field during August, 2003. 

2.3 Growth observation 

The growth observations were recorded after every six months in the third year (2006) 
and fourth year (2007). Observed parameters include total height, DBH (Diameter at 
breast height), CD (Collar diameter) and canopy of the each individual tree. A total of 
882 trees from the complete experiment were measured for the third year. However, 
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from fourth year only 862 trees were marked from the whole experiment for further 
growth observation. 

Time schedule for growth observation in third year. (2006). 

1. June-2006 

2. December-2006 

Time schedule for growth observation in fourth year. (2007). 

1 . June-2007 

2. December-2007 

2.4 Volume-biomass studies 

Eight trees (two from each clone) were selected from compact block plantation of 
spacing 2.5x2.5m and other eight were selected from agrisilvicultuxe (new plantation at 
the spacing of 5x4m). 

Time schedule for harvest data 

1. September-2006 

2. March-2007 

3. September-2007 

4. March-2008 

2.4.1 Above Ground Biomass components 

Trees were felled 2-3 cm above the ground. The diameter at breast height (DBH) was 
recorded by measuring vernier calipers. The total height of felled tree was measured by 
tape in meter. The felled tree was divided into one meter long pieces starting from the 
toe to the top of the tree. The leaves were separated and the fresh weights were recorded 
in kg/tree for leaves, branches and main bole for every one-meter piece. Fresh samples, 
were collected for leaf, branch and bole, oven dried at 70°c till constant weight and 
finally dry matter percentages were computed. Fresh weights were multiplied by dry 
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matter percentages to obtain the oven dry biomass of leaf, branch and bole separately. 

The three components were added to yield the total above ground dry biomass (AGB). 

2.4.2 Below Ground Biomass components 

For below ground observations, root system were harvested carefully and biomass of 
main root, primary root, secondary root and fine root were measured separately. Fresh 
samples for these components were collected and oven dried at 70°c till constant wei^t 
and finally dry matter percentage were computed. Fresh weights were multiplied by dry 
matter percentages to obtain the oven dry biomass of all components separately. The 
total components were added to get the total below ground dry biomass (BGB). 

2.4.3 Over-dry Weight 

Over-dry weight was obtained by putting sample in an oven at 70° Celsius for 48 hours 
until there were no more changes in weight. Over-dry weight to fresh weight ratio was 
used to convert fresh weight to oven-dry wei^t. 

> Foliage Sample: All foliage (tree wise) were separated and weighed for fresh weight. A 
representative sample for each harvested tree, were collected and sealed in a plastic bag. 
This sample was taken to the research laboratory to obtain oven-dry wei^t. These 
provided an oven-dry weight to fresh wei^t ratio. This ratio was used to obtain oven- 
dry weight of foliage fresh wei^t for harvested trees. 

> Branch Sample: All branches of a tree were weighed to obtain fresh weight. A 
representative sample for each harvested tree, were collected and sealed in a plastic bag. 
This was taken to a research laboratory to obtain over-dry weight, which provided an 
oven-dry weight to fresh weight ratio. This was used to obtain oven-dry weight of 
branch fresh weight of all harvested trees. 
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> Bole Sample". For dry biomass estimation purposes, felled trees were divided into one- 
meter logs and these meter logs were weighed to obtain fresh weight for each tree. 
Samples were collected from the top, middle, and bottom sections for each tree. Finally 
the dry matter percentage was multiplied by fresh weight to obtain oven dry weights for 
each tree. 

2.4.4 Tree Volume Estimation 

Tree volumes can be estimated in a number of ways. Most method estimate different 
portions of the tree volume using different formulas. Fig. 1 below indicates the most 
commonly used formula for each portion of a tree. 



Fig.l : Volume formulas suggested by portion of the stem. 

Each of these sections of the stem can be estimated using one of four basic equations; 
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Cylinder - The cylinder formula is the simplest of all the formulas. It requires only the 
area of one end and the length of the cylinder. 


Cy{irei«’ 



Vaa Abh 


Fig. 2: The formula for the area of a cylinder 





Paraboloid 


m^Ms and Meiiads 



Paratoloid 



2 


’"=''6(A^t4A 


''-|(a.*a,) 

V- h(A.) 


Smalfeuife formula 
Hufaer'a formula 


A)-h (Newtons Formula) 


Fig. 4: The formula 


for a paraboloid 


and the frustum of a paraboloid. 


• Neiloid 


Neiloid 



Fig. 5: The formula 


for a neiloid and the 


frustum of a neiloid. 
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For estimating the volume, observation of lower 
for every one-meter piece of bole of tree. 


upper diameter was recorded 


In present study we used Newton’s Fomiula, (by „ ^ 
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or V - + Artrl + nrl)*h 


= ^(3.14*(-^)' + 4*3.14*(^)^ + 3.14*(^)")*A 


1 314 314 314 

_(_^*^2 ^ 4 h =_^*^2 ^ [-.jj, in meter, 'dbh' in cm] 


^ Adi ^ dl)*h 


1 d d d 

- - * 0.784 * ((^)^ + 4(^)^ + )*h [ Now 'h' and 'dbh' , both are in mts] 


1 df + Adi + dl 

- * 0.784 * (— = * h 

6 100*100 ’ 


[lm= 100cm] 


: — * *(df + Ad^ + d^')*h 

6 100*100 ^ ^ 


- * 0.0000784 * (d; + Adi + dl)*h 
6 


= I * 0.0000784 * idf + Adi + dl) 

6 

2 

= — * 0.0000784 * {df + Adi + df) [for our case h = 2m, i.e log of 2mts] 
0 

= ^ * 0.0000784 * {df + Adi + ) 

V = 0.0000261 * {df + Adi + df ) [where d„ = 

Where V-volnme of the log (to be calculated) 
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where height in meter, diameters (i.e. d’s) in centimeter and resulting volume in m^ 
in m^ (Cubic meter) 



Fig.6: Sectional piece of the main bole for computation of volume 

For calculating the volume of upper end of tree the formula for parboiled were used 
which is Fo/M?Me= (0.00003927) 

Volume of Paraboloid = ^*A*l 
2 

= [Area=a:r"] 

1 2 

= ~(0.784*^ */) [/ inmts, dincms] 

= |(0.784*(:j^)^*Z) [lm=100cm] 

= I (0.0000784 

Where d is the lower diameter and 1 is the length of remaining part of the tree 
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2.5 Index/Productivity 

The formula for index calculation 

Index = (volume of one tree) *(number of tree in one hectare) 
since lhactare = 10,000 m^ 

Number of tree m one hectare = 1 0,000/(spacing m^) 


Volume Index = 


volume{m^ / tree) * 10,000 , 

: — r- ; — -(unif.m ! ha I yr) 

spacing{mxni)* ageiyrs) 


j . . biomasslkg! tree)*\OfiQQ 

Biomass productivity = — — 

spacingimxm) * age(j;rj) * 1000(to convert into tons) 

Biomass(kgl tree)*\0 . ,, , x 

= : — — ; — - {unit: tons / ha / yr) 

spacmg{mxm) * age{yrs) 

2.6 Statistical Analysis 

The statistical analysis of data were performed using SYSTAT software version 
12.01.04 (Wilkinson and Coward, 2007 [225]). The statistical analysis includes 
descriptive statistics (viz mean, standard deviation, coefficient of variation, skewness, 
kurtosis etc), probability plots, scatter diagrams, fitting of the functions (parameter 
estimates with confidence interval, asymptotic standard error, value). Linear and 
non-linear models of regression analysis were used for fitting various functions on the 
observed data. Residual diagnostics (mean residual, mean absolute residual) and various 
residual graphs pertaining validation procedures (viz the plot of residuals against their 
expected values for testing the normality of residuals, die auto correction plot of 
residuals for testing the independence of residual, the plot of residual vs. independent 
variate to confirm that the residual are not continuously being over/under estimated and 
the plot of residual against estimate to confirm that the residual have constant variance) 
were prepared using various graphical options in SYSTAT 12 (systat Inc@2007, CA, 
USA). 
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2.7 


I*h.a Thesis 


Field experiment layout 


The layout was as below. 
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Collection and compilation of existing models on Eucalyptus spp. in 
India and their critical analysis. 

3.1 Introduction 

The estimation of stem volume and tree biomass is needed for both sustainable planning 
of forest resources and for studies on the energy and nutrient flows in ecosystem. 
Determination of biomass and growth in vegetation is also an essential component for 
the estimation of carbon sequestration. Biomass is often estimated using allometry. For 
scientific purposes standing biomass is a fundamental state variable in several 
ecological and eco physiological models. Traditionally the determination of 
aboveground tree biomass has been related to ensure sustainable planning of forest 
resources. The forester’s apply different methods to obtain such estimations. 
Undoubtedly, tide most commonly used mathematical model for biomass studies takes 
the form of the power function Y = a*(X)^ where ‘a’ is scaling coefficient and ‘b’ is the. 
power function, Y is total biomass or volume and X is the Diameter at breast height (D) 
or total height (H) or D^H, as the case may be. In most cases, variability of Y is largely 
explained by the variability of X. The standard method to obtain estimates for the 
coefficients a and b is by the least square regression of the raw or log-transformed data 
for X and Y measured from destructive sampling of trees that represents the diameter 
range within the stand under investigation, This is laborious and time consuming 
process and it would be difficult to implement it as such at the national level until and 
unless we make a thorough study of all available equations. Accordingly this work was 
initiated for an in-depth study at the broader level. 

Looking to the complexity of problems for obtaining the data of tree harvesting for 
different parts of the country, existing equations need to be compiled and evaluated to 
facilitate identification of the gaps in the coverage of the equations and more, so that the 
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ziirT r - 

ones, as weU as to validate the pn^oess based models. 

ZZ ” '"=■ “0 

Whole tree h' tree species growing and researched in India. The 

whole tree biomass as well a«! * • , 

consideted fot the shtdy, l,e Imlrlr ™ 

Of these equations This dat h ^ S^ide to the ongmal publications 

informsH v effective exploitation of existing 

Le of 11 . nteasntesproposed win also be useful in 

case of growth and yteld models for other species in addition m Eucal^ms. 

3.2 Data 

the available ' f “flctto compile 

India 1 “ ““ '-“lyP'-a in 

specilsH^d"” information like state, 

. measured parametem, uni, of panmtete, DBH and heigh, range, spacing age 

used equation, equation type, number of tiees used for estimating model, em. During 

r:, 

spec.es or sties, for aU cases included in tins meta - database. 

Hitt 

^ways tire dtarnemr a. breast hei^t (D), the tiee heiglti (H) or combination of 

Ic!' 1 “’’“ve ground 
(Bruce and Schumacher, ,950[24]. fhese two variables (D and H) are tite most 
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variableMik!*^ “-iependent vanables, but uqaatiom with several other independent 
eouations ^ 

oontparisorn “ 

ZZ an? r“ 

to biomass and volume separately. 

3.3 Results 

3.3.1 Biomass Equations 

^ found a total number of 1 1 biomass equations from different journals covering five 
tes Of euoalypms (^. rer«. £ HytH,. K ,U,mus. E. ^an,is and 

Most Z T 
very l e T" 

very equattons wem available for tte root biomass. In om compilation, we found 

_d but for other five equations D^H were also measured. Por five equations we 

seven ” f spacing of tree plantation, 

equanons wem allomefrio, dnee were linear and one equafion was 

molecular, ^me equations mentioned age of tees and for two equations we could 

not &d number of tees used for model estimation. A brief description of foese 
equations is as follows. 


J. « al (2006X6], fitted non linear functions for AGB to foe obsewed data an. 

otoughly evaluated foe model using forougb residual diagnostic which mvealed foa 

e aUometac equations: Biomass=0.264*(DBH)^ >’s Biomass=0.421‘(DBH)‘ * 

were ^ttcally valid for a wide range of DBH (140) under energy and boundai, 

P an a OM respectively in semi arid conditions of Uttar Pradefo for estimation of abov. 
ground plant biomass. 
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Tandon et al, (1988)[197], established linear biomass equation for Kerala state at 3, 5, 

7, 9 and 1 1 year old trees of Eucalyptus grandis. Regressions were calculated between 
independent (DBH and height.) variables. Total Biomass= -3.2911 + 0.0207*D^H 
produced the best results (and required the least data), which were significant at 1% 
probability level. The total biomass ranged between 64.797 kg/tree at 3 year age to 
300.564 kg/tree at the age of 1 1 years. Validation of model has not been performed. 

Sharma, (1978)[176], developed linear equation for Uttar Pradesh of 8 to 10 year old 
Eucalyptus hybrid and reported the equation as Biomass= -6.10+0.02*(D^H). A linear 
equation was also evolved for representing the relation between biomass and volume of 
tree. 

Tandon et al, (1993) [198], fitted many biomass equations for Haryana state in 4, 6, 8 
and 10 year old plantation of Eucalyptus hybrid. Biomass estimates were made using 
the stratified tree technique; a total of 17 sample trees were harvested, and 
measurements were made on above-ground biomass viz, stem, bark, leaves, twigs, 
branches). Various equations were used for the biomass estimates, and those using 
diameter at breast height as independent variable were selected. The reported equation 
was Biomass = 0.1353357*(DBH)^ '^^^®^ . Density of the plantations were 1805, 1116, 
844 and 955 stems/ha at 4, 6, 8 and 10 years of age, respectively, and total aboveground 
standing biomass was 20.1, 34.6, 88.7 and 137.7 t/ha, respectively. DBH range was 3.3- 
20.29 cm. Validation of model has not been performed. Biomass components also 
increased with age, but while stem biomass increased with increasing diameter and age 
firom 54 to 78%, that of the other biomass components decreased accordmgly. 

Negi and Sharma, (1987)[132], has reported biomass equation for 6 plantations of E. 
globulus (5-16 year old) and 14 plantations of Eucalyptus hybrid [E. tereticomisj (5-20 
yr old) m Tamil Nadu. In each plantation, 4-5 sample trees with diameter and height 
near the average values for each diameter class were felled and values obtained for 
above-ground biomass. Root biomass for E. tereticomis was estimated separately from 
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data of Uttar Pradesh. Total biomass per ha was calculated by summing the calculated 

biomass per diameter class. Regression equations were given for relations between D^H 

and dry weight. Reported equation was Log (total biomass)= -1.2335+0.09041*Log 
(D^H). 

Rai, (1 984)(a), developed allometric regression equations viz 
Log (Biomass)- -0.5181+0.8195*Log (D^H) and Biomass = -325. 69+3 .97*(DBH) for 
Karnataka state. Equations were presented for the relation between biomass of various 
components of 12 MPTS {Eucalyptus longana one of them) and various logarithmic 
functions of D and D^H. Log-log relationships gave the best fit. The best estimator of 
bole biomass was D^H, and for other components it was D. 

Other Biomass Equations 

Kushalappa, (1984)[99], presented regression equation W/D^= 0.1 139 + 1.562*(1/D^) + 
0.00367 (H/D) + 0.01949*H (D-DBH, W-biomass and H- height) for above-ground 
dry weight, in reUtion to DBH and total height, based on samples from a 12-year old 

plantation of Eucalyptus tereHcomis with spacing 1.2 x 1.2 m for BCamataka state, the 
DBH range was 2-26cm. 

Rai, (1984) (b)[150], developed linear biomass regression equation with respect to 
vanable D^H for bole, leaf, branch and root biomass separately, no equation was found 
for total aboveground biomass, the study was conducted for Eucalyptus Longana. 

Dhawan, (1990)[45], developed correlation between leaf area and leaf dry mass and to 
leaf number, based on data collected from seedlings of Eucalyptus hybrid [E. 
tereticomisj over the course of a year (starting at the 14-16 leaf stage). The equations 
were Y = 333.258 +92.277*X and Y = 166.626 + 20.808 which were highly significant 
(r = 0.99 and 0.98, respectively). A multiple regression equation relating leaf number to 

leaf area and leaf dry weight was also given. 
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Dogra,and Upadhyay, (2005)[48], studied for Eucalyptus hybrid in Punjab, measured 
data from 133 pennanent experimental plots with age and top heights ranging from 3.9 
to 1 1 years and 6 to 26.2 m respectively and developed various equations that depicted 
the relation between height and age of tree. Reported equation was Height = 3 1 .5249 - 

61.0317/A™withR2 = 0 . 47 . 


Pal and Raturi, (1991)[136], evaluated the growth parameters and frequency 
distribution of trees of different categories in an energy plantation of Eucalyptus hybrid 
at the age of three years, and observed that the growth of individual trees in the 
plantation was very heterogeneous. Variable height was 550-990 cm and girth at breast 
height (g.b.h.) 5-35 cm. The larger trees (average height 990 cm and g.b.h. 29.0 cm) 
were much fewer in number (IP/o of the total) than the medium sized and smaller trees. 

It is to be mentioned here that all the equations discussed above are not directly 
comparable (at a starch in one go) in the sense that they used different dependent and 
independent variates, varying ages and tree spacings etc. however they can be compared 
by classifying them on the basis of the type of function used. Accordingly attempts 
have been made in this study to include maximum possible similar information and 
details about each published model on biomass and its components (Table-1). 

The DBH and hei^t range in the biomass models varied from 2-40 cm and 10-30 m 
respectively with Rvalues greater than 0.95 for majority of models. However for four 
models moderate values of R^ were reported ranging from 0.6 to 0.8 . 

The three models at Sr. No. 1,2 and 5 used power equation of the form Y=aX'* with ‘a’ 
values rangmg from 0.14 to 0.42 and ‘b’ values from 1.96 to 2.14. Ter- Mikaelian and 
Korzukhm, 1997[200] have also reported almost similar ‘b’ values ranging from 1-3.3 
for 65 North American tree species. All these three corresponds to above ground 
biomass from two states namely Uttar Pmdesh and Haryana and both these areas fall 
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under the semi and climate accordingly taking the average of the values of ‘a’ and ‘b 
parameter may hold good for prediction of AGB under semi arid conditions. 

Table-1: Published biomass equations in India at a glance 


s. 

No 

1 

Species 

State 

Variables 

used 

Unit of variables 

Observed range 

Biomass 

DBH 

Height 

DBH 

Height 

X 

Eucalyptus 

tereticomis 

Uttar 

Pradesh 

DBH, 

Biomass 

Kg 

cm 

m 

1-40 

10-30 

Z 

'i 

Eucalyptus 

tereticomis 

Uttar 

Pradesh 

DBH, 

Biomass 

Kg 

cm 

m 

WO 

10-30 

j 

A 

Eucalyptus 

tereticomis 

Uttar 

Pradesh 

DBH, 

Biomass 

Kg 

cm 

m 

1-40 

10-30 


Eucalyptus 

tereticomis 

Uttar 

Pradesh 

DBH, 

Biomass 

Kg 

cm 

m 

1-40 

10-30 


Eucalyptus 

hybrid 

Haryana 

DBH, 

Biomass 

Kg 

cm 

m 

3 . 3 - 
21 (1 


O 

Eucalyptus 

globulus 

Eucalyptus 

Hybrid 

Tamilnadu 

1)BH, D^H 
Biomass 

Kg 

cm 

m 


- 

7 

Q 

Eucalyptus 

globulus 

Eucalyptus 

hybrid 

Tamilnadu 

DBH, D^H, 
Biomass 

Kg 

cm 

m 


- 

o 

Q 

Eucalyptus 

grandis 

Kerela 

"OBH, D^H, 
Biomass 

Kg 

cm 

m 

- 

- 

y 

1 A 

Eucalyptus 

hybrid 

Uttar 

Pradesh 

“dbh, D^H, 

Biomass 

Kg 

cm 

m 

- 

- 

10 

t t 

Eucalyptus 

longana 

Karnataka 

dbh, d^h. 

Biomass 

Kg 

cm 

m 

“ 

- 

11 

Eucalyptus 

longana 

Karnataka 

DBH, 

Height, 

Biomass 

Kg 

cm 

m 


- 


Table continued on next page 
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S. No, 

Age (yrs) 

Spacing 

Cm) 

Equation 

for 

No of trees 

Reference 

1 

11 

2x2 

AG 

86 

Ajit et aL, 2006[6] 

2 

11 

2x2 

AG 

86 

Ajit et aL, 2006[61 

3 

11 

3 

AG 

86 

Ajitefa;.,.2006r61 

4 

11 

3 

AG 

86 

Ajit et aL, 2006[6] 

5 

4, 6, 8, 10 

- 

AG 

17 

Tandoneta/., 1993ri981 

6 

5-16 & 

5-20 

- 

AG+BG 

28&53 

N^ and Shanna, 

1987[132] 

7 

5-16 & 

5-20 

- 

AG 

28&53 

Negi and Shanna, 
1987[1321 

8 

3-11 

- 

AG+BG 

- 

Tandoneta/., 1988ri971 

9 

8-10 

- 

AG 

- 

Shanna, 1978[1761 

la 


lOxlOm 

2x2m 

AG 

433{bole) 

221(branch.) 

20iaeaf) 

Rai, 1984 (a)[149] 

11 


lOxlOm 

2x2m 

AG 

433(bole) 

221(brandi) 

2010eaf) 

Rai, 1984 (a)[149] 


AG- Above ground, BG- Below ground Table continued below 


S.No 

Equation 

Functions 

Parameter values | 


a 

b 


1 

Y = a*(DBH)‘' 

Allomertic 

0.26 

2.17 

0.70 

2 

Y = a*bBH)“ 

Allomertic 

0.42 

1.96 

0.82 

3 

rY = a + b*(DBH) 

Linear 

-149.08 

15.99 

0.80 

4 

Y = a + b*{DBH) 

Linear 

-376.31 

28.05 

0.60 

5 

Y = a*(DBH)‘’ 

Allomertic 

0.14 

2.41 

0.97 

6 

Log(Y)=a+b*Log(D^H) 

Allometric 

1.23 

0.09 

0.97 

7 

Log(Y)=a+b*Log(iyH) 

Allometric 

-1.04 

0.84 

0.98 

8 

Y=a + b*(D"H) 

Linear 

-3.29 

0.02 

0.99 

9 

Y=a + b*(D'‘H) 

Linear 

-6.12 

253.24 

0.99 

10 

Log(Y)=a+b*Log(D^H) 

Allometric 

-0.52 

0.81 

0.99 

11 

Y=a+b*(DBH) 

Linear 

-635.41 

55.78 

0.95 
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3.3.2 Volume equatioiis 

A total of twenty-six volume equations could be traced from the published literature 
(Table-2), covermg five different Eucalyptus species (E. Hybrid. E. Globulus, E. 
Grandis, E. Camendulansis, E. Piularis) for six states (Tamilnadu, Karnataka, 
Uttrakhand, Madhya Pradesh, Punjab, and Rajasthan) of India. 

In sixteen of die compiled volume equations, the independent variable were diameter at 
breast height, total height and volume over bark, and six of them use various 
ematical combinations of those like and D^H. A total of four equation measured 

volume over bark as well as volume under bark. In seven equations we could not find 
the DBH and height range. Mostly all ihe compiled equations for volume were linear, 
except seven equations those were parabolic and one equation was allometric. Two 
p p mentioned age of trees and only one paper mentioned spacing. Dogra and 
Sharma, (2003)[47], does not reported the number of trees used in model estimation. A 
brief description of these equations is presented below. 

Jain et al, (1993a)[80], prepared volume tables (diameter at breast height (D) versus 
volume over and under bark) and the equations were, (Volume) over back = 0.00183 + 
0.00003309*D"H and (Volume) „„derb^k= - 0.121986 + 0.0344338*D foxE. globulus in 
the Nilgiri HiUs of TamU Nadu, based on combined data from eight sample plots in 
plantations and coppice forests, and using best fit regression equations and fiimival 
mdex to calculate volume as a function of D^H and the DBH. The DBH and height 

range was 6-30 cm and 8-32 m respectively. 

Smgh et al, (1995)[179], prepared volume tables using data from roadside plantations 
of Eucalyptus hybrid [E. tereticomis] in Uttar Pradesh. Data, were coUected from 85 

trees at 3 X 2m spacing, and different regression models were tried, such as 
(Volume)overba±~ 0.001188 + 0.00003108*D^H and 
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(VolumeWr bark = -0.090281 + 0.00080164*dI In each model, diameter at breast 
height (D) was taken in cm, height (H) in meters and total wood volume overbark 
(VOB) and volume under bark (VUB) in cubic meters. At the age of five years the 
range of height and diameter was 8-26 m and 9-39 cm respectively. 

( 93b)[81], evolved volume tables, for total volume over- and under-bark, 

together with the regression equations, 

(Volume) over bark ~ " 0-00458+0.000036 1*D^H 

(Volume) u:^erbark= - 0.0461 76+0.0008743 8 *d' containing tree height as a variable, 
used m their calculation, based on data from 40 trees in Dehradun, Uttranchal state. The 
observed DBH and height range was 6-30 cm 8-26 m respectively. 

al., (2001)[205], proposed volume equations for Eucalyptus camaldulensis at 
the Indira Gandhi Nahar Pariyojana Area in Rajasthan. Six different models were fitted 
to each data set to compute the single and double entry volume equations. Compiled 
equations for over bark volume with respect to variables D^H and D were 

Tvob = - 0.00308 + 0.0000333*D^H 
Sqrt (Tvob) = - 0. 1 1 887 + 0.03 1 077*D. 

Jam et al., (1991)[79], developed regression equations (Volume)over bark = - 0.00458 + 
0.0000361 *d2h and (Volume) ^er bark = - 0.046176 + 0.00087438*D" representing the 
relations between volume and diameter at breast height. The independent variables 
considered were D^H and D^. They also reported that it is often difficult, and even 
erroneous, to fit a single equation to the whole range of diameters, and in these cases it 
is advisable to divide the diameter range into small segments and derive a model 
(equation) for each separately. This technique was applied to diameter and volume data 
from two types of Eucalyptus plantations in the Punjab: block unirrigated plantations 
and canal side strip plantations. The 'piecewise' models were shown to give more 
accurate results than a single model, and based on these, local volume tables were 
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presented for both plantation types, giving total volumes over and under bark over a 10- 
40 cm DBH range for both types of plantation. 

Pandey and Jain, (1976)[139], presented tables derived from measurements of 233 E. 
grandis trees from plantations of varying density in Kerala and Tamil Nadu, for 
different age class (6, 8, 10, 12 and 14 years) and developed linear equation of the type 
V= a + b* D^H + c* D^H*Log A + d*H for volxime, for the DBH range 5-34cm, and 
height range 6-42m (here A is the basal area). 

Surey et at, (1999)[195], developed volume tables for Eucalyptus Hybrid for Madhya 
pradesh by calculating best fit regression equation using data of 993 sample trees. The 
best fit regression equation was of the form Log (YOB) = a + b*Log(D)+ c*Log(H) and 
Log(VOB) = a + b* D+ c* (D^H)+( D^H) ^ 

Dogra and Sharma, (2003)[47], presented linear volume equation with variable D and 
D^H separately. Equations were of the form Log V=a + b*Log D + c*Log H, 

Log V= a + b*Log D and V=a + b* D^H for Eucalyptus Hybrid in Punjab. A total of 
111 felled trees ranging in diameter from 7.1cm to 55.7cm and height from 10.0-36.4m 
at the age of 9-1 1 years were used for the study. 

Chaturvedi and Khanna, (1982)[36], presented linear and parabolic equations of the 
type V= -0.0015 + 0.24*D^H and Vv = -0.084 + 2.5*DBH fot Eucalyptus hybrid, with 
D H and D as an independent variable. 

Tewari and Kumar, (2003a)[201], proposed different volume functions viz V= -0.00514 
+ 0.000033*D^H, V = 0.000169* for Eucalyptus camaldulensis at the Indira 
Gandhi Nahar Pariyojana Area in Rajasthan, India, and concluded that the combined 
variable equation performed well in both, the fitting and vahdation process. 
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Prasad, (1984)[147], studied on three fertilizer treatments of two varieties of 
Eucalyptus, 'F.R.L-4' and T.R.I.-5', planted at Gamharia Research Centre (Bihar). The 
reported equation was Vv = -0.0868 + 2.8335*DBH. 

Other volume equations 

Pande and Chaturvedi, (1972)[138], developed correlations between volume over bark, 
volume tmder bark with top height, age, mean diameter, top height, on the basis of 284 
sample plots from all over India for Eucalyptus hybrid and Eucalyptus tereticomis. 

Chaturvedi, (1973)[35], developed linear regression equations, relating mean annual 
increment to age, basal area and top height for Eucalyptus hybrid. Data were collected 
from Punjab, Haryana, Uttar Pradesh, Bihar, West Bengal, Madhya Pradesh, 
Maharashtra, Karnataka, Tamilnadu, Andhra Pradesh, Kerala etc. Spacing was varying 
from 1x1 to 4x4m and age varied from 4 years to 20 years. 

Pande, (1974)[137], developed a curvilinear regression equation of stand volume as a 
function of top height and stem numbers per ha. The equation was V = -62.453599 + 
13.488326*X + 1 .972614*X^ for 9-12 year old E. grandis in Kerala and Tamil Nadu. 


Other articles on volume components, but without equations 

Choubey et ah, (1990)[40], determined the sample size for estimating the volume of 
Eucalyptus in an agroforestry situation, where diameter, height and volume were 
measured in 1986 in each of five field boundary plantings of Eucalyptus hybrid in 
Ambala District of Haryana. The plantations were of four types; single sided (1 row 
along 1 field boundary, 6 yr old); double sided (rows along 2 sides of the field 
boundary, 6-8 yr old); triple sided (rows along 3 sides of the field boundary, 2 
plantations of 6-8 and 12 yr old); and a square plantations (along all 4 sides of the field 
boundary, 6 yr old). It has been reported that althou^ variability is not uniform for a 
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giv^ populaaon toriMios fr„„ p,a.e to place and for varying ages, toe 
icient of vanation remained fairly stable for population of differing locations and 

ages. 


Smgh et al, (1995)[179]. prepared volume tables using data from roadside plantations 
of Eucalyptus hybrid [B. tereticomisj in Uttar Pradesh. Single tree data were coBected 
from 85 trees planted at a spacing 3x2m. Different regression models were attempted. 
In aU repotted models, diameter at breast height (DBH) is taken in cm. height (H) in 
meters and total wood volume over baric (VOB) and volume under baric (VUB) in cubic 

s. At the age of 5 years, the average height and diameter was 11.1m and 11.4cm 
respectively. 


Prajapati, (1993)[144], reported that a local volume equation of the form V=d*b has 
been found to be the best fit for a Eucalyptus pole crop (planted in 1983 and felled in 
1993) in Dharwad, Karnataka. The coefficient of correlation was 0.9, and the 
coefficient of determmation was 0.8, showing that the 80% variation in volume was 
explamed by the variation in diameter at breast height. 

Agnihotn er al., (1989)[2], studied on 20-25 year old E. alba. E. resinifera. E 
tereticomis, E. hybrid [E. tereticomisj. E. crebra, E. citriodora, E melanophloia, E. 
gomphocephala, E rostrata and E robusta for Chandigarh planted at a spacing of 
2x2m and reported that the volume was best estimated by a regression equation 

mcluding basal girth and height as variables, although an equation using diameter at 
breast height alone was also satisfactory. 

Kulkami and Gore, (1988)[94], studied on a 20 year old Eucalyptus hybrid [E. 
tereticomis] trees in plantation in West Nasik Divsion, Maharashtra, and derived 
multiple regression equations for volume prediction based on girth, height, girth + 
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height, and girth + height + girth 
respectively). 


at 1, 2, 3, 4 and 5 m ht. (gl, g2, g3, g4 and g5 


As 3.11 theses volume eematinno x i 

1 ■ -1 . j, compared directly, therefore maximum 

T bl ;> iT ^ ” ““nation and details of twenty six equations have been compiled in 

Table-2. The reported DBH range in most of fte volume models is 5-30em where as in 

ew eases it was up to 45cm. The reported height range was 8-32 meters, surprising 

only m one ca^ it ™s r^orted to attain a height of 42 meters. The values for aU the 

mode s are exhemely satisfectoty ranging flom 96o/. to 99«/. employing drereby that die 

1 3 d 0 iTi **” Observed values of harvested data. Equations 

1 f; 

ra an . ttiakhand, Rajasthan used similar function i.e of the form volume- a+b 

4 fr oqnations, 1,3, 6.9, and 21 although developed 

or tfferent states having varying climatic conditions, different age and spacing) were 

a most of comparahle parameters i.e. -a' and -b' values are more or less alike infeiring 
diemby drat a common model of volume (over-bark by avemging a and b values) 

enve outofthemmay hold good at national level. 

However critical comparison of the biomass and volume models can be made on the 
bams of functional fom. of dre equation developed viz linear.' allometric/power, 
sigmoid, parahohe etc.) Accordingly an analysis of the feasibiUty, limitations and 
pre iction capabilities of die equation has been made below as per the functional form. 
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TabIe-2: P„bUak.d votame eqaatiom In India at a glance 

ffi Tc— ITT- 


s. 

No 

i 

Species 

State 

V^ariabies used 

Unit dr V 

ariahles 


Observed Range 

2 

■t^ucalyptus 

globultis 

Eucctlvnfiiv 

Tamil Nadu 

DBH, Height, Volume 
(OB.UBI 

Volume 

m"* 

DBH 

cm 

Ht 

m 

DBH 

6-30 

Height 

8-32 

3 

^globulus 

lamil Nadu 

DBH, Height, Volume 
(OB, US’) 

—3 

cm 

m 

6-30 

8-32 

4 

hybrid 

Uttar Pradesh 

UBH, Height, Volume 



cm 

m 

9-39 

8-26 

“s 

hybrid 

Uttar Pradesh 

OBH, Height, Volume 


cm 

m 

9-39 

8-26 


hybrid 

Uttar Pradesh 

UBH, Height, Volume 


cm 

m 

9-39 

8-26 

6 

7 

Eucalyptus 

camaldulensis 

Uttranchal 

DBji, Height, Volume 

li? — 

cm 

m 

^30 

8-26 

8 

camaldulensis 

Uttranchal 

UBH, Height, Volume 

li? — 

cm 

m 

6-30 

8-26 


camaldulensis 

Uttranchal 

OBH, Heigjit, Volume 

1? — 

cm 

m 

6-30 

8-26 

9 

10 

Eucalyptus 

camaldulensis 

Kajasthan 

■UBH, Height, Over 
bark volume of timber 
& small wood 

li? — 

cm 

m 


11 

camaldulensis 

Kajasthan 

UBH, Height, Over 
bark volume of timber 

small wood 

— 

cm 

m 

- 

- 

12 

species 

EliCdl'un'hj^ 

Punjab 

tiHH, Height, Volume 

^ 

m 

cm 

m 

7-32 

- 

13 

grandis 

EucdTvTifij c 

Kerala, Tamil 
Nadu 

a^hjh, Height, Volume “ 
^Height, DBH ~ 

— ■ 

cm 

m 

5-34 

6-42 

14 

yj^ius 

pilularis 

i rt 

Madhya 

Pradesh 

_iOB, UB) volume 

li? - 

cm 

m 

- 

- 

15 

j-jui^uiypms 

pilularis 

Madhya 

Pradesh 

Height, DBH 
(OB, UB) volume 

m^ 

cm 

m 

- 

- 

16 

mjucui ypius 
hybrid 

Punjab 

DBH, Height, Volume, 
tfH 


cm 

m 

7.1^3.9 

10.0- 

17 

^uLcuypcus 

hybrid 

Kut^nl'MTt'N # CT 

Punjab 

DBH, Height, Volume, 
D'H 

m 

cm 

m 

7.143.9 

32.8 

10,0- 
‘7'y o 

18 

j^ui^uL yptus 
hybrid 

Punjab 

DBH Height, Volume, 
D'H 

~3 

m 

cm 

m 

7.143.9 

32.0 

10.0- 

19 

^i4^K^U-lypHlS 

hybrid 

/t >1^ /f # CT 

Uttaranchal 

DBH, Height, Volume, 
D^H 

1? — 

cm 

m 

- 

32.8 

20 

^uL^uiypjus 

hybrid 

Uttaranchal 

UHH, Height, Volume, 
D'H 

— 

cm 

m 

- 

- 

21 

^uLUiypius 

species 

iT't # av 

Uttaranchal 

UHH, Height, Volume 

1? — - 

cm 

m 

- 

- 


t^ucaiypius 

camaldulensis 

Rajasthan 

DBH, Height, Over 
bark volume of timber 

^ — 

cm 

m 

“5-52 

- 

22 

Kucctlvntiiv 


& small wood 






< 

camaldulensis 

Rajasthan 

DBH, Hei^t, Over* 
bark volume of timber 

— 

cm 

m 

5-52 

- 

Table coi 

itinued on next p 

>age 

& small wood 
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23 

24 

25 

Jtiucalyptus 

hybrid 

Eucalyptus 

hybrid 

Rurnl^mfii o 



X uiJL|ap * 

t^jab 

Volume 

DBH, Height, Volume 

li? — 

rrf 

cm 

cm 

m 

m 

14^6 

10-16 

- 

26 

Table c 

hybrid 

Eucalyptus 

hybrid 

ontinued below 

Punjab 

Punjab 

UBH, Heigbt, Volume 

DBH, Heigjht, Volume 

"s? — 

cm 

cm 

m 

m 

10-32 

2«2 

- 



(Voluiiie)i)veiharir~3.+ b*D + c*D^ 


(V Plume) o« ri,rir= a + 


Linear 


Linear 


-0.0208 


0.006597 


0.0009 


0.96 


7 

9_ 

10 , 

11 


12 


13 


14 


15 


V= a + b*D^ 


Linear 


^=a + b*D + c*iy H 
V= a + b*l5^ 


Linear 


-0.0045 


0.000036 


Linear 


-0.0461 


0.000874 


0.0234 


^(Twob)=a + b*D 


Linear 


-0.009753 


-0.0030 


0.0011 


Vv-a + b*DBH 
V=a+b*D^H+c*D'‘H*logA+ ' 
d*H 

Log (VOB)= a+b*Log(DBH) 
c*Log(Heieb.tl 


Linear 


0.000033 


-0.1188 


Parabolic 


0.031077 


Linear 


-0.1545 


-0.0377 


0.035370 


0.216865 


0.0392 


0.003 


Linear 


0.2000 


2.203000 n 0.000 


Log (VOB)=a+b*DBH-(- (?D% 


.LogV-a + b*Log D + c*Log H 
LogV=a + b*LogD 


Linear 

Linear 

Linear 


0.0018 


-1.3060 


-0.011086 I 0.3614 


1.921000 


0.99 


0.98 


0.99 


0.98 


0.97 


0.97 
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p ‘son, analysis and critical evaluation of various mathematical 
functions being used in tree growth modeling 

3.4.1 Linear functions 

growth models employ linear regression to describe the change in the size/ 
response variate with respect to explanatory variate. The response variate may represent 
tree volume, biomass etc. and the explanatory variate include age, total height (H), 
diameter at breast height (DBH), or a combination of these variables says D^H, to 
descnbe the change in response variable 00 with respect to explanatory variable (X). 

of linear model is very common in forestry/ agroforestry (Chaturvedi and 
Khanna, 1982[36]; Jam c( al., 1993(a & b)[80 &8I]; Tewari, 1996[202]; Kumar etoi, 
1996[96]; Jam et al, 1998[83]; Chaturvedi and Sood, 1995[37]). BasicaUy, these 
models provide the mathematical budding blocks of simple/complex biological growth. 

Problem identification and solution 

The statistical model employed in linear regression is of the form Y=a + b*X+ e, where 
Y IS the response variate, X the explanatory variate, ‘a’ and ‘b’ are regression 
coefficients and e’ the error term. Statistical procedure to estimate ‘a’ and ‘b’ assume 
error term is normally distributed with mean zero, constant variance and 
mdependent between cases (Drapper and Smith, 1982[49]). The corresponding 
mathematical equivalent of simple linear model is the straight line Y=a + b*X, where 
‘a’ denotes the intercept on Y-axis, and ‘b’ tiie slope of the line (b=tan 0), where 0 is the 
angle which this line makes with the +ve direction of X-axis). The method of least 
squares is used to fit this straight line and to estimate the parameters ‘a’ and ‘b’ in linear 
egression using the given data set. There are numerous graphical and statistical 
methods to support the fitting of the equation and consequently model vaUdation viz the 
plot of residuals against their expected values; the autocorrelation plot of residuals; the 
plot of residuals against the estimate; the plot of residuals against the independent 
vanate (Janssen and Titus, 1995[84]; Reynolds et al, 1988[163]; Mayer and Butler, 
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[ 8], Soares et al, 1995[187]). All of these are based on the comparison of 
observed and estimated values and the diagnosis of residuals. Once these statistical 
assumptions are folfiUed, the linear model may be used for predictive purposes. 

However m true sense the biological growth process is not linear, which is the basic 
assump ion in linear models. Due to this reason, simple linear growth models 
sometimes fail to precisely predict the change in the response variable, particularly for 
he lower range of explanatory variable. For this range, these models sometimes 
estimate misleading values and even vague viz-negative values of response variable 

( g 1). This problem associated with linear tree growth models has been addressed as 
negative estimation’ of tree size. 


Thus, the Imear model reflects a shghtly unrealistic picture of tree growth, which is 
sigmoid in nature. Growth of the tree results from two opposing factors: firstly the 
positive component associated with biotic potential, photosynthesis activity, absorption 
of nutrients, constractive metabolism etc. and secondly the negative component 
representing the restraints imposed by competition, limited resources, stress, respiration 
and aging factors. These two components when combined together results in the 

sigmoid shape of tree growth. 

The apparent reason for ‘negative estimation of size’ with linear models (Singh et al, 
1998[178]; Jain etal, 1993(a &b) [80 & 81]; Tokyo etal, 1994[212]; Raj and Cheema, 
1990[151]; Jametfl/., 1998[83]; Gupta eta/., 1996[67]; JaineiaZ., 1996[82], even with 
satisfactory values is that the data set utilized to fit the regression equation represents 
that range of explanatory variable, belonging to the second phase of growth and which 
is linear in nature. However, when we use this linear equation for interpolation 
/extrapolation particularly for lower range the results of prediction are not encouraging 
and some times vague viz. negative estimation of size. In the last phase (third) they lead 
to o verestimation of size. 
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The apparent reason for high values associated with linear models is that the data set 
usually belongs to the second phase of tree growth and which is linear in nature. 
Howe er, reckoning into the prediction capabilities of these models, the results are not 
encouraging. For the initial (first) phase of tree growth they result in underestimation 
and even negative estimation and for the last phase (third) they lead to overestimation 
of size. It is therefore recommended that linear models, in context of tree growth, 
should be used with caution and in case the fitted simple linear equation results in 
negative value of a and positive value of ‘b’, then the lower bound of explanatory 
vanable should be computed as proposed (here viz X > a/b) to avoid the negative 
estimation of tree size (Table-3). Thus when the phenomenon is inherently undesirable, 
as the negative estimation of size in simple linear models for lower range of explanatory 
variable, to put imphed lower bound on the raw values is most appropriate. 
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Table-3; Different eucalvotus linear 

avoid the negative estimation probleiT ^^ equations along with the proposed lower bound to 


S.No 

Model Description 

Reference 

Lower Boand 
proposed 

1 

AGB- - 3.291 1+ 0.0207*D‘‘H 

R=0.99 

Tandonefa/., 1988[197] 

iyH=158.9 

2 

Volime- - 6.10877 + 0 253*D^H 

R'=0.99 

Shamia, 1978[176] 

D^H= 0.024 

3 

(Volumc)ovi!rbiirk— 0.121986 + 0.0344338*D 
R^=0.98 

(Volume)overtark= - 0.0083+0.0000309*D*H 

R =0.99 

Mnetal. 1993(a)[80] 

D=3.54 

4 

(Voliime)„™b^=-0.09028 1+0.0008016^*5^ 

R=0.96 

Singh era/., 1995[179] 

D^=l 12.64 

5 

(Volume)ovcAark«**0.046 1 76+ 0.0008743 8*D^ 
R = 0.98 

(Volume)oveibaric“-0. 0045 8+ 0 . 00003 6 1 *D^H 
R‘^=0.99 

Jainefa/., 1993(b)[81] 

D'^I=127 

D^=52.80 

6 

Twob - - 0.00308 + 0.0000333 *D'‘h 

R^=0.99 

Squa(TWob) = - 0. 1 1 887 + 0.03 1 077*D 

R =0.99 

Tewari et a/.,2001 [205] 

D^H=92.49 

D=3.825 

7 

Biomass = -635.41+55.78*DBH 

R^=0.95 

Rai, 1984 (b)[150] 

D= 11.39 

8 

Volume^ -0.00514 +0.000033*0^^11 

R^=0.99 

Tewari and Kumar, 2003 
(a)[2011 

D'‘h= 155.75 
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Fig. 1: Function graph depicting problem of negative estimation 



A 

Jl.1^ 



1 

I 

1 

1 

/ 

/ • 



■ 










4 ! i s - ■ I 

0 100 200 300 400 500 


CFH 

Biomass = -3.291 1+0.0207*D^H 
R^=0.999 (Tandoneta/., 1988 
[197]) Fig 1(b) 


PIlD. Thfisis Ott “DocHmentation, Analysis and Modelling of Eucalyptus...at National Level” by NIGHAT JABEEN (2008) 


50 


Chaptcr-3: CoiectlM aad CoaipiJjitfea cxistirag modds on Ewalyptas spp.... 


3.4.2 Concept of Non Linear function: 

Functions other than linear (linear means X to the power one only) may be broadly 
classified as non linear models. Non-linear models are in general more powerful than 
published rion-linear Eucalyptus models can be broadly categorized 
three main groups, namely allometric, parabolic and sigmoid. Accordingly 
nonlmear equations has been classified and presented below into three categories. 

3.4.2.1 AUometric/Power function 

Allometry is used to describe the morphological evaluation of species, and is based on 
the relation between a tree biomass, volume and of any part of the tree component The 
general form of the allometric equation is Y = a*X" where y = measure/process in 
question viz volume/biomass, X is size (usually DBH), b is the allometric exponent 
(which tells about the relationship between x & y), and a = a constant (the allometric 
coefficient). Allometric equation, that compensates for nonlinear functions is: log (Y) = 
log (a) + b [log (X)]. Allometric modeling may be theoretically, physiologically, and 
statistically superior to alternative methods of modeling used in tree growth estimation. 
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Tabk-4: Different eucalyptus allometric model equations 


S.No. 

Model Description 

lifefiprptir’A 

1 

momass - 0.1353357*(DBH)‘‘‘'‘'^“' 

R^ = 0.97 

Tsxidmetcd., 1993[198] 

2 

=-1.2335+0.0Wl*Log,(D^H) 

Log(bioniass)= -1.0376+0.84*LogfD^ID 

R^=0.98 ^ 

Negi and Sharma, 

1987[132] 

3 

= -0.5181 -K).8195*Log(D''H) 

l^i, 1984(a) [149] 

4 

tsioinass = 0.264*(DBH)'^‘‘'^ ' — 

R^ = 0.71 

Biomass = 0.42 1 ’"(DBE')*’*®^ 

R' = 0.82 

Ajit et al, 2006 [6] 

5 

Volume = 0.000 169*(DBH)^'"^ 

R^ = 0.99 

Tewari and Kumar, 2003(a) 
[201] 


Limitations of allometric models 

These allometric equations (Table-4) behave reasonably well within the observed range. 
However they lead to steep overestimation (Fig. 2) just outside the observed range of 
mdependent variable. This overestimation intricacies has been clearly depicted below 
through the extrapolated graphs of some of the pubUshed Eucalyptus aUometric models. 
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Fig. 2: Allometric function graph depicting problem of overestimation 



Bioinass= 0.1353357*(DBH)^‘‘’®^*^ 
= 0.99 (Tandon et al, 1993[198]) 
Fig 2(a) 



[201]) Fig2(b) 
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QiapteiuS: Cdi«cti<^ sud Ciw^ktiiHi @f cxiilliig iWHidte on Eucalyptus spp.... 


3,4.2.2 Sigmoid Function 

In forestry/agroforestry studies, growth functions have been used for many years, 
y o provide a statisticaVmathematical summary or prediction for time-course data 
on growdi of trees. lUe term growth function is frequently used to denote die rate of 
change m tree size that may generally be condensed in a single equation. Tie size may 
p sent the growth attributes of tree viz wood volume, wood biomass, etc. as response 

variate and age, total height (H), diameter at breast height (D) or a combination of these 
say D H as explanatory variable. 

The development of growth functions traces back to 18*^ century: the prominent and 
worth mentiomng are the equation of exponential decay (Gompertz, 1825[63]) 
descnbmg the age distribution in human population and Logistic equation (Verhulst, 
1838[218]) which are probably the most famous equations in ecology. Though the 
Gompertz equation was initially developed for population studies, a century later it was 
applied as a growth model (Winson, 1932[227]). The Gompertz equation was found by 
Causton and Venus, 1981 [32], and many other researchers (Laird et al, 1965[101]; 
Zweifel and Laskar, 1976[236]; Nokea, 1978[134]; ZuUinger et al. 1984[235]) to be 
more appropriate in biological work than others. An outstanding exposition of Logistic 
equation is given by Hutchinson (1978)[77]. The early papers of Richards (1959[165], 
1969[166]) also contributed to the development of this area and gave a thorough 
account of growth functions known as Richards functions. The monomolecular 
equation, also known as the law of diminishing returns in agriculture and economics 
and law of mass action in chemistry, has not been widely used in forestry. Putter 
(1920)[148] and Weber (1891)[222] have used this equation in the study of tree growth. 

Growth functions, frequently encountered in forestry/agroforestry modeling, may be 
classified either on the criterion of number of parameters involved (viz two, three or 
four or higher parameter families) or on the basis of general functional form (viz 
exponential, hyperbolic, power, monomolecular, etc.). The classical functions like 
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Logistic, Gompettz. Richards etc. belong tt, 3-parametet family and the recent 
developments viz Slodoba, 1971[182]; Zeide, 1992[231], appettnins to 4-parametet 
family. Amongst the exponential form the prominent contributors are Ratkowsky, 
1990[1S6], Huang and Titus, 1992[74], Buford, 1986[25]; Mofiat et al, 1991(123], 

while Tang, 1994(199]; Curtis, 1967(42), etc. augmented 4e development of 
hyperbolic class of families. 


Majonty of the equations described above share a peculiar feature that their growth 
are S shaped or Sigmoid (Table-5). Recently, various forms of sigmoid 
functions (simple, reverse, rotated, etc.) have been employed in tree modelling viz stand 
growth (Feng, 1997)[56], volume growth (Shibuya etal, 1997)[177], survival function 
(Mason and Whyte, 1997)[117], DBH distribution (Kunisaki and hnada, 1996)[97], 
fruit development and ripening (Nerd et al, 1998)[133] and disease infection levels 
(Wiensczyk et al, 1997)[224]. This section discusses the feasibility and limitations of 
the S shape growth functions and their suitability in prediction of tree size. The 
apparent flaw of sigmoid functions has been primarily addressed to as ‘Constant 
imation of size. However, one of the possible solutions to overcome the ‘constant 
estimation’ mtricacy is to use the concept of simulation to generate the data outside the 
observed range and then fitting the sigmoid function to the generated d ata 

While fitting the sigmoid growth models to the data (Table-6), we may sometimes 
encounter the problem of ‘constant estimation of tree size’ (Fig. 3) depending upon the 
observed range of explanatory variate. Smaller the observed range of explanatory 
vanate, greater is the severity of problem i.e. more is the error of prediction in 
extrapolated zone. The term constant estimation of size’ means that when we use the 
fitted sigmoid models for extrapolation, particularly for higher range of explanatory 
variable, the predicted values of the response variable comes out to be almost uniform. 
In other words, as the value of explanatory variable increases, the value of the predicted 
response increases, but the pace of increment is almost negligible and it ultimately leads 
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to nearly constant estimation of the predicted response. The apparent reason for 

s ant estimation problem ts the small observed range of explanatory variate. When 

the observed range of explanatory variable is large enough, there is no such problem of 

constant estimation of size with sigmoid functions and due to this reason, for tree 

growth models with height or DBH as response variate and time as explanatory variate, 

we may avoid this problem since it is possible to record such data for a long duration. 

, for tree yield models which necessitates destructive sampling with timber 

volume or biomass as response variate and DBH or height as explanatory variate, one 

would have to wait for years together to harvest trees at different intervals for obtaining 

a large observed range of explanatory variable in systematically planned experiments 

h studies. In usual plantation forestry/agroforestty experiments, the trees are in 

general more or less of the same age group, and that is why in a single harvesting, there 

IS practicaUy no variability in the observed valnes of the explanatory variate and which 

nlnmately leads to ‘constant estimation' problem while fitting sigmoid tactions with 
such data. 
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Table-5: Some prominent 


sigmoid functions with their integral forms 


SNo. 

i 

Reference 

general mtegral functional 
form 

Parameter interpretation 

1 

'j 

Logistic 

Y - a (1 + exp ^''’Y 

a= Asyn^tote 

b= Rate of approach to 

the asyn^tote 

p= No direct 

interpretation 

JL 

-1 

Chapman-Richards 

Y = a (1-exp (-b*x)) 

a= Asymptote 

b= Rate of approach to 

the asyn^tote 

c— Shflnf* naramfitftr 

A 

Gompertz 

Y = aexp(-b‘-'^‘'^-^"*^») 

a= Asymptote 

b= Rate of approach to 

the asyrqptote 

No direct 

interpretation 

4 

Slodoba 

Y=aexp[-b.exp{-cjc''}] 

a= Asymptote 

b= Rate of approach to 

the asyii^tote 
c= Shape parameter 

^ No direct 

interpretation 

5 

Weibull 

Y=a exp [l-b.exp{-c.x“}] 

a= Asymptote 

b= Rate of approach to 

the asymptote 

c= Shape parameter 

d= No direct 

interpretation 
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Tabled: Different eucalyptus sigmoid model equation 


S.No. 

1 

2 

Model BescriDtioti 

V louo(i-e3q) (-0.0073A»‘‘'“^ 

PAI= 14.47(6X0 (-0.0112A'i'> 

PDlAj ~420[1 

Reference 

Rayner, 1992[158] 


Reed era/., 2003[161] 


|^-0.7723+0.3334H -0.00043 6 IHN/B 

Garcia and Ruiz, 2003[58] 


PAI Periodic and mean annual increment 
A = Basal area 

PDIAj - The potcmial diameter at the j* measwemeiit data 

mpi -He cumulative number of degme da^ iomplautiugumil the f measured data 

V/B = Volume —basal area ratio 

H, N and B = State variables (Height, Number of trees and Basal area) 


Fig. 3: Function graph depicting problem of constant estimation of size 


2Dnnh- 



Basal 

Voliime=l 605(l-exp (-0.0073A))‘““ 
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Biological Mterpretafi,, mathematical function parameters 
The sigmoid growth functions have a lot of potentiaUties that reflect their snitabiUty in 
tree growth models. To mention a few: most of the paiameteis of these fimctions 
(Table-5) have biological inteqiretation viz the parameter ‘a’ represents the asymptotic 

value ofthe response variable, ‘b' narrates the rate of approach to the asymptote, ‘c’ is 

the shape parameter. However, the parameters -p’ and ‘q’ have as such no direct 

biological mteipretation. GeneraUy, these functions lead to statistically unbiased and 
efficient models of tree growth 


Advantages 

They mimic, the complex tree growth process, more precisely and more accurately and 
parameters of these models have mathematical description and biological interpretation. 
They genetaUy lead to statisticaUy unbiased and efficient models. 


Disadvantages 

These modeis are not transportable to other data sets for extrapolative predictions since 
the asymptotic value of predicted size is approximately equal to experimentally 
observed maximum value of size. For fliese models theoreticany 

y -C. the asymptotic value as x =o. but practicaUy y ^ even sHghtly outside the 

actual observed range ofx and much before the infinite value of X. ' 


Proposed mathematical/statistical solution to the problem 
It is worth mentioning here that liuear/curvilineai'allomettic models some times may 
lead to unpractical predictions (even negative estimation) for lower range of 
explanatory variable and sigmoid models result in unrealistic predictions (constant 
estimation) for higher range. The proposed approach, a combination of curvilinear and 
sigmoid functions, overcomes the drawback of both the models to a certain extent and 
results in more realistic predictions for lower and higher ranges. The reason being that 
the generated data set now includes the feasible 'asymptotic value' of size, which a tree 
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W*... 

may attain. Regarding the selection of sigmoid function, the choice may depend on the 
relative ease of achieving convergence, the equations mathematical properties, and 
possible biological/physical interpretations of parameters (Fig. 4 ). 

Fig. 4: Graph depicting the mathematical/statistical solution to the problem 



Where 

N - Observed data (i.e. the basic observed data set) 
N+Ni - Intermediate simulated data 
N + Ni + N2- Ultimate simulated data 


Let (Xi, X2) represents the observed data range containing ‘N' points. Simple linear 


model is fitted to the observed data range of explanatory variable, and say the 


resulting 
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equation is Y= a + b*X, where Y is the response variate/size; X is the explanatory 
variate, a, b are constants. This developed linear equation is used to simulate/generate 
additional data points (‘ Ni' in number) for the DBH ranged by Xi to X3. The additional 
generated Ni data points are added to the 'N’ observed data points, thus the 
intermediate-simulated data set contains (N + Ni) point. 

Allometric model is fitted to the intermediate-simulated data set consisting of (N+Ni), 
and say the resulting equation is Y=a*X*’, where Y is the response variate/size: X is the 
explanatory variate: a is multiplication constant and ‘b’ is the shape parameter. The 
developed allometric equation used to simulate/generate additional data points (‘N2' in 
number) for the DBH range Xi to X4. The generated additional ‘N2' data points are 
again added to the intermediate-simulated data points, thus the ultimate-simulated data 
set contains (N + Ni + N2) points. Suitable sigmoid model is than fitted to the ultimate- 
simulated data set. 
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3.4,2.3 Parabolic Equations 

A power function with degree two is known as parabolic function. Parabolic function 
has also been attempted to model volume growth for Eucal 3 ^tus species (Table-7). The 
problem observed with parabolic functions of the type Y = bX^ is again illogical 
estimation intricacies as shown in the graphical representation of the Eucalyptus 
volume equation proposed by Jain et al. 1991 [79] (Fig.5). For the initial range, as the 
DBH increases the biomass/volume first decreases from a certain medium values and 
then it again increases with increasing DBH. In fact, the positive constant term (+a) and 
the positive linear term (+cX) result in minimizing/decreasing the estimation error term 
where as the negative constant term (-a) and the negative linear term (-cX) leads to 
maximizing/increasing increases the estimation/prediction error. This phenomenon has 
been clearly depicted in the concerned graphs of parabolic functions (Fig’s. 6, 6a, 6b, 
6a-l, 6a-2, 6b-l, 6b-2). 


Table-7: Different eucalyptus sigmoid model equations 


S.No. 

Model description 

Reference 

1 

= - 0.1545 + 0.035370*DBH 

V = 0.024 - 0.01 1*DBH + 0.001238*(DBH)^ 
R^=0.97 

1991 [79] 

2 

Vv = - 0.0840 + 2.51PDBH 

V = 0.0070 - 0.4228*(DBH) + 6.31*(DBH)^ 

R^ = 0.99 

Chaturvedi and Khaima, 
1982[36] 

3 

VV = -0.0868 + 2.8335*DBH 

V = 0.0075 - 0.492*(DBH) + 8.02*(DBH)^ 

Prasad, 1984[147] 
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Fig. 5: Function graph depicting problem of constant estimation of size 
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Fig.6a: Effect of inclusion of 
negative constant term (-a) in the 
basic parabolic function Y = bX^ 



Fig.6b: Effect of inclusion of 
positive constant term (+a) in the 
basic parabolic function Y = bX^ 
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Fig. 6a: Equation graph of the parabolic function Y = -a + bX^ 
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Fig. 6b: Equation graph of the parabolic function Y = +a + bX= 
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3-5 ^®™P^rison of published models using the same type of equations 
In this section, we have attempted to have a comparison of similar models. In fact the 
comparison is that of the model parameters for example if the equation is of the type 
Biomass a (DBH) then the parameter ‘a’ and ‘b’ are to be compared for different 
published model. Initially a table of parameters is prepared (Table-8), then these 
equations are plotted in a single frame (Fig.Sa) and finally we have fitted a model (Fig. 
8b) in between the two parameters ‘a’ and ‘b’ to find out how the value of these two 
parameters are interrelated. The schematic representation of the comparison made (in 
the following section) is shown in Fig. 7. 
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3.5.1 Comparison of Biomass models (using same type of eqnation) 

In this section, we have compared biomass models. From the published bterature, two 
types of models for biomass could be traces (i) Allometric (ii) linear. Moreover, some 
of the reported models were based on DBH and others on D^H. Accordingly, they have 

been presented below by specifying the type (allometric or linear) and variable used for 
prediction (DBH orD^H). 

Table-8; Allometric equation Y = a*(DBH) " for total above ground biomass vs. DBH 
(parameter comparison) 


S.No. 

Reference 

Parameters ! 


a 

b 

1 

Tandon et 
a/.,1988ri971 

0.135 

2.42 

0.971 

2 

Ajitcra/., 2006, 

-(PP)[6] 

0.264 

2.17 

0.82 

3 

Ajit et al, 2006, 
-IBP)[6] 

0.421 

i.97 

0.70 

4 

Ritu, 2006 
(AS)[1671 

0.227 

2.15 

0.99 

5 

Ritu, 2006 
(CB)ri671 

0.265 

2.08 

0.97 
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Fig.Sa. Pictorial view of above ground biomass allometric equations vs. DBH (plotting in a single 
frame) 



Big. 8b: Relationship between the parameters of allometric above ground biomass equations vs. 
DBH (fitting of model between the parameters ‘a’ and ‘b’ ) 
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Table-9: AUometric equation Y 
(Parameter comparison) 


a*(D^H) 


for total above ground biomass vs D^H 


S.No 

Reference 

Para 

a 

meters 


1 

Ante/ a/., 2006, 

iEP)[61 

"003 

0.94 

0.79 

z 

'1 

Ajit etal., 2006, 
-®[6] 

13.01 

1.08 

0.71 

J 

A 

Kitu, 2006 

~0.102 

0.82 

0.99 

4 

Ritu, 2006 
(CB)ri671 

"0.226 

0.71 

0.96 

5 

Negi and Shanna, 
1987(132') 

"0.09i ^ ■ 

0.84 

0.98 

O 

Rai, 1984(a)[149] 

0.303 

0.82 

6.99 
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Fig. 9a. Pictorial view of above groirnd biomass allometric equatioES vs D^H (pIottiEg ib a siEgle 



Fig. 9b: Relationship between the parameters of allometric above ground biomass equations vs. 
DBH (fitting of model between the parameters ‘a* and ‘b’) 
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Table-10: Linear equation Y = a + 
comparison) 


b*(DBH) for total above ground biomass vs DBH (parameter 


S.No 

Reference 

Parameters I 

R' 

a 

b 

1 

Ajit et al, 2006 
iEP)[6] 

-149.08 

15.99 

0.79 

2 

Ajit et ah, 2006 
-(BP)[6] 

-376.31 

28.05 

0.71 

3 

Ritu, 2006 
1AS)[1671 

-9.815 

4.132 

0.91 

4 

Ritu, 

2006(CB)ri671 

-6.025 

il78 

0.93 1 

5 

Rai, 1984 
-(a)[1671 

-635.41 

55.78 

0.95 
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Fig. 10a. Pictorial view of above ground biomass linear equations vs. DBH (plotting in a single 
frame) 



Fig. 10b: Relationship between the parameters of linear above ground equations vs. DBH (fitting 
of model between the parameters ‘a’ and ‘b’) 
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Table-11: Linear equation Y = a + b*tfH 
comparison) 


for above ground biomass vs D^H (parameter 


S.No. 

Reference 

Parameters I 


a 

b 

1 

Shaima, 

1978fl761 

-6.108 

6.0253 

0.99 

2 

Tandon et al, 
1993ri981 

-1.233 

0.090 

0.97 

3 

Ajit et al, 2006, 
(EP)r61 

T.976 

0.02 

0.79 

4 

Ajiteffl/., 2006, 
-(BP)[6] 

-3.272 

6.025 

0.71 

5 

Ritu, 2006 
-(AS)[1671 

3.169 

0.026 

0.93 

6 

Ritu, 2006 
m[16T] 

0.86 

6.03 

0.98 
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CIiapter-3: CoiectS<m '(if escifi&g wmMst #« Iiacal^^ptw a^ip.... 


frame) *" ground biomass Unear equations vs tfH (plotting in a single 



Fig. 11b: Relationsliip between the parameters of linear above ground equations vs. D^H (fitting 
of model between the parameters ‘a’ and ‘b’) 
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3.5.2 Comparison of Volume models (using same type of equation) 

In this section, we have compared volume models. From the published literature it is 
inferred that in majority of cases only linear models are reported for volume 
predictions. Moreover, these reported models were based on D^H only, and accordingly 
they have been covered in a single heading. 

While making the comparison of published models using the same types of equations, it 
was observed that for the power equation (Y=a*X’’ there was a clear cut negative 
correlation between these two parameters i.e. as the value of ‘a’ increases, the value of 
‘b’ decreases and vice versa. The same trend was observed in biomass Vs DBH (Fig. 
8b) and biomass Vs D^H (Fig. 9b). Almost similar trend was observed for the linear 
equations models (Y=a+b*X) for Biomass Vs DBH (Fig. 10b) and Biomass Vs D^H 
(Fig. 1 lb). The same trend was apparent for the linear models of volume Vs D^H (Fig. 
12b). 


Table-12: Linear volume function = a + b*(D^H) for volume vs. D^H (parameter comparison) 


S.No. 

Reference 

Parameters 


a 

b 

1 

Jain et ai, 1993 
(a)[801 

0.00183 

0.00003309 

0.99 

2 

Sin^ et al, 

199511791 

0.001188 

0.00003108 

0.99 

3 

Jain et al, 

1993(b)[811 

-0.00458 

0.0000361 

0.99 

4 

Tiwari et al, 
20011205] 

-0.00308 

0.0000333 

0.99 

5 

Dogra and Sbarma, 
2003147] 

0.0405 

0.0000285 


6 

Wtu, 2006 
.(AS)[167] 

0.0015 

0.0000375’ 

0.99 

7 

Ritu, 2006 
(CB)tl67] 

0.0011 

0.0000375’ 

0.97 


Mean 

0.005494 

0.0000339 
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Fig. 12a: Pictorial view of volume linear equations vs. D^H 



Fig. 12b: Relationship between the parameters of volume linear vs. vs. D^H (fitting of model 
between the parameters ‘a’ and ‘b’) 
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3.6 Concluding remarks 

The majonly of the published equations on biomass and stem volume of Eucalyptus 
species can be broader classified into four functional types. Linear models, in general 
suffered from the problem of negative estimation of size (Fig. 1) particularly for the 
lower range of independent variate. It is therefore advisable to put a lower bound on the 

of explanatory variate i.e. the lower limit below which the function should not be 
used for extrapolative predictions. 

on linear category of functions (in diis case allometric /power, parabolic, sigmoid) 
usually do not suffer firom the negative estimation problems. However they are 
associated with other kind of estimation intricacies. In particular, the aUometric/ power 
functions leads to extreme over estimation of size, just outside the observed higher 
range of independent variate (Fig. 3). An amicable solution has been suggested to 
over ome the problem of constant estimation of size by combining the curvilinear and 
sigmoid functions after simulation of data in a feasible manner (Fig. 4). This approach 
of simulation of data after initial stage and then applying suitable sigmoid function for 

fittmg at the second stage will lead to a model that may avoid the negative and constant 
estimation problems. 

The parabolic functions, on the other hand, suffers firom illogical estimation problem 
(Fig. 5) m the sense that for the initial range, as the DBH increases the biomass/ volume 
first decreases up to a certain medium value and then it again increases with increasing 

DBH. Accordingly, it is not advisable to use parabolic functions for extrapolative 
prediction purposes. 

Thus looking to the negative, constant, over and illogical estimation intricacies- 
associated with the extrapolative predictions by linear, sigmoid, allometric and 
parabolic functions respectively. It is recommended that while reporting a specific 
equation for prediction of tree wise volume and biomass, not only the value and 
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residuals diagnostics are important but equally important is the behavior of the 
proposed model outside the observed range of independent variate to avoid the various 
t5q)es of estimation problems that are otherwise generally overlooked. 
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Modelling and comparison of growth attributes (height and dbh) of 
Eucalyptus tereticornis under different systems, spacing and clones. 

4.1 Introduction 

In forestry, mathematical equations form the basis of many prediction studies. Output 
from these equations (height, biomass etc.) based on easily quantifiable metrics such 
as dbh (1.37cm above ground level) have long been employed by foresters and plant 
scientists for use in ecological research (Ter-Mikaelian and Korzukhin, 1997[200], 
Anderson et al, 2006[10], Kajimoto e? al, 2006[88], Delphis and Levia, 2008[44]). 
Such models gain importance as modem scientific management relies heavily ©n well 
defined growth and yield functions that can be used to access the status of stands at 
any point of time. Eucalyptus tereticornis being exotic to India, its cultivation has 
spread to nearly the whole of the Indian subcontinent. Large-scale plantations of E. 
tereticornis were established between 1970 and 1985. Nearly 2% of the cultivable 
land in northwestern India is planted with E. tereticornis. The average productivity of 
E. tereticornis plantations is 10 m^/ha/year on forestlands and 15-10 m^/ha/year on 
farmlands. The average productivity of commercial clones is about 20-25 m^/ha/year 
and many farmers have achieved up to 50 m^/ha/year (Lai, 2007[103]). Infact, for the 
computation of productivity of such plantations, precise measurements of growth 
attributes (height and diameter) is essential. It has been noticed during the course of 
experimentation that it becomes extreihely tedious to accurately measure the height 
of large trees (more than 10 meter), whereas the dbh (diameter at 1.37m height) can 
be easily measured and accordingly height-diameter models provide a handy tool for 
prediction of hei^t on the basis of dbh values. Recently height-diameter curves have 
been frequently developed in forestry studies (Sharma and Parton, 2007[175]; 
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Lootens « 2007[m]; Castedo « ,006p01; Reed . «/., 2003[161]; 

e tata o, 0O5[121]; Inoue and Yoshida. 2004[78]; Sochacki etai. 2007[I88]). 
Present y, no work has been reported on height-diameter models for E.tereticomis in 
xcept a few reports like Tewan et ah, (2007)[206], Tewati et al., (2002)(207] 
and Soaroes and Tome, (2002)[185]. on modelling of eucalyptas nnder pure 
plantations. Although growth performance of £ tereticomis has (hei^t, dbh) has 

!r.roT'“‘'. Kumar and Bangarwa, 

], Prajapati et al., 2006[145]; Tandon, et al., 1993[198]; Sundararaju, et al., 

1995[194]; Aulokh and Sandhu. 1990[14]), but most of these studies have not 

considered the role of spacmg and systems on growth attributes which is iiqiortant in 

agroforestry. In the present study attempt have been made to evaluate the growth 

performance of four clones of E. tereticomis planted at eight spacing in three 

agroforestry systems and to develop height diameter models separately under three 

agioforestry systems, for reasonable predictions of hei^t values without destractive 
sampling. 


4.2 Data 


4.2.1 Study Site and Experimental details 

The experiment was cairied out at the Central Research Farm of National Research 
Centre for AgroforesSy. Jhansi which is located at an elevation of 300 m above sea 
level between 24°11’ N latitude and 78”17- E longitude in tropical semi arid climate 
with mean annual tainM of 900 mm. More than 75% of the rainfeU is received 
during monsoon season (last week of June to first week of September). Mean 
maximum tensierature ranges from did” C (June) to 23.5» C (Deceniber) and mean 
minimum temperature flom 27.2° C (June) to 4.1° C (December). The soil of the 
experimental area is sandy loam with low organic carbon, nilrogen, phosphorous and 
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TTr RT, /I ^ ^ 'Tbe experiment was initiated in August 2003 with four 

ITC Bhadrachalam clones, viz Cj r r r< ji - i- , 

3j 6> ^ 7 , Cio and laid in three systems of Eucalyptus 
lereticontis namely, agrisilvicultee (AS), block plantation (CB) and bonndaty 

pl^tton (BP). Tree spacings considered in agrisflvicnlture were 5x4 m (1), 10x2 m 
ni (3), 8x4 m (4) and 5x5 m (5); compact block are 3x3 m (6) and 2.5x2.5 
m (7) and boundary plantation is 2.5 m (8). Wheat and Black gram were grown in the 

mterspaces dnrmg -RaW and ‘Kharif seasons respectively under boundary 
plantations and agrisilviculture systems. 


4.2.2 Data and observations 

The available data includes the tree wise observations on height in meter taken from 
are ground level m tip of the shoot of die main bole and dbh (diameter at breast 
hei^t) m centimeter taken at a height of 1.37m at the nnun bole) for a period of 52 
months. Durmg the fot year, observations were recorded four times in a year on 
1300 frees distributed over die three agmforestry systems viz agrisilviculture, 
compact block and boundary plantation. Whereas from second year onwards, 

observations were recorded twice a year on 862 selected frees distributed 
proportionally over the three systems. 


4.2.3 Model fitting and validation 

The complete dataset consisted of the free wise observation on height (m) and dbh 
(cm) covering an age span of 6- 52 months. Three such datasets were available, one 
for each of the three systems. Each of the three datasets (AS, CB, BP) were randomly 
split into two parts; one to estimate the model parameters (estimation dataset) and the 
other to vahdate the performance (vaUdation dataset). This split sample approach 
used 80% of the otigiual samples for model development while 20 % of the sample 
was reserved for vahdation (Cannela et aL. 2007[29]). The system wise models, 
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developed using die estimation dataset, were used to compute die predicted hei^t in 
the validation dataset. The residuals (predicted hei^t- observed hei^t) were also 
computed and then tested for their normality and other regression requisites. For 
model validation, linear regression was fitted between the observed and predicted 
height of validation dataset (predicted hei^tr=a+b* observed height); for the perfect 
fit model, the standard value of ‘a’ must be zero and ‘b’ must be one, along with a 
perfect umty R value. However, in real experimentations, the value of ‘a’ should be 
approaching zero and ‘b’ should be approaching one. 

4.3 Statistical analysis 

Systat-12 statistical software (Wilkmson and Coward, 2007[225]) was used for 
computation of descriptive statistics (mean, standard deviation, skewness, kurtosis, 
Shaphiro Wilk’s statistic for normality testing etc), and fitting of different equations 
(estimates of model parameters, asyn^totic standard error of estimate, confidence 
interval, R values etc) and plotting of various graphs pertaining to residual diagnostics 
(probability plot of residuals, auto correlation plots of residuals, plot of residuals 
against their expected values, plot of residual against independent variate). Different 
functions, frequently used to model height-dbh curves viz schumacher, chapman- 
richards and allometric were atten^ted (Table-1) and choice of the best one was based 
on validation criterions as well as the prediction behaviour of the function within and 
outside the observed dbh range. 

Table-1: Functional fonn of various functions used for fitting height-dbh curves 


Name of function 

Functional form 

Parameter interpretation 

Allometric 

Hei^t = a*(DBH)'’ 

‘a’ - Scaling parameter 
‘b’ - Power parameters 

Richards 

”Y=a*[l-esp (-k*x)] 

a = Asymptote 

k = Rate of approach to the asymptote 
m = Shape parameter 

Schiimacher 

Height=l .37+a*exp(-b/DBH) 

a = Asymptote 
b= Shape parameter 
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4.4 Results and Discussions 
4.4.1 Growth attributes 

Growth attributes were statistically at par in different systems and spacings, clonal 
differences were significant. At the age of 22 montiis, maximum values of height and 
dbh were recorded for the clone C-3 in 8x4 spacing and similar trend was observed at 
52 months (Table 2). At the age of 22 months, the MAI (mean annual increment) of 
height was 3.96 m in agrisilviculture, 4.1 m in contact plantation and 3.98 in 
boundary plantation. But with increase in age at 52 months, there was a decreasing 
trend in MAI and the value was 3.27 m in agrisilviculture, 3.14 m in case of compact 
block and 3.62 m for boundary plantation, respectively. Chandra et al, 1998[34] 
reported average MAI for plant height as 2.21m in high rainfaU area (Rudrapur, 
Nainital), MAI for plant hei^t was less than 3.00 mas reported in many studies from 
different parts of the country viz, Tandon, et al, 1993[198]; Sundararaju, et al, 
1995[194]; Aulokh and Sandhu, 1990[14] and Venkatesh and Shanna, 1977[216]. At 
the age of 22 months, the MAI of dbh was 3.78cm in agrisilviculture, 3.99 cm in case 
of compact plantation and 3.86 cm in case of boundary plantation. At the age of 52 
months, the MAI of dbh was 3.41 cm in agrisilviculture, while it was 3.22 cm in case 
of compact block and 3.88 cm for boundary plantation respectively. Mean annual 
increment for dbh in this trial was quite good when compared to results of other trials 
conducted at various places (Tandon, et al, 1993 [198]; Aulokh and Sandhu, 
1990[14]). Chandra et al, 1998[34] reported average MAI of 2.27 cm of best 
provenance in 25 provenances trial of Eucalyptus species at Rudrapur, N ainital which 
is hi^ rainfall area. Promising hybrids FRI-4 and FRI-5 had MAI of 1.67 cm at 
Dehra Dun (Venkatesh and Shanna, 1977[216]). Ajit et al, 2000a [3] reported MAI 
of 3.22 m yr'^ for 3.5 yrs of Eucalyptus tereticomis clones obtained from TERI, New 
Delhi. Annandale (2002)[12] reported MAI of height as 2.5m yr’^ M K cloezicma aX 
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4.4 years, 3.3ni yr ^ for E. urophylla; more than > 2m yr’* for E. cloeziana and E. 

robusta; 1.7 m yr’* for E.peUita; 3.20m yr’* of height and 3.02 cm yr'^ for K 

camandulansis in different locations of north Queens land, Australia. Similarly Jara 

(1991)[85] reported that MAI of height ranged between 1.67-1.82 m yr’^ and MAI of 

dbh ranged between 1.20-1.72 cm yr ^ at the age of 6 years for different provenances 

ofE. saligna m Australia; MAI of height as 3.18m yr’* and MAI of dbh as 3.22 cm yr’ 

at 3 years m Austraha has been reported by Henson et al, 2007[73]; MAI of height 

was 1.04m for E. marginata clones in Australia as reported by Stukely et al., 
2007[192]. 

4.4.2 Clone evaluation 

The survival of clone C-3, C-7, C-10 was 100% in all the three systems whereas that 
of C 6 was 72.73 ^ at the age of 52 months. Hi^y significant differences in growth 
attributes were observed among the four clones, at all stages of observations, under 
compact block plantations i.e tree control without crops. However, the clonal 
differences were not significant (p=0.11) in some cases of intercropped systems 
specifically the growth attributes were statistically at par at 22 and 52 months age in 
case of boundary plantation and at wider spacings (10x2, 8x4, 5x5m) in 
agrisilviculture at 52 months of age (Table-2). Maximum growth was recorded for C- 
10 and C-3 followed by C-7 and least in C-6 (Fig. 1). In majority of cases the clonal 
differences were highly significant for growth attributes. Mact the differences were 
not significant only for intercropped systems. This may be attributed to the fact that 
in the intercropped systems the tree components indirectly receives the benefits of 
irrigation, fertilization and intercultural operations meant for crop con^onent and this 
exerts minimum competition for light, nutrients and moisture, and thus clonal 
variation is suppressed to a certain limit, whereas the clonal variations exhibits 
dominance in pure tree treatment (compact block) where there is stiff competition for 
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growth stimulators viz light, moisture and nutrients. The results obtained from this 
study indicate that genetic differences exist among clones of E. tereticomis. Kumar 
and Bangarwa, 2006[95] also reported clonal differences for all growth characters for 
E. tereticomis m Haryana. Similarly Rao et al, 2005[155] reported both inter and 
intra clonal variations among four years clones of E. tereticomis in Andhra Pradesh. 

4.4.3 Role of spacing 

Growth performance was evaluated at eight spacing levels spread over three systems. 
While comparing the average values of growth attributes (height and dbh) under 
spacing levels, 8x4 performed better compared to other spacings in agrisilviculture 
system and interestingly it was clone 3, which performed best at both ages. Like wise, 
3x3m spacing lead to better growth attributes, at all ages, as compared to 2.5x2.5 in 
compact block plantations. At two years age, the mean height was hi^er in CB 
(3x3m) because initially trees compete for light and therefore the height grows more 
than crown and diameter but at the age of 52 months optimum values of height was 
observed in the intercropped systems viz 8x4 in agrisilviculture, 2.5 spaced in 
boundary plantations. This may be attributed to the fact that wider spacing, both 
between row to row and plant to plant, result in better growth due to reduced 
competition for the want of li^t, moisture and nutrients. The analysis of growth data 
over 52 months suggests the superiority 8x4 and 10x5m in agrisilviculture over other 
spacing in terms of both dbh and height (Fig. 2). Over all, significantly hi^er tree 
growth was recorded under wider spacing in all the three systems viz, 8x4 and 10x5m 
in AS, 3x3 in CB and 2.5 in boundary plantation. Several authors reported better 
plant growth under wider spacing (Swamy et al, 2007[196] for Gmelina arborea; 
Khan and chaudhry, 2007[91] for Populus deltoides; Nasir et al, 2006[130]for Citrus 
paradisi). Nevertheless, the results are at early age of plantations in agroforestry and 
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p block suggests need of wider spacing and cultivation inputs for optimum 
growth. 

4.4.4 System Comparison 

The growth performance was statistically at par in all the three systems during the 
ini 1 period upto three years of age, however significant and distinctly higher 
growth was recorded nnder the two agroforestry systems (AS, BP) as compared to 
monocroppmg system (CB) from the third year onwards (Fig. 3). Gill and Ajit 
(2004)[60] also reported better growth of trees in agrisilviculture compared to 
silviculture. Dhyani and Tripathi, (1998)[46] also reported similar results under 
agrisilvicultural practices in northeast India. Patil et al, 2000[141] reported best 
height growth of Eucalyptus in agroforestry. This may be attributed to the feet that 
during cultivation of crops, various inputs like fertilizers, irrigation, given to the crop 
,seed bed preparation etc. indirectly boost the growth of tree component also. 
Chauhan et al., 1997 [3 8] have reported that Eucalyptus growth was improved by file 
intercrops, which closely relates to the finding of the present investigation. 

4.4.5 Development and vaUdation of growth models 

The descriptive statistics for various growth attributes pertaining to the complete 
observed dataset for the age group of 1-4.5 years, have been compiled in Table-3. 
The dbh (cm) values of trees ranged from 0.20 to 22.61 cm and height from 1.08 to 
20.04m. Since no other independent data set was available for validation of the 
models, the original data set itself was randomly splitted into two mutually exclusive 
and pseudo independent data sets containing 80% and 20 % observations, used for 
model estimation and validation respectively. The smaller p-values (<0.05) of the 
Shaphiro-Wilk’s statistic confirms the normality of these dependent variates. To get 
an idea of the shape of the function to be fitted on the data, a scatter plot of height vs 
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dbh was initially drawn. It was clear from this scatter plot that the candi^te functions 
y dopted for modeling hei^t-dbh curves viz richards, Schumacher and 
allometric/power will fit well the observed pooled data set (AS+CB+BP). The plots 

se fitted functions along with the observed data points have been corralled in 1 1 
Fig It is clear from the above computations and the graphs of the fitted functions v 
that the R values are comparable for aU the three functions and was maximum for 
allometric one. However, as has been pointed out by other researchers also in the 
recent articles (Ajit et at, 2006[6]; Prajneshu and Chandran, 2005[146]), that 1 
value alone should not be used to judge the best fitting function and it is equally 
important to consider the validation and more importantly the behavior of the fitted 

function within and outside the observed range of the independent variate. ;| 

Accordingly, to judge the prediction capabilities of these functions in the extrapolated 
region, the function curves for the three-fitted function were drawn (Fig. 5) for a ^ 

wide range of dbh values. These comparative extrapolated prediction graphs of the | 

three fitted functions suggest that the allometric function results m reasonably | 
acceptable estimations, even quite outside the observed range also, whereas the other 
two namely nchards and Schumacher leads to merely constant estimation of size for 
the extrapolated range. Therefore, the allometric functions, which meets both, the ? 
criteria s of hi^ R -value and reasonably acceptable extrapolated predictions, was | 

preferred over the other two. Accordingly, allometric models were developed 1 

separately under the three agroforestry systems on estimation dataset (Fig. 6). The j 
descriptive statistics for various growth attributes pertaining to both estimation and 
validation dataset for three agroforestry systems for the age group of 1-4.5 years, 
have been compiled in Table-4. Parameter estimates of the allometric model 
(Heigiit=a.(dbh) ) on estimation dataset (80%) under different agroforestry systems 
has been compiled in Table-5. To judge the prediction capabilities of these functions, 
the system wise models developed in estimation datasets were used to compute the 
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predicted height in the vaUdation dataset for each system. All these models were 
evaluated for their statistical vaUdation through residual diagnostics (Table-6). The 
error of prediction, termed as residual, is computed as the difference in the observed 
and predicted values. Theoretically, the residual should be independently and 
normally distributed with mean zero and constant variance. These assun^tions were 
evaluated through pertinent graphs (Fig. 7,8 and 9). The plots of residuals against the 
explanatory variate (dbh) ensures that the residuals are not being continuously over or 
under estimated. The plot of residuals against their expected values clearly portrayed 
the normality of residuals. The plot of residual against es timat e indicated that the 
residual have constant variance. The auto correlation plot of the residuals ensures that 
the residuals are not correlated. Thus the proposed aUometric equations fulfilled the 
regression requisites. Validation of the model was also ensured by fitting a linear 
equation (Predicted Height=a+b* observed height) between observed and predicted 
hei^t (Fig. 10). Theoretically, for a perfect fit model the value of ‘a; and ‘b’ must be 
zero and one respectively, along with a perfect unity R^-value. The value of ‘a’ 
approaching to zero and ‘b’ approaching to one ensures the good ness of the fitted 
model (Table-7). Moreover, the predicted hei^t values constantly increases with 
increasing dbh values. The proposed models can be used for predicting the hei^t of 
standing trees at any stage where height estimation is not possible using instruments, 
particularly in dense plantations. 

4.5 Conclusion 

The analysis of growth data over years suggests the superiority of wider spacing in 
terms growth increments for both dbh and height. Overall performance of clone- 10 
and clone-3 for all three systems was better then other clones. The growth 
performance was statistically at par in all the three systems during the mi tial period 
upto three years of age, however significant and distinctly hi^er growth was 
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der the two agroforestry systems (AS, BP) as compared to monocropping 
system (CB) from the third year onwards. This may be attributed to the fact that 
’ g ivation of crops, various inputs given to the crop like fertilizers and 
irrig n, indirectly boost the growth of tree conqponent. In model estimation, all 
us frmctions (allometric, richards and Schumacher function) based on the criteria of 
( s pred.) values, gave comparable results. However, the aUometric function 
led o the best predictions when these functions were validated on the criteria of 
extrapolated predictions within and outside the observed range. Allometric models 
were therefore separately developed for three agroforestry systems (AS, CB, BP). 
However, since model selection using nonlinear regression is an inherently subjective 
process, and somewhat data dependent, basing the final evaluation of the model on 
the validation procedure throu^ pertinent graphs and extrapolated prediction seems 
to be a viable alternative when the prime objective of modeling exercise is prediction. 
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Table-2: Growth attributes of the clonal Eucalvutus 
systems and spacings. 


plantations 


under dilSerent agmforestry 


Age 

Syste 

m 

Spacing 

Variables 

Clone 1 

Mean 

p-Yalne 


C6 

C3 

CIO 

C7 

22 

months 

AS 

5x4 

Dbh 

7.01 

6.51 

6.82 

8.22 

7.14 

0.00000 


Height 

7.76 

7.59 

i55 

8.14 

7.76 

0.019670 

10x2 

Dbh 

8.59 

7.44 

8.28 

8.43 

8.18 

0.000119 

Height 

8.71 

7.87 

8.31 

8.48 

8.34 

0.039424 


8x4 

Dbh 

7.90 

9.55 

7.58 

7.80 

8.21 

0.000455 


Height 

8.14 

8.86 

7.28 

7.63 

7.98 

0.000705 

5x5 

Dbh 

7.37 

8.75 

7.11 

7.90 

7.78 

0.004457 


Height 

8.08 

8.55 

7.27 

7.83 

7.93 

0.000067 

10x5 

Dbh 

6.89 

5.01 

6.36 

8.28 

6.66 

0.000165 

Height 

7.80 

6.96 

6.72 

7.91 

7.40 

0.000018 

CB 

3x3 

Dbh 

8.64 

8.23 

7.70 

8.53 

8.28 

0.000000 

Height 

8.41 

8.31 

8.05 

8.42 

8.30 

0.000333 

2.5x2.5 

Dbh 

6.46 

7.40 

6.65 

7.16 

6.89 

0.000024 

Height 

7.95 

8.10 

7.34 

8.02 

7.84 

0.000001 

BP 

2.5 

Dbh 

7.66 

7.80 

7.69 

7.73 

7.72 

0.918885 

Height 

7.81 

8.16 

7.86 

8.04 

7.97 

0.110120 

52 

months 

AS 

5x4 

Dbh 

13.22 

14.86 

14.95 

12.45 

13.90 

0.000260 

Hei^t 

13.16 

14.39 

14.45 

12.58 

13.66 

0.000418 

10x2 

Dbh 

ii70 

16.51 

17.38 

16.36 

16.58 

0.027328 

Height 

15.06 

15.67 

16.29 

15.61 

15.72 

0.038648 

8x4 

Dbh 

16.43 

20.70 

17.19 

17.70 

18.18 

0.013902" 

Height 

15.52 

18.28 

16.15 

16.53 

16.73 

0.038387 

5x5 

Dbh 

17.51 

17.75 

17.00 

15.44 

16.77 

0.151145 

Hei^t 

16.34 

16.63 

16.07 

14.87 

15.87 

0.136345 

i0x5 

Dbh 

16.45 

12.02 

14.73 

20.30 

15.94 

0.000000 

Height 

15.63 

12.24 

14.33 

18.28 

15.18 

0.000000 

CiB 

3x3 

Dbh 

13.94 

15.19 

15.35 

15.11 

14.88 

0.000096 

Height 

13.73 

14.68 

14.81 

14.61 

14.45 

0.000096 

2.5x2.5 

Dbh 

12.19 

13.13 

15.16 

12.49 

13.27 

0.000004 

Height 

12.38 

13.10 

14.66 

12.60 

13.21 

0.000004 

BP 

2.5 

Dbh 

17.62 

16.75 

18.61 

17.46 

17.49 

0.030314 

Height 

16.40 

15.81 

17.04 

16.37 

1 16.32 

0.056657 
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Chapter-4.- Modcffiag ami attrihrtes (he^W: <liii).... 


Table-3. Summary characteristics of growth attributes for complete dataset. 


Measured 









variable 

Cases 

Average 

(Min, Max.) 

Standard 

Deviation 

Skewness 

(SEof 

Skewness) 

Kurto^s 

(SEof 

Knrtods) 

SW 

Statistic 

SWP- 

Value 

Heignt (m) 

4667 

8.85 

(1.08-20.04) 

4.72 

-0.0447 

(0.0358) 

-1.376 

(0.0716) ' 

0.9297 

0.0000 

Ddii (cm) 

4667 

8.26 

(0.20-22.61) 

5.31 

0.0680 

(0.0358) 

-1.225 

(0.0716) 

0.9441 

0.00000 


PIiJ>. Thesis OD “DocunieBtafioii, Analysis and Modelling of Eucalyptus... at National Level” by NIGHAT JABEEN (2008) 


92 
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Table-4: Descriptive statistics 
three agroforestry systems 


of model estimation (80%) and validation (20%) datasets under 


Number of 
observations 

Measured 

variable 

Average 
(Min, Max.) 

Standard 

Deviation 

Skewness 

(SEof 

Skewness) 

Kurtosis 
(SE of 
Kurtosis) 

*SW 

Statistic 

SWF- 

Valiie 

Agrisilviculture (ctods +t 

rees within field 

i) 





Estimation dataset 

1908 

xieignt 

8.62 

(1.08-19.38) 



-1.344 

(0.1120) 

0.93427 

0.000 1 


dbh 

1.02 

(0.29-21.66) 

5.33 

0.168 

(0.056) 

-1.204 

(0.1120) 

0.94110 

■i 


Validation da 



Ltaset 




478 


8.95 

(1.08-18.48) 


lau^H 


0.93581 

0.000 


abii 

8.35 

(0.30-20.70) 

5.49 

0.130 

(0.111) 




— : — Compact Block Plantations 

(pure tree cent 

rSS^B 




^ Estimation dataj?et 


Height 

8.82 

(1.40-17.48) 

4.52 


-1.411 

(0.135) 

0.910 

0.000 


dbh 

8.14 

(0.20-18.47) 

4.94 


Qgnn|| 




Validation’ dj 

ltaset 

IHBIIIi 


■■■■ 

355 


9.36 

(1.40-16.35) 

4.53 



0.898 

0.000 


fismiin 

8.76 

(0.23-17.83) 

4.97 

ElMi 

-1.27 

(0.258) 

0.928 


Boundary Plantation (crops+ trees 

on field bound 

aries only) 


1 


Estimation dataset 

494 

Height 

9.60 

(1.40-20.04 

5.26 


-1.475 

(0.219) 


0.000 


dbh 

9.28 

(0.28-22.61) 



-1.368 

(0.219) 

0.929 



Validation dataset 






7.94 

(1.43-18.00) 

4.90 

0.314 

(0.219) 

-1.270 

(0.436) 


0.000 


dbh 

7.38 

(0.35-20.38) 

5.50 


wm 


0.00001 


*SW- Shaphiro wilk’s statistic for normality testing 
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Cliapter-4: ModeQing and comparison of growtfe attrltjntes <heiglit and 


Table-5: Parameter estimates of the aflometric 
(80%) under different agroforestiy systems. 


model (Height=a.(dbh)*^ on estimation dataset 


System 

Parameter 

Estimate 

Asymptotic 

Standard 

Error 

Wald Confidence Interval 

95% 


Lower 

Upper 

AS 

a 

1.84 

0.013 

1.814 

1.868 

0.98 

b 

0.76 

0.0029 

0.756 

0.768 

CB 

a 

1.80 

0.024 

1.761 

1.855 

0.97 

b 

0.77 

0.0054 

0.762 

0.783 

BP 

a 

1.76 

0.0271 

1.709 

1.816 

0.98 

b 

0.77 

0.0058 

0.768 

0.791 


T^le-6: Descriptive statistics of residuals computed on validation dataset (20%) 
different agroforestry systems. 


System 

No. of cases 

Average 

Standard 

Deviation 

Skewness 

(SEof 

Skewness) 

Kurtosis 

(SEof 

Kurtosis) 

sw 

Statistic 

SWP- 

Value 

AS 

478 

0.16 



1.437 

(0.222) 



CB 

355 




3.142 

(0.258) 

0.854 

0.000 

BP 

121 

0.21 

0.614 

-0.643 

(0.219) 

1.116 

(0.436) 

0.968 

0.005 


Table-7: Validation of allometric model by fitting a linear equation between predicted and 
observed height values (predicted height=a+b*observed height) on validation dataset (20%) 
under different agroforestry systems. 


System 

Parameter 

Estimate 

Asymptotic 

Standard 

Error 

Wald Confidence Interval 

95% 


Lower 

Upper 

AS 

a 

-0.14 

0.051 

-0.24 

-0.04 

0.98 

b 

0.99 

0.0050 

0.98 

1.00 

CB 

a 

0.09 

0.085 

-0.070 

0.266 

0.97 

b 

0.97 

0.0082 

0.957 

0.989 

BP 

- i 

a 

-0.01 

0.105 

-0.22 

0.189 

0.98 

b 

0.97 

0.011 

6.953 

0.998 
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Chapter-4; Modelling and comparison of growth attributes (he^ht and dbh). 


Fig. li Clone wise comparison of growth attributes over an age sequence. 
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Fig. 2: Spacing wise comparison of growth attributes over an age sequence. 
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Fig. 3: System wise comparisou of growth attributes over an age sequence, 
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Fig. 5: Extrapolated predictions of growth with different functions fitted on pooled dataset. 
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Fig. metric function fitted on individual datasets under the three agroforestry systems. 
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Chapter-4.- ModeUiag asd comparisoD of growth attriMca (hdgkt aad dMi).. 


Fig. 7: Various plots of residual diagnostics 
function under agrisilviculture system. 


on the validation 


dataset for the allometric 
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Fig. 9 : Various plots of residual diagnostics 
function fitted under boundary plantations. 


on the validation 


dataset for the allometric 
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Fig. 10: Graphs of the linear functinn u ^ 

(predicted height^ a+b*observpH ^ between the predicted and observed values 

model, ‘a’ should be zero and ‘b’ shonW agroforestiy systems. (For a perfect 

should be unity along with absolute one R^-value). 
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tereticornis 



Chapter-S: ModeUing and Comjwison rfyidd attrifcnte* (bienass and valnrae)..- 


Modelling and comparison of yield attributes (biomass and volume) of 

Eucalyptus tereticornis under two contrasting systems [agrisilviculture 
vs. compact block] 

5.1 Introduction 

Eucalyptus is used in short-rotation forestry for pulp production with a rotation length 
of 10-12 years and a productivity ranging from 10 m' ha-' year-' to 40 m' ha-' year-' 
(mean annual increment of total over bark volume at harvesting) (Antonio aL, 
2007[13]). The measurement of tree and stand biomass is of considerable interest 
today for the study of ecosystem productivity, energy flow, or nutrient and carbon 
cyclmg (Turner and Lambert, 2008[214]; Munoz et al, 2008[126]; Mani, and 
Parthasarathy, 2007[115]; Brown, 2002[23]; Snowdon et al, 2002[184]; Jenkins et 
al, 2003[87]; Hart et al, 2003[71]; Specht and West, 2003[189]; 7ianis and 
Menuccmi, 2003 [233]; Zhang et al, 2004[232]; Montagu et al, 2005[125]). Further, 
renewed interest and increasing attention is being directed towards the need to use 
woody biomass as a feedstock for bioenergy and byproducts (Volk et al, 2006[221]). 
Relative to this current fossil fuel feed stocks, biomass production systems can meet 
the growing raw material and fuel needs, creates less pollution, have positive 
environmental impacts, and maintain national security (NRC, 2000[127]). In future, 
short rotation woody crops (SRWC) will be significant source of bio mas s as more 
bioenergy and bi byproducts are produced (Bemdes, 2001 [17]). On the other band 
vegetatively propagated clonal plants developed from a single mother tree having 
most desirable genetic qualities, are uniform and true to type with aU of the genetic 
qualities of the mother plant. Comparative genetic, superiority of each clone and 
adaptability to different agro climatic and site conditions needs to be confirmed 
through scientific research trials giving equal conditions of growth and development 
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of all clones (Lai, 2007[103]). Eucalyptus clones like C3, C6, Cl, CIO and C27 
developed at Bhadrachalam formed the basis of initial clonal plantations since 1992. 
Continuous research and development support is most essential for long term 
development of clonal plantations because of many reasons. Each clone has a specific 
genetic make up and capacity to tolerate adverse conditions. Because of genetic 
uniformity of each clone, there can be serious risks of epidemic diseases and insect 
pest attacks in future if the genetic base of clonal plantations is narrow. Therefore, it 
is essential to widen the genetic base and diversity of clones by developing a large 
number of fields tested geneticaUy superior clones for commercial scale plantations 
(Lai, 2003 [102]). The area under agroforestry system of clonal Eucalyptus plantations 
is expanding fast in many states as formers are expanding production of timber and 
pulpwood to meet the growing demand. Therefore, agroforestry research with 
Eucalyptus as the prmcipal crop needs to be strengthened. Clonal Eucalyptus based 
agroforestry plantations continue to expand in many states like Andhra, Orissa, 
Maharashtra, Punjab, Haryana and U.P. Many intercrops both during the Rabi and 
Kharif season, like oats, barseem, cowpeas, ginger, turmeric and millets have been 
successfully raised by Pragati Biotechnologies at their research centre at Semi 
(Punjab) even during 5* year of clonal Eucalyptus plantations. However, further 
research on better adapted and more shade tolerant crop varieties will be very helpful 
to farmers growing Eucalyptus under agroforeshy system. 

Eucal 3 ^tus modeling has been used as an indispensable tool for prediction of volume 
and biomass in forestry and agroforestry studies. A number of studies has been 
conducted throughout world for estimation of volume growth and biomass 
conqionents of Eucalyptus species using prediction equations viz. Sochacki et ah, 
2007[188] for AustraHa; Andre er ah, 2005[11] for Congo; Soares and Tome, 
2003[186], Yuancai et ah, 2000 and 2001[229 & 230] for Portugal; Ferreira and 
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Trevisan, 2001 [571; Laclau et al ^nncriAAi ^ ^ 

’ ^ 2000[21]; Saporito, 

2003 [172] for Italy. Reeressinn a«oi ■ • 

ysis IS commonly applied for developing 

biomass ftnctioas. Tie predictor variables often nsed, are diametor a. breast heigh. 

(Malimbwi er a/., 1994[U4]; Ter-Mihaelian and Korzuhhin, 1997[200]; Brown, 
2002 [23], Landsberg and Waring 199711 ^7^^ /i,'o 

&5 'LI”7]), diameter at stem base or ankle height 

(Vertannen et al, 1993[220]; Grundv 199^r^i^^^ 

j, vjrunay, 1995[65]), crown diameter (Olsson, 

1985[135]; Hellden, 1987[72]; Tietana, I993[211]), and tree height (Stromgaard, 
1985E191]; loomis er al. 1966[110]; Peatoon er al. 1984[142]; Bartolinh .1966[16]; 
Reed and Tome, 1998[160]; Monserud and Marshall, 1999[124] ). hr most cases ftre 
non-lmear allometric equation Y=a*7^ is used for modeling tire relationship between 
tire biomass (Y) and tire predictor variable. Here ‘a’ and ‘b’ are recession 
coefficients estimated using metirod of least squares on observed data. In India, some 
woric has been recently canied out on Eucalyptus biomass production (Luna and 
Chanroh, 2008[112], Prajapatt et al. 2006[145]; Faiz et al. 2003[54]; Nagaveni and 
Vijaylahshmi, 2002[128]; Hunter, 2001[75]; Gill et ul., 2001[62] and modeling for 
Eucalyptus spp. (Tewari and Kunrar. 2003(a) [201]; Tewari et al. 2002[207]; Thaker 
and Singh, 2001[208]; Ajit et al, 2000 (a) [3]). However, all these studies have been 
conducted for pure Eucalyptus plantation only. To tire best of out knowledge, no 
report could be traced on development of growth and yield models for agrisilviculture 
and their conqrarison Sum pure plantations, before tire initiation of this study. 
Accordingly, attempts have been made hr (his study to develop the prediction 
equation for estimating above and below ground biomass, total biomass, volume and 

index of Eucalyptus tereticomis for compact block and agrisilviculture plantations 
under semi-arid conditions. 


This chapter concentrates on the twin objectives of quantifying the volume and 
biomass com^omivAs of Eucalyptus tereticomis over an age series, and to develop 
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and compare non-linear models for prediction of volnme , a w 

the basis of dbh. and biomass components on 


5.2 Materials and methods 

5.2.1 Study site 

The eitperiment is being carried out at the Central n 

n 1, r- +_ Research Farm of National 

Research Centre for Agroforestty, Jhansi which is locateH * , ■ 

, 1 1 j • . located at an elevation of 300 m 

above sea level and is situated between 24°ir N , 

, . . , . ^ w latitude and 78n7’ E longitude 

havmg tropical semi arid climate with mean annual mw n ^ ^ 

. , - , ■ 1? II ■ einfell of 900 mm. More than 

75% of the ramfeU is received during monsoon season n . , n 

, n o - n s of ft™® to first 

week of September). Mean maximim temperatiire ranges from 47.4“ C (June) to 

23.5“ C (December) and mean nnnimum temperature from 27.2“ C (June) to 4.1“ C 

(December). The soil of the experimental area is sandv i i 

loam with low orgamc 

carbon, mtrogen, phosphorous and medium to high in potassium 

5.2.2 Experimental details 

Four clones of Eucalyptus tereticornis (C3, Ce, C7 and r ^ . 

V 3, M and Cio) were obtamed from ITC 

Bhadrachalam, Andhra Pradesh. The field experiment 1 -ri . - a 

F uucnt was laid out durmg August, 

2003 m different systems and spacing, namely agrisilviculture (5x4 10x2 10x5 8x4 

5x5 m). compact block (3x3, 2.5x2.5 m) and boundary phmtetimis (^dm). Duel 
systems (agrisilviculture, compact block and boundary phmtation) were employed 
for growth smdies, however frees were harvested only from compact block 
plantation, because if flie frees wem harvested from agrisilviculture for growth/yield 
studies, the design of agrisilviculture study would have changed and accordingly a 
sepamte block of agrisilviculture were laid out wifli (he of destructive 

harvesting only. Thus eight frees each from con„a« block as well as from 
agrisilviculture were harvested after every six months, starting from the age of six 
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months upto four years age, for growth and biomass studies. Data on above ground 

and below ground growth attnbutes were recorded on each harvested plant. The 

height reported is in meters DRH in a. , • 

2 ’ ti m centimeter, biomass oven dried in kg/tree, 

volume m m/tree and volume index in mWyear and biomass productivity in 
tons/ha/year. Smce the volume and biomass are not directly con^arable under these 
two systems, because one of the variable, namely tree spacing is also involved in 


measured data. In fact, both systems have different spacing viz, 5x4m for 
agrisilviculture and 2.5x2.5m for compact block plantation. To overcome this 
problem, a more generalized measure caUed index for volume and biomass 
productivity has been used. Index is a parameter lhat does not depends on spacing 
and all other variables and accordingly these two different systems can very well be 
compared with this single measure. By the end of four years, a total of 128 cases 
were processed for total biomass and volume representing 64 harvested trees from 
both compact block and agrisilviculture respectively. Both datasets were further 
randomly split into two parts. 80% to estimate the model parameters (estimation 
dataset) and other 20% to validate the performance (validation dataset). The 
vahdation dataset acts as a pseudo independent dataset. The examinations of basic 
descriptive statistics of the 3deld (above and below biomass conqionents), volume 
and growth attributes (height and diameter at breast height) revealed that DBH and 
total biomass were approximately normally distributed. The total height was not 
considered for developing the model since usually it is difficult to measure the height 
of the standing tree with definite accuracy while one can accurately measure the 
DBH of the standing tree and thus can easily compute the biomass of the standing 
tree using the developed model. Moreover, DBH alone explained significant amount. 

of variation in volume and biomass conqionents. 
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5.2.3 Observational methods 

5.2.3.1 Above ground biomass and volume calculation 

Trees were felled 2-3 cm above the ground. The diameter at breast height (DBH) was 
recorded by measuring vernier calipers. The total height of felled tree was measured 
by tape in meter. The leaves were separated and the fresh weigjits were recorded in 
kg/tree for leaves, branches and main bole for every one-meter piece. Fresh samples 
were collected for leaf, branch and bole, oven dried at 70*^0 till constant wei^t and 
finally dry matter percentages were confuted. Fresh wei^ts were multiplied by dry 
matter percentages to obtain the oven dry biomass of leaf, branch and bole separately. 
The three components were added to yield the total above ground dry biomass. For 
the computation of timber volmne, stem from bottom to tip of the tree was 
considered. Diameter was measured on every meter hei ght along the tnaiu bole. 
Newton's formula V=(db+4dm+d,)/6)*h where db denotes diameter at base, dm denotes 
diameter at the middle, dt denotes diameter at the top and h is the length of the log, 
has been used to con:q)ute the volume of meter log and finally added to get the total 
volume (m^) of mdividual tree. 

5.2.3.2 Below ground biomass 

Althou^ some roots may extent to great depths, the overwhelming proportion of 
total root biomass was generally found within 75 cm of soil surface. For below 
ground observations, root system was harvested carefully and biomass of main root, 
primary root, secondary root and fine root were measured separately. Fresh samples 
for these components were collected and oven dried at 70°c till constant wei^t and 
finally dry matter percentage were computed. Fresh wei^ts were multiplied by dry 
matter percentages to obtain the oven dry biomass of aU con^onents separately. The 
component wise biomass was added to get the total below ground dry biomass. 
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5.2.33 Model fitting and Model vaUdation 

Non-linear regression was merl fr.r • 

j . ^ vanous functions on the estimation dataset. 

For valrdauc. ™ models can be .eated in vado. ways. Firsdy. tbe resanrpling 

approac may e suggested (Bi and Hamilton, 199[20], Bi. 1999[18], 2000[19), 

y, Iterative vaUdation procedure may be considered (Williams, 1997[226]). 

Tlnrdly, rf one has two independent datasets fiom die same area or population; one 

set can be used for fitting and the other to vaHdation or, alternatively, if one has only 

one dataset, it can be divided into two psedo independent sets through random 

procedura to flie purpose (Tiwari and Krmiar. 2003 (a)[201]). Hera the last approach 
has been adopted. 


Geisser’s (1975)[59] ‘data sphtting- pmcedure h^ been applied to vaUdate the 
model. The complete dataset under each system consisted of 64 points out of which 
4e random sample of 20% data pointe (13 trees) was selected and rescued as test 
dataset for vaUdation procedure while the rest 80% of data points (51 trees) were 
enyloyed to model selection and parameter estimation procedure. NormaUty of the 
response variable was verified using Shaiphiro-Wilks test (Rao et al., 1985[154]). 
Lmear regression was initially carried out for estimating biomass/volutiie through 
Dbh as explanatory variate, however these models suffers fiom the problem of 
negative estimation (Ajit er al., 2000 (b)[8]) of size and therefore non-linear models 
were attempted and the selection of the model was decided on validation as well as 

the behavior of model within and outside observed range. * 

5.2.4 Statistical analysis 

The statistical analysis of data was performed using SYSTAT software version 

12.01.04 (Wilkinson and Coward, 2007[225]). Analysis of variance (ANOVA) was 
used to evaluate the differences between clones and systems. The descriptive 
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Statistics (range, mean, skewness, kurtosis and Saphiro wilk statistics) of the growth 
attributes were nutiaUy computed. Non-linear regression was enqjloyed to develop 
relationship between the growth attributes and diameter at breast height as the single 
predictor variate. Validation of the model was performed initially through residual 
diagnostics, regression assumptions were verified by plotting graphs between 
observed, predicted and raw residual values. To further strengthen the validation 

procedure, simple linear regression was fitted between the observed and predicted 
values. 

5.3 Results and discussions 

5.3.1 Biomass production and volume growth 

Agrisilviculture system 

The mean dbh value in the complete harvested data set was 10.13 cm vaiydng from 
0.95 to 19.43 cm whereas the mean height value was 10.67m varying firom 2.34 to 
18.35 m (Table-1). The tree wise AGB (above ground biomass) ranged firom 0.266 to 
152.82 kg/tree and the BGB (Below ground Biomass) firom 0.11 to 46.52 kg/tree. The 
average BGB to AGB ratio was 0.17 ranging firom a mininnim of 0.01 to a TnaTcimum 
of 0.42. Looking to the component wise (bole, branch, leaf) breakup of AGB, the 
contribution of bole biomass increased firom 52% to 85% firom 1®^ to4* ycsr? whereas 
the contribution of leaf decreased firom 34% to 6%. The average contribution of leaf 
was nearly 10% in the AGB and it almost remained constant over the years (Table 3). 
The trend of biomass allocation over years has been depicted in Fig. 1. Tree wise 
volume (m^/tree) ranged firom 0.00033 to 0.221 with a mean value of 0.064. Almost 
aU the growth (height, dbh) and yield attributes (AGB, BGB, total tree biomass,, 
volume etc.) were normally distributed (normality was tested using shaphiro wilks 
criterion). 
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Compact block plantations 

The mean dbh value in the complete harvested data set was 7.52 cm varying from 
0.95 to 13.69 cm whereas the mean height value was 9.52m varying from 2.41 to 
18.30 m (Table-2). The tree wise AGB ranged from 0.30 to 70.05 kg/tree and the 
BGB from 0.10 to 36.76 kg/tree. The average BGB to AGB ratio was 0.17 ranging 
from a minimum of 0.08 to a maximum of 0.33. Looking to the component wise 
(bole, branch, leaf) breakup of AGB, the contribution of bole biomass increased from 
56% to 86% from 1® to4** year, whereas the contribution of leaf decreased from 28% 
to 4 %. The average contribution of leaf was about 10% in the AGB and it almost 
remained constant over the years (Table 3). The trend of biomass allocation over 
years has been depicted in Fig 1. Tree wise volume (m^/tree) ranged from 0.00034 to 
0.1161 with a mean value of 0.034. Almost all the growth (height, dbh) and yield 
attributes (AGB, BGB, total tree biomass, volume etc) were normally distributed 
(normality was tested using shaphiro whks criterion and smaller p-values (<0.05) 
indicates normality). 

5.3.2 Comparison of yield attributes (biomass and volume) under the two 
contrasting agroforestry systems 

Tree wise average abovegroxmd biomass was more (almost double) in agrisilviculture 
(45.85 kg/tree) as compared to compact block plantation (23.54 kg/tree). A similar 
trend was observed for volume also(0.034 mVtree in CB and 0.064 mVtree in AS. 
However, if we wish to compare the yield under the two systems at the wider level 
^lot/field), this tree wise con^jarison is not sufficient since the tree spacing under the 
two harvested systems was different (2.5x2.5 in CB and 5x4 m in AS). Thus to 
compare the systems yield at plot level (per hectare basis), we have to bring yield 
attributes at the same level (i.e. by eliminating the role of tree spacing). Biomass 
Index (tons/ha/yr) and Volume Index (m^/ha/yr) are two such measures for comparing 
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the systems. The Biomass Index varied from 1.02 to 42.72 tons/ha/yr (with a mean 
value of 18.09) under compact block plantation and from 0.45 to 24.06 tons/ha/yr 
(with a mean value of 10.78) under agrisilviculture system. A similar trend was 
observed for Volume Index ranging from 1.08 to 53.08 m^/ha/yr (with a mean value 
of 19.44) under compact block and from 0.33 to 27.82 m^/ha/yr (with a mean value of 
11.30) under agrisilviculture system 

At the age of 4 years, CB showed higher productivity (28.37 m^/ha/yr and 29.20 
tons/ha/yr) as compared to AS (20.20 m^/ha/yr and 17.87 tons/ha/yr). Few 
comparable studies have been published on the productivity of E. tereticornis like 
Sudha et al, 2007[193] reported MAI of 10.17 t/ha/yr of aboveground biomass for a 
four year rotation cycle of eucalyptus in Khamam district of Andhra Pradesh. 
Considering the allocation of biomass in the below ground component, although 
some roots were extended to great depths, the overwhelming proportion of the total 
root biomass was found within 30 cm of the soil surface. Measuring the amounts of 
biomass in roots and their turnover is an extremely time consuming and costly 
exercise. Therefore, regression equations are often used to extrapolate aboveground 
biomass to whole-tree biomass (Kurz et al, 1996[98]; Cairns • et al, 1997[27]). 
However in this study, the trees were completely excavated and therefore exact 
belowground biomass values were available for development of root biomass models 
along with aerial biomass models. 

5.3.3 Development and validation of yield models 

Various non linear functions(allometric and variants of allometric) were attempted to . 
fit estimation data set and the values of were almost comparable and in some cases 
R^ values were on higher side for allometric fimctions (Y=a*X^, but if we compare 
the prediction capabilities of these fimctions outside the observed dbh range, variants 
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of allometric were well behaved and resulted in reasonable estimation, whereas 
allometric functions lead to steep over estimation outside the observed dbh range. 
Accordingly, the parameter estimates of the fitted function (exp(a+b/dbh^l.5) under 
agrisilviculture and compact block has been compiled in Table-4 and Table-8 
respectively. The summary characteristics of the predicted values from the fitted 
model and observed values in the vahdation data set under agrisilviculture and 
compact block have been compiled in Table-5 and Table-9 respectively, whereas the 
statistics of the residuals under agrisilviculture and compact block are listed in Table- 
6 and Table-10 respectively. The lower mean value of residuals (1.77 kg/tree for total 
biomass in agrisUviculture and 0.68 kg/tree in CB respectively) supported by the 
normality of residuals confirms the appropriateness, of the models in statistical terms. 
All these models were evaluated for their statistical validation throu^ residual 
diagnostics for agrisilviculture (Fig. 3), as well as for conpact block plantation (Fig. 
5). Validation of the model was also ensured by fitting a linear equation (Predicted 
biomass=a+b* observed biomass) between observed and predicted component 
(Table-7 and Table-11 for agrisUviculture and compact block respectively). 
Theoretically, for a perfect fit model the value of ‘a’ and ‘b’ must be zero and one 
respectively, along with a perfect unity R^-value. The value of ‘a’ approaching to zero 
and ‘b’ approaching to one ensures the good ness of the fitted model (Fig. 2 and Fig. 
4 for agrisUviculture and compact block respectively). 

The predictive abUity of these models was assessed through validation dataset by 
computing various statistics based on residuals. The models obtained from estimation 
dataset was applied to the validating dataset for computing the predicted values. The 
bias is computed as die difference in the predicted and observed values in the 
validation dataset. The bias mean gives the accuracy of prediction whil e the variance 
provides information regarding precision of the prediction. Fig. 2 and Fig. 4 in AS 
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and CB respectively shows the plots of residual/bias (observed- predicted values) 
against the predicted values. No clear trend was seen in the residuals, hence it can be 
said they are not being continuously over or under estimated. The R^-values for the 
fitted models were more than 0.79 (for AGB, BGB, total biomass and volume) but on 
the basis of alone, the equations should not be recommended for use because they 
may or may not behave properly if used for prediction outside the observed range. 

This emphasizes the importance and need of validating properly a model prior to its 
use. The vahdation process is necessary so that model can be used with some 
confidence (Goulding 1979[64], Reynolds and Chung, 1986[162]). 

Further more, it would be is desirable a model (preferably allometric) for estimating 
root biomass directly from tree diameter at breast hei^t. This putative relationship is 
based on the hypothesis that the growth of roots depends on stem diameter and that 
the above and belowground stump development of the tree has an allometric balance 
(Kostler et al, 1968[93]; Santantonio, 1990[171]). The results of this study 
confirmed that dbh can be used as a simple variable to estimate the root biomass. This 
facilitates the efforts to determine root biomass without the high cost of excavating 
the root system. 

5.3.4 Comparison of yield index models (biomass and volume) under the two 
contrasting systems 

The recorded trend of the observed Index values over an age sequence has been 
portrayed in Fig. 7, while the future predictive estimation trend of the fitted model of 
Index values over the wide range of dbh-values has been depicted in Fig. 6. Index 
values (for biomass conqionents and volume) were on higher side in CB as compared 
to AS, however the difference were small in initial years, whereas as the dbh 
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increases (or age advances) the Index values were significantly greater in CB as 

compared to AS. 

In this index, the role of spacing on yield and growth has been eliminated, since 
spacing is the major determinant factor in biomass accumulation. Index values 
(biomass or volume) were on higher side in CB as compared to AS, and it was so 
simply because the number of trees were more in CB (2.5x2. 5m spacing amounting 
to 1600 trees/ha) than AS (5x4m spacing or 500 trees/ha). Although, this shght 
decrease in average Index values in AS system, were con^ensated duly by the 
intercrop biomass. Infact, if we overall calculate the economics, under the two 
systems, agrisilviculture system is more viable, economic and sustainable than CB. 
Moreover, fi:om fanners point of view, agrisilviculture system is more adaptable and 
acceptable since it starts giving cash dividends (tibrou^ crop yield) from the first 
year itself whereas in case of CB the farmers have to wait for a longer duration (at 
least up to completion of first rotation cycle). Additionally, the AS is sustainable in 
the sense that the fertilizers and irrigation apphed for the crop component sustains the 
fertility of soil up to a certain extent, whereas in CB it is only the leaf litter which add 
to the soil fertility, though very little. 
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Table-1. Summary characteristics of growth attributes for complete harvest dataset under 
agrisilviculture system. 


System 

Measured 

variable 

Average 

(Minimum, 

Maximum) 

Standard 

deviation 

Skewness 

(Standar 

d error of 

skewness) 

Kurtosis 

(Standard 

error of 

kurtosis) 

SW 

Statistic 

SWP- 

Value 

"as 

Height (mt) 

10.67 

(2.34-18.35) 

4.454 

-0.2722 

(0.299) 

-0.8875 

(0.590) 

0.960 

0.036 


Dbh (cm) 

10.13 

(0.95-19.43) 

4.728 

-0.30132 

(0.299) 

-0.5973 

(0.5904) 

0.964 

0.064 


Cd(cm) 

33.44 

(2.49-78.00) 

25.12 

0.220 

(0.299) 

-1.562 

(0.590) 

0.866 

0.00001 


Dry bole wt 

(kg/tree) 

37.30 

(0.101-134.16) 

35.75 

1.120 

(0.299) 

0.4677 

(0.590) 

0.872 

0.00001 


Dry branch wt 

(kg/tree) 

4.55 

(0.030-40.23) 

5.481 

4.564 

(0.299) 

28.46 

(0.590) 

0.602 

0.00000 


Dry leaf wt 

(kg/tree) 

3.99 

(0.135-12.44) 

2.467 

6.6717 

(0.299) 

0.9347 

(0.590) 

0.960 

0.0397 


Above ground 

biomass 

(kg/tree) 

45.85 

(0.266-152.82) 

40.36 

1.0621 

(0.299) 

0.452 

(0.590) 

0.887 

0.00003 


Belowground 

Biomass 

(kg/tree) 

14.15 

(0.113-46.52) 

11.72 

0.550 

(0.299) 

0.343 

(0.590) 

0.925 

0.00084 


Total Biomass 

60.01 

(0.44-184.34) 

50.42 

0.848 

(0.299) 

0.0312 

(0.590) 

0.911 

0.00022 


Volume 

(m^/tree) 

0.064 

(0.00033- 

0.2219) 

0.0592 

1.028 

(0.299) 

0.340 

(0.590) 

0.8894 

0.00003 


Volume index 

(m^/ha/yr) 

11.30 

(0.334-27.82) 

7.13 

0.429 

(0.299) 

-0.238 

(0.590) 

0.960 

0.0366 


Biomass index 

(t/ha/yr) 

10.78 

(0.450-24.06) 

10.78 

6.123 

(0.299) 

-0.549 

(0.590) 

0.9717 

0.0148 
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Table-2. Summary characteristics nf in- \<.i. 

compact block. * attributes for ♦ u 

•«es tor complete harvested dataset under 



0.00002 


0.00105 


0.00040 


0.00715 


0.00789 
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Table-3: Component wise allocation of biomass over an age series for both systems 
(agrisilviculture and compact block). 


Syste 

ms 

Year 

% Breakup of 

total biomass 

% Breakup of total 
Aboveground biomass 

BGB- 

AGB 

iN^et Biomass I 

Productivity 

"^et Volume 

Productivity 



AG 

BG 

Bole 

Branch 

Leaf 

ratio 

(tons/ha/yr) 

(m^Aia/yr) 

AS 

1 

7.419 

(76%) 

2.291 

(24%) 

4.160 

(56%) 

1.206 

(16%) 

2.053 

(28%) 

0.15 

4.85 

3.87 


2 

37.73 

(69%) 

16.85 

(31%) 

27.32 

(72%) 

4.938 

(13%) 

5.471 

(15%) 

0,22 

13.64 

13.12 


3 

50.75 

(68%) 

24.86 

(32%) 

44.16 

(86%) 

3.29 

(7%) 

3.29 

(7%) 

0.28 

15.60 

16.30 


4 

116.57 

(82%) 

26.12 

(18%) 

99.94 

(86%) 

12.21 

(10%) 

4.41 

(4%) 

0.14 

17.87 

20.20 

CB 

1 

1.840 

(79%) 

0.495 

(21%) 

0.960 

(52%) 

0.2672 

(14%) 

0.613 

(34%) 

0.11 

3.73 

3.81 


2 

22.59 

(80%) 

5.272 

(20%) 

15.823 

(70%) 

2.2887 

(10%) 

4.481 

(20%) 

0.11 

22.29 

22.17 


3 

40.48 

(75%) 

13.41 

(26%) 

31.55 

(78%) 

4.19 

(11%) 

4.73 

(12%) 

0.20 

28.74 

27.74 


4 

51.33 

(72%) 

22.03 

(28%) 

43.76 

(85%) 

4.40 

(8%) 

3.16 

(6%) 

0.22 

29.2 

28.37 
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Table-4: Values of coefficients for different biomass components obtained for jBtfing dataset of 
agrisilviculture Y=e (a+b/x • ) (Y-component biomass (kg/tree) and X-dbh (cm) 


Biomass 

Component 

(ag) 

Parameter 

Estimate 

Asymptotic 

Standard 

Wald Confidence Interval 

95% 


Error 

Lower 

Upper 

Bole 

a 

5.99 

0.069 

5.85 

6.13 

0.95 


b 

-94.32 

3.902 

86.49 

102.15 


Branch 

a 

4.14 

0.351 

3.44 

4.85 

0.47 


b 

-109.74 

20.63 

68.34 

151.41 


Leaf 

a 

1.80 

0.118 

1.567 

2.04 

0.37 


b 

-8,20 

3.717 

0.739 

15.66 


Total 

a 

6.05 

0.064 

5.92 

6.18 

0.96 

Aboveground 

b 

-87.52 

3.551 

80.39 

94.64 


Total 

a 

4.21 

0.129 

3.95 

4.47 

0.79 

Belowground 

b 

-55.28 

6.318 

42.60 

67.96 


Total biomass 

a 

6.15 

0.056 

6.04 

6.26 

0.96 


b 

-78.97 

3.024 

72.90 

85.04 


Volume 

a 

-0.50 

0,069 

-0.641 

-0.362 

0.95 


b 

-87.12 

3.806 

79.48 

94.76 


Index volume 

a 

3.70 

0.081 

3.53 

3.86 

0.87 

(m%a/yr) 

b 

-43.32 

3.756 

35.78 

50.86 


Index biomass 

a 

3.37 

0.084 

3.20 

3.54 

0.82 

(t/ha/yr) 

b 

-31.68 

3.612 

24.43 

38.93 
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Table-5: Summary statistics of the predicted and observed values in the validating dataset 
(agrisilviculture) 


Measured variable 

(pooled) 

Average 

(Min, Max.) 

Standard 

Deviation 

Skewness 

(SEof 

Skewness) 

Kurtosis 

(SEof 

Kurtosis) 

sw 

Statistic 

SWP - 

Valne 

dbh 

8.59 

(0.95-16.88) 

5.84 

0.214 

(0.687) 

-1.612 

(1.334) 

0.916 

0.332 

Obs. aboveground 

38.44 

(0.26-121.42 

43.69 

1.00 

(0.687) 

-0.296 

(1.334) 

0.838 

0.042 

Pred. aboveground 

37.94 

(0.002-120.06 

48.66 

0.851 

(0.687) 

-1.100 

(1.334) 

0.779 

0.008 

Obs. belowground * 

“12.37 

(0.18-41.86) 

14.62 

1.07 

(0.687) 

0.021 

(1.334) 

0.821 ^ 

0.026 

Pred. belowground 

11.08 

(0.003-30.35) 

12.75 

0.562 

(0.687) 

-1.662 

(1.334) 

0.801 

0.015 

Obs. total biomass 

50.81 

(0.44-149.14) 

55.65 

0.770 

(0.687) 

-1.062 

(1.334) 

0.836 

0.0395 

Pred. total biomass 

49.04 

(0.001-150.10) 

61.39 

0.783 

(0.687) 

-1.246 

(1.334) 

0.785 

0.0096 


Table-6: Descriptive statistics of residuals in the validation dataset (agrisilviculture). 


Residuals 

Average 

(Min, Max.) 

Standard 

Deviation 

Skewness 

(SEof 

Skewness) 

Kurtosis 

(SE of 

Kurtosis) 

SW 

Statistic 

SW P- 1 

Value 

Aboveground 

0.49 

(-15.40-12.85) 

8.72 

-0.539 

(0.687) 

-0.377 

(1.334) 

0.954 

0.723 

Belowground 

1.35 

(-3.18-7.6) 

109 

0.478 

(0.687) 

1.017 

(1.334) 

0.950 

0.670 

Total biomass 

1.77 

(-21.25-9.71) 

9.18 

-1.89 

(0.687) 

4.655 

(1.334) 

0.765 

0.0050 
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Cliapter-5: Modeling and Comparisoii of yield attributes (biomass volum ) 


Tablc-7: Component wise summary characteristics of model parameters (agrisilvicultur ) 


.. — ; — r 

1 

psysteBa 

(Linear) 

"Totef^ ~ 1 

aboveground 

Parameter 

a 1 

Estimate 

-4.35 

Asymptotic 1 

Standard 1 

Error j 

3.45 t 

Wald Confid 

95% 

Lower 

-12.31 

Upper 

3.61 

"0.97 

6.73 

fb 1 

1 "i 

1 AQ 

0.061 j 

3.434 

0.95 

-6.42 1 

1.24 

9.41 

Total 

belowground 

a 

b 

1 L.^y 

0.98 

0.209 

0.49 

1 1.46 

0.98 \ 

^Xotal biomass 

a 

6.68 

2.990 

-0.20 

13.58 


r 

0.89 

0.039 

0.80 

0.99 

J 1 

1 — — — 

J — « — 



Fig 1; Year wise aUocation of biomass in different above 
umier two agroforestry systems. 


and below ground components 
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Chapter-S. Modelling and Comparison of yield attritjutes (biomass and volume).... 

Fig. 2: Graphs of the linear function u 

(predicted component^ a+b*observed between the predicted and observed values 

component) under agrisilvicnlture. 
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and valmine)..,. 


Table-8: Values of coefficients for different hiAmo», 

compact block (Y=e (a+b/x^^ (Y-comoonent hh?m obtained for fitting dataset of 

F MX component biomass (kg/tree) and X-dbh (cm) . 


Biomass 

Component 

Parameter 

Estimate 

Asymptotic 

Standard 

Error 

Wald Confidence Interval 

95% 

■r' 




Lower 

Upper 


Bole 

a 

5.09 

0.1116 

4.87 

^32 

0.92 


D 

-56.37 

3.929 

-64.26 

-48.47 


Branch 

a 

2.14 

0.252 

i.633 

2.64 

0.45 


b 

-30.85 

8.015 

-46.96 

-14.75 


Lcai 

a 

T89 

0.185 

1.521 

2.26 

0.48 


b 

-16.81 

5.236 

-27.33 

-6.29 


Total 

Aboveground 

a 

U ■" 

5.11 

0.1178 

4.876 

5.34 

0.91 

b 

-49.41 

4.054 

-57.56 

-41.26 


Total 

Belowground 

a 

t- 

4.57 

. 0.2210 

4.131 

5.01 

0.79 

b 

-68.15 

8.040 

-84.31 

-51.99 


Total biomass 

a 

■u 

5.53 1 

0.125 

5.281 

5.78 

0.90 


b 

-53.76 

4.365 

-62.53 

-44.99 


Volume 

a 

• * 

-1.13 

0.084 

-1.30 

-0.96 

0.96 


b 

-58.88 

3.000 ~ 

-64.91 

-52.85 


Index volume 

a 

4.35 

0.108 

4.139 

4.57 

0.87 

ftia/yT) 

b 

-32.25 

3.483 

-39.26 

-25.25 


Index biomass 

(t/ha/yr) 

a 

3.99 

0.110 

3.773 

4.217 

0.81 

b 

-23.24 

3.339 

-29.25 

-16.53 
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Table-9: Summary statistics of the predicted and observed values in the validation dataset 

(compact block). 


variable 

Average 

(Min, Max.) 

Standard 

Deviation 

Skewness 

(SE of 

Skewness) 

Knrtosis 

(SEof 

Knrtosis) 

sw 

Statistic 

SWF- 

Value 

"dbh 

5.77 

(1.11-12.74) 

4.68 

0.266 

(0.637) 

-1.975 

(1.2322) 

0.813 

0.0135 

■^bs. aboveground 

18.84 

(0.33-67.42) 

23.06 

0.994 

(0.637) 

-0.0918 

(1.2322) 

0,807 

0.0110 

Pred. aboveground 

17.58 

(0.003-55.89) 

21.37 

0.651 

(0.637) 

-1.348 

(1.232) 

0.787 

0.0068 

Obs. belowground 

5.42 

(0.10-16.08) 

6.44 

0.748 

(0.637) 

-1.181 

(1.232) 

0.787 

0.0068 

Pred. belowground 

5.98 

(0.001-21.57 

7.74 

0.901 

(0.637) 

-0.568 

(1.232) 

0.783 

0.0061 

Obs. total biomass 

24.26 

(0.52-79.38) 

28.79 

0.781 

(0.637) 

-0.895 

(1.232) 

0.809 

0.0120 

Pred. total biomass 

23.58 

(0.002-77.30) 

29.10 

0.708 

(0.637) 

-1.187 

(1.232) 

0.787 

0.0067 

. 


Table-10 1 Descriptive statistics of residuals in the validation dataset (compact block). 


Residuals 

Average 

(Min, Max.) 

Standard 

Deviation 

Skewness 

(SEof 

Skewness) 

Kurtosis 

(SEof 

Kurtosis) 

SW 

Statistic 

SWP- 

Value 

Aboveground 

1.25 

(-7.57-11.52) 

4.48 

0.544 

(0.637) 

2.815 

(1.232) 

0.873 

0.0731 

Belowgroimd 

-0.561 

(-9.61-4.19) 

3.30 

-1.917 

(0.637) 

5.574 

(1.232) 

0.757 

0.0032 

Total biomass 

0.681 

(-3.51-5.54) 

2.40 

-0.161 

(0.637) 

1.440 

(1.232) 

0.873 

0.0723 


PhJI. Thesis on ^^Oocnmenti^oii, luni Mo4el^ of £a€afyptiis...at Nation^ Level” by NIGHAT JABEEN (200$) 




Chapter-5: MotteOiog and Comparlsou of j^eldl attributes (biomass and volra«).... 


Table-11: Validation of non-linear model by fitting a linear equation between predicted and 
observed values (predicted component =a+b*observed component) on validation dataset 
(20%) under compact block plantation. 


System 

(linear) 

Parameter 

Estimate 

Asymptotic 

Standard 

Error 

Wald Conlldence Interval 

95% 


Lower 

Upper 

Total 

a 

0.42 

1.586 

-3.105 

3.963 

0.96 

aboveground 

b 

0.91 

0.054 

0.788 

1.032 


Total 

a 

0.07 

1.311 ! 

-2.846 

2.998 

0.82 

belowground 

b 

1.08 

6.159 

, 6.733 

1.445 


Total biomass 

a 

-0.85 

0.967 

-3.015 

1.295 

0.99 


a 

1.00 

0.026 

0.948 

1.066 
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Chaptcr-5: Modelling and Compamon of jdeM attribiifes (biomass and volume)..,. 


Fig. Graphs of the linear function fitted between the predicted and observed values 
(predicted coniponent= a+b*observed component) nnder compact block. 
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Fig. 5 : Various plots of residual diagnostics on the vaKdation dataset for the non-Unear 
function under compact block plantation. 
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Chapter-5: Modelling and Compadsoa of yield attributes (biomass and volume)-.. 


Fig. 6. Predicted frend of Eucalyptus volume and biomass yield over a wide range of dbh 
under agnsilviculture (AS) versus compact block (CB) systems. 
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Chapter-5: Modelling and Con^arison of yield attributes (biomass ami votame).... 


Fig* 7: Observed trend of Eucalyptus volume and biomass yield under agrisilviculture (AS) 
versus compact block systems (CB) over an age series. 
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Chapter-6; Devdopment af Comoidated yidd no4els (voIimc ajad biomiss)..-. 


Development of consolidated yield models (volume and biomass) of 
Eucalyptus spp. at national level by employing primary, secondary 
and simulated data. 

6.1 Introduction 

Among the fast growing species in India, eucalyptus have assumed great importance 
in large scale plantation programme to meet the rising demands of wood for local and 
industrial consumption (Singh et a/., 1992[180]). Eucalyptus species has been widely 
accepted by the farmers in India not only in pure plantations but also even with more 
enthusiasm along with crops i.e. in agroforestry because of its short rotation cycle, 
rapid growth and excellent fibre qualities. Eucalyptus species has also been 
recommended among the six tree species to be used for agroforeslry in India as 
identified by the Planning Commission through its “Task Force on Greening India” 
that pubhshed its report in June, 2001. Fast-growing plantations will be required to 
fill in die shortfall of supply from natural forest. The yields that have been obtained in 
existing plantations (often with species exotic to the site) have typically been many 
times greater than those in natural forests. Thus, plantations have, additionally, the 
possibility of satisfying the current demand for wood from a much smaller land area 
than at present. Rapid and easily implemented methods are needed for the assessment 
of standing biomass to estimate the degree of carbon sequestration by forest 
ecosystems. Equations, that provides accurate predictions of volume, without local 
bias over the entire range of diameter are one of die basic building blocks of a forest 
growth and yield simulation system (Bi and Hamilton, 1998[20]). Development of 
sound management practice is one of the major priorities of the forestry sector. The 
most often used, for determining individual tree volume or biomass are the allometric 
relationships (Causton, 1985[31]; Reed and Green, 1985[159]; Thomley and Johnson, 
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1990[210]; Reed and Tome, 1998[160]; Whittaker and Woodwell, 1968[223]). 
Normally, the volume or biomass of a tree is predicted as an equation of some easily 
measured variable, such as diameter at breast height (dbh or D) or height (H) or a 
combination of them say D H.Whenever there is need for estimating the biomass of 
individual trees, the abundance of currently available predictive equations provides 
an alternative to the destructive sampling of trees for developing local equations. 
Comprehensive collections of stand-specific biomass equations are available in the 
literature (for North America Ter-Mikaelian and Korzukhin 1997[200], Jenkins et al, 
2004[86]; for Australia Eamus et al, 2000[51], Keith et al, 2000[90], Snowdon et 
al, 2000[183] and for Europe Zianis et al, 2005 [234]). Most published biomass 
equations were developed using trees sampled from specific study sites or from sites 
that represent small regions only. As a result, use of existing volume or biomass 
equations with forest inventory data at large spatial scales is unreliable because the 
equations of previous studies may be site-specific, often disorganized and sometimes 
inconsistent (Pastor et al, 1983/1984[140]; Jenkms et al, 2003[87]; Wirth et al, 
2004[228]). Furthermore, unless an equation was developed exclusively for the 
species and study region of interest under conditions typical for the study site, it is 
impossible to know which equations to choose for a particular species and site. 
Reliable models of tree-level biomass to be used at large spatial scales are available 
mainly for Scandinavia and were compiled by Marldund (1987)[1 16]. Thus there is a 
need to develop generalized biomass and volume equations of eucalyptus, that will 
hold good at a wider spatial level in the country. 

There may be additional value in deriving generalized equations for comparing 
different regions, since it is important to know that inconsistencies between regions 
are not due solely to the application of different regression equations, which yield 
contrasting values even for the same dataset but do not significantly vary among 
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themselves (Pastor e? at, 1983/1984)[140]. Despite these inconsistencies, or perhaps 
because of them, the need is clear for a consistent method and for generalized 
equations to estimate forest biomass at large scales (Jenkins et at, 2003)[87]. 

Biomass equations or volume equations are now widely developed in forestry and 
agroforestry for both scientific and industrial purposes. Crow and Schlaegel 
(1988)[41] cautioned that there is an inherent danger in using site specific allometric 
equations as they are generally not applicable to other locations. Thus, site-specific 
allometric equations are of limited utility beyond the site for which they were 
developed as they cannot be extrapolated to other locations or throughout the 
biogeography range of a particular tree species with a high degree of confidence (Law 
et at, 2001)[108]. Pastor et at, (1983/1984)[140] developed statistical approach to 
develop generalized allometric equations that are of greater utility to forest ecologists 
and plant scientists because they are usually valid for a large portion or whole of a 
tree species’ biogeographic range. The development of generalized allometric 
equations is essential to the work of other forest ecologists, plant scientists, and 
biogeographers, especially for wide-ranging and economically important tree species, 
because they allow for comparisons of biomass and production among sites (Pastor et 
at, (1983/1984)[140] and provide die requisite data for model input. Since most 
scientists have not published the raw data from which their volume or biomass 
equations were developed. However, several authors reported that generalized 
regressions developed from field data can reasonably predict the biomass or volume 
of trees from other sites (Landsberg and Waring, 1997[107]; Schmitt and Grigal 
1981 [173]; Wirth eta/., 2004[228]). 

The purpose of this study was to develop a generalized allometric equation for 
estimating and predicting the aboveground biomass and Volume of Eucalyptus 
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species on the basis of dbh and validate the derived equation using a separate set of 
destructively sampled Eucalyptus trees. To date, no known study in India has derived 
a generalized allometric equation predicting eucalyptus volume or biomass on the 
basis of dbh from site-specific regressions throughout its biogeographic range, 
despite its importance for other work in the forestry and the plant sciences. 
Accordingly attempts has been made in this chapter to develop and validate 
generalized equations diat are applicable for a wide range of site and climatic 
conditions. The resulting equation will be of interest and use to researchers, planners 
and mangers to estimate volume accumulated by the standing eucalyptus trees at any 
stage/ age by simply measuring the dbh alone. 

6.2 Data 

Site-specific aboveground biomass (dry weight) and volume equations have been 
developed for eucalyptus species in different parts of the country. Equations relating 
volume of the stem or biomass of the tree component (stem, branches, fohage and 
total aboveground) to dbh were compiled (Table-1) from the comprehensive hterature 
published in journals pertainiug to Indian context only. Several equations based on 
independent tree samples from different sites were reported, all were included in this 
study. In the present study we analyzed the regression equations from following 
states of India: Uttar Pradesh, Uttrakhand, Karnataka, Rajasthan, Punjab and Tamil 
Nadu. In accordance with the methodology set forth by Pastor et al., 1984, five 
equally spaced points were generated from each of the compiled regression (Table-1) 
and volume/AGB was computed. These points covered the range of dbh values 
specific to each regression, and were then used to generate a new simulated dataset. 
Thus for each study, five volume/AGB poiats were computed, yielding a total of 1 12 
points in the simulated dataset from which the generalized allometric regression 
equation was derived. This method ensures that each site-specific allometric equation 
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was weighed equally in the formulation of the generalized equation. Validation of the 
model was performed using a geophysically and statistically independent datasets 
both for volume and aboveground biomass. The validity of the mathematical and 
physical assumptions used in developing a model and estimating the coefficients is 
less open to question if die model gives accurate predictions of new data In effect, 
the collection of new data provides an overall check on the entire model construction 
process. 

6.3 Statistical analysis 

Systat-12 statistical software (Wilkinson and Coward, 2007[225]) was used for 
computation of descriptive statistics (mean, standard deviation, skewness, kurtosis, 
Shaphiro Wilk s statistic for normality testing etc) and fitting of different equations 
(estimates of model parameters, asymptotic standard error of estimate, confidence 
interval, R values etc) and plotting of various graphs pertaining to residual 
diagnostics (probability plot of residuals, auto correlation plots of residuals, plot of 
residuals against their expected values, plot of residual against independent variate). 
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Table-1: Site-specific published equations in India for Eucalyptus species relating dbh to 
volume 


's. 

No 

Species 

State 

Equation 

a 

b 

c 

F 

DBH 

Reference 


'~E, globulus 

Tatmlna 

du 

CVolmne) oveAark” 0 

+ b*D 

-0.1219 

0.034433 

- 

0.97 

6-30 

Jain eta/., 1993(a)[80Tl 

"2 

'E. hybrid 

Uttar 

Pradesh 

(Volume)overt.Mk=a+ 
b*D + c*D^ 

-0.0208 

0.006597 

0.000 

9 

0.96 

9-39 

Sin^cf ai., 1995[179] 

3 

~E. hybrid 

Uttar 

Pradesh 

(Volume) oveibaik~ a 

+b*D* 

-0.0903 

0.000801 

- 

0.96 

9-39 

Sin^ctaZ., 1995[179]i 

4 

'IT 

camoldulen 

sis 

Uttranc 

hal 

V= a + b*D^ 

-0.0461 

0.000874 

- 

0.97 

6-30 

Jainef a/., 1993(b)[81] 

5 

E 

.camoldulen 

sis 

Rajasth 

an 

V(Twob)=a + b*D 

-0.1188 

0.031077 

- 

0.98 


TewariefaZ., 

2001[205] 


E. species 

Puniab 1 

Vv = a+b*DBH 

-0.1545 

0.035370 


0.97 

7-32 

Jainet aL 19911791 

s 

E. hybrid 

Puniab 

Vv = a + b*DBH 

-0.2297 

0.03417 


0.89 

14-26 

Jain efaZ., 19911791 

E. hybrid 

Punjab 

Vv = a + b*DBH 

-0.3015 

0.04032 


0.90 

10-16 

Jainef a/., 1991 [79] 

m 



# = a + b*DBH 

-0.2582 

0.03606 


0.96 

10-32 

Jain ctaZ., 19911791 

d 

■ 

WMSBM 

lEsmsi 

IE^BEEmiiSSH 

llllitfiWIM 

0.04013 


0.67 

24-32 

Jain eZaZ., 19911791 

E. hybrid 

Punjab 

LogV=a + b*Log 

D 

Raw 

dataset 


" 

- 

7.1- 

43.9 

Dogra and Sharma, 
2003[47] 


E, hybrid 

Uttranc 

hal 

Vv = a + b*DBH 


2.511 

- 

0.99 

10-40 

Chaturvedi and 

Khanna, 19821361 

13 

E.species 

Uttranc 

hal 

VV = a + b*DBH 

-0.0868 

2.8335 

- 

0.99 

1040 

Prasad, 1984[147] 

14 

E, 

camoldulen 

sis 

Rajasth 

an 

V=a*D'’ 

0.00016 

2.412 


0.99 

_ 

5-52 

Tewari and Kumar, 
2003(b)[204] 

_ 


Table -2: Site-specific published equations in India Eucalyptus species relating dbh to above 
ground biomass 


s. 

No 



Equation 

a 

b 


DBH 


■ 

E. tereticomis 

Uttar 

Pradesh 

Y = a + b*(DBH) 

-149.08 

15.99 

0.80 

140 

Ajit et ah, 2006[6] 

9 



Y = a + b*(DBH) 

-376.31 

28.05 

0.60 

140 

Ajit et al., 2006[6] 

9 



migggiiiiii 


2.17 

0.70 

140 

Ajit ef a/., 2006[6] 

■ 

E, tereticomis 

Uttar 

Pradesh 


0.42 

1.96 

0.82 

140 

Ajit ef a/., 2006[6] 

m 

E. hybrid 

Haryana 

imgmuii 

0,14 


1^ 

M 

Tandon ef fl/., 1993[198] 



Kamatalca 


-635.41 

55.78 

0.95 

1 342 

■ iiiiiiliiiliiii— 

■ 

E. globulus 

Tamil 

nadu 

Raw data 

Negi and Sharma, 1987[132] 

Ll_ 

K tereticomis 

Karnataka . 

Raw data 

Knshalapna, 1 9841991 
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6,4 Results and Discussions 
6 . 4.1 Volume simulation 

Xhe simulated data set, covering the pseudo generated harvested data pertaining to 
five states (Uttar Pradesh, Tamil Nadu, Uttranchal, Punjab, Rajasthan), was generated 
from the site specific equations (secondary data) listed in Table-1. Additionally, a 
primary data set (Ajit et al., 2000c[4]) consisting of thirteen harvested Eucalyptus 
trees was also included in the simulated data set. The dbh and tree wise volume 
ranged from 5.1 to 52.2 cm and 0.0078 to 2.203 m^/tree respectively (Table-3) in the 
simulated dataset. Perfect normality was observed for volume, whereas dbh was 
approximately normally distributed (sw. p. value=0.10). The parameter estimates 
along with other related statistics, for the generalized equation relating volume to 
dbh, has been compiled in Table-4. The fitted generalized allometric equation (Fig. 1) 

is 

Y=0.00093*(X)^-®^ (1) 

where Y= volume (m^/tree) and X=dbh (cm). 

Fig.4. depicts a case by case comparison of observed values generated from the site 
specific allometric equations, predicted values from the d.erived generalized 
allometric equation (i.e. Eq. (1)), and the raw residuals. Positive raw residuals 
indicate over prediction by the generalized allometric equation whereas negative raw 
residuals indicate under prediction by the generalized allometric equation (Fig. 4). 
The relatively small residual values depicted m Fig. 4 confirms the appropriateness 
and utility of generalized allometric equation to predict the volume of standing 
Eucalyptus trees over a very wide range of spatial (location) and temporal (age) 

range. 
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Xhe validation dataset, used for validation of the generalized volume equations, 
pertains to sequential harvesting of trees (tissue culture raised clones obtained from 
Tata Energy Research Institute, New Delhi) from a designed experiment conducted at 
NRCAF, Jhansi encompassing a dbh range of 4.05 to 9.30 cm and volume of 0.0033 
to 0.0335 m /tree (Table-5) witib perfect normality of both tree variate i.e dbh and 
volume (sw. p values < 0.04). 

The residuals (defined as the difference of predicted values from generalized 
regression and observed values in the validation data set) were distributed normally 
with a mean volume values of 0.00205 m^/tree (Table-6). The validation of the model 
was further strengthened by fitting a linear equation (predicted values= a+b* 
observed values) between the predicted and observed values (Fig. 2). For a perfect fit 

I 

model the ‘a’ value should approach to ‘zero’ and ‘b’ value should approach to ‘one’ 
with a perfect unity R^- values. In this case, the value of ‘a’ was observed to 0.00685 
and ‘b’ as 0.564 with R^- value of 0.91(Table-7). The error of prediction, termed as 
residual, is computed as the difference in the observed and predicted values. 
Theoretically, the residual should be independently and normally distributed with 
mean zero and constant variance. These assumptions were evaluated through 
pertinent graphs (Fig. 3). The plots of residuals against the explanatory variate (dbh) 
ensures that tiie residuals are not being continuously over or under estimated. The 
strong agreement between the derived generalized allometric equation and validation 
with new data set is further sustained by low mean absolute error of 0.00414 m /tree. 


A 


PhJX Tbesis OB “Socaaoaatidiea, AmOyils and Mod^i^ (MFETCalyi>tiis...at National LeveT by NIGHA.T JAMEEN (2008) 


140 


Chaptcr-6; Development of Consoiilale^yldi^ models (volmne aad biomass).... 


Table-3: Summary characteristics of growth attributes for volume-dbh simulated dataset. 


variable 

No. of 
Cases 

Average 
(Min, Max.) 

Standard 

DeviatioE 

Skewness 

(SEof 

Skewness) 

Knrtosis 

(SEof 

Knrtosis) 

*sw 

Statistic 

SWF- 

Value 

■dS 

112 

23.74 

(5.1-52.2) 

9.642 

0.1475 

(0.2284) 

-0.5263 

(0.4530) 

0.9807 

0.10693 


112 

0.548 

(0.00780- 

2.203) 

0.421 

0.93849 

(0.22843) 

1.13038 

(0.45309) 

0.93226 

0.00003 


*SW- Shaphiro wilk’s statistic for nonnality testing 


Table- 4; Allometric model fitted on the Volume-dbh simulated dataset. 


Parameters 

Estimate 

Standard 

Confidence Interval 

R* 


Value 

Error 

Minimum 

Maximum 



0.00093 

0.00027 

0.00038 

0.00147 



1.96 

0.0837 

1.8024 

2.13449 

0.87 


Table-5: Summary Characteristics of grovrth attributes for the independent validation 
dataset. 


Measured 

variable 

Average 
(Min, Max.) 

Standard 

Deviation 

Skewness 

(SEof 

Skewness) 

Knrtosis 
(SE of 
Knrtosis) 

sw 

Statistic 

SWP- 

Value 



1.54 


-0.562 

(1.090) 

0.948 




0.00955 



0.918 

0.0161 


Table-6: Descriptive statistics of residuals for the volume-dbh generalized regression 
equation. 


Average 
(Min, Max.) 

Standard 

Deviation 

Skewness 

(SEof 

Skewness) 

Knrtosis 

(SEof 

Knrtosis) 

SW 

Statistic 

SWP-Value 

0.00205 

(-0.00541-0.00833) 

0.00446 

-0.4138 

(0.564) 

-0.9779 

(1.090) 

0.935 

0.029733 


PIlD. niesis cm ‘‘Docimmatatioii, Andj^sis md ModeKi^ ofEucstfPt®—®* National Ijeve!” by NIGHAT JABEEN (200$) 


141 












Cliapter-6: Development of Com:^idbte<i yield models (vcdmne awi Momass)-... 


Table-7: Validation of allometric model on new independent dataset by fitting a linear 
equation between predicted and observed volume values (predicted volunie=a+b*observed 
volume). 


Parameter 

Estimate 

Asymptotic 

Standard 

Error 

Wald Confidence Interval 

95% 


Lower 

Upper 

a 

0.00685 

0.0010 

0.00467 

0.00903 

0.91 

b 

0.564 

0.0452 

0.4677 

0.66186 


Fig. 1* Points generated from individual regressions for Eucalyptus volume pertaining to 
various states of India and the resulting generalized volume regression line. 



Fig. 2. Graphs of the linear function fitted between the predicted and observed values 
(predicted volume= a+b*observed volume) for Eucalyptus species in India. (For a perfect 
model, ‘a’ should be zero and ‘b’ should be unity along with absolute one R^-value). 
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Fig. 3: Various plots of residual diagnostics on the validation dataset for the generalized 
volume regression allometric function. 



Autocorrelation Plot 



Lag 
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Fig- Case wise plot of the predicted and observed values with raw residuals pertaining to 
the validation dataset, for volume-dbh generalized equation. 



Independent validation dataset 


Case Number 


0 

1 
1 
I 

3 

■© 

> 


obsvol 
predvol 
i residual 
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6.4.2 Above Ground Biomass simulation 

Xhe simulated data set was generated from the site specific equations covering four 
states namely Uttar Pradesh, Haryana, Karnataka and Tamil Nadu as listed in Table- 
2. Additionally, a primary data set (Negi and Sharma, 1987[132]; Kushalappa, 
1984[99]) consisting of total fifteen harvested Eucal 3 q)tus trees was also included in 
the simulated data set. The dbh and tree wise aboveground biomass (AGB) ranged 
from 1 to 42 cm and 0.260 to 778.83 kg/tree respectively (Table-8) in the simulated 
dataset. Perfect normality was observed for biomass, and dbh (sw. p values < 0.00). 
Xhe parameter estimates along with other related statistics for the generalized 
equation relating AGB to dbh, has been compiled in Table-9. The fitted generalized 
allometric equation (Fig. 5) is 

Y= 0.74*(X)^-^^ (2) 

where Y= Aboveground Biomass (kg/tree) and X=dbh (cm). 

Fig. 8 depicts a case by case comparison of observed values generated from the site 
specific allometric equations, predicted values from the derived generalized 
allometric equation (i.e. Eq. (2), and the raw residuals. Positive raw residuals indicate 
over prediction by the generalized allometric equation whereas negative raw residuals 
indicate under prediction by the generalized allometric equation. The relatively small 
residuals values depicted in Fig. 8 confirms the appropriateness and utility of 
generalized allometric equation to predict the aboveground biomass of standing 
Eucalyptus trees over a very wide range of spatial (location) and temporal (age) 
range. 

The validation of the model was performed using an independent data set (pertaining 
to an experiment on MPTs for energy and boundary plantation conducted at NRCAF, 
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Jhansi) of harvested trees encompassing a dbh range of 9.55 to 26.1 1 cm and biomass 
of 39.30 to 275.92 kg/tree (Table-10) with normality of both tree variates viz dbh and 
AGB (sw. p values < 0.06). 

The residuals were distributed normally with a mean values of -1.00 kg/tree (Table- 
1 1). The plots of residuals against the explanatory variate (dbh) ensures that the 
residuals are not being continuously over or under estimated (Fig. 7). The validation 
of the model was further strengthened by fitting a linear equation (predicted values= 
a+b* observed values) between the predicted and observed values (Fig. 6). For a 
perfect fit model the ‘a’ value should approach to ‘zero’ and ‘b’ value should 
approach to ‘one’ with a perfect unity R^- values. In this case, the value of ‘a’ was 
observed to be 3.17 and ‘b’ as 0.98 with R^- value of 0.98 (Table-12). Although, the 
positive values of ‘a’ as 3.17 indicate slight over estimation by the derived 
generalized biomass equation. However, if we convert this absolute residual values to 
relative residual value: 

Relative residual value= (absolute residual/observed biomass)* 100 

These values (relative residuals) are less than 2% indicating there by the prediction 
errors are relatively very small and thus the derived generalized equations can be 
used with a accuracy prediction level of 98% or more. The strong agreement between 
the derived generalized allometric equation and validation with new data set is further 
sustained by low mean absolute error of 6.348 kg/tree. 

While site-specific allometric equations are likely the best method to calculate 
biomass/ volume, the destructive sampling necessary to derive them is prohibitive in 
terms of time and expense. Logic dictates that the application of site-specific 
allometric equations to sites where they were not developed would not be prudent. 
Thus, the generalized allometric equation derived and validated in this study may be 
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used with confidence by foresters, forest ecologists, and other scientists seeking 
estimates and predictions of volume and aboveground biomass for Eucalyptus in 
India. 

Table-8: Summary characteristics of growth attributes for ABG-dbh simulated dataset. 


Measured 

variable 

No. of 
Cases 

Averse 

CIVIin,Max.) 

Standard 

Deviation 

Skewness 

(SEof 

Skewness 

) 

Kurtosis 

(SEof 

Kurtosis) 

*sw 

Statistic 

SWP- 

Value 

'nbh 

55 

YlAl 

(1.0042.00) 

10.33 

0.8537 

(0.3217) 

0.2313 

(0.6335) 

0.9161 

0.00096 

AGB 

55 

163.63 

(0.260-778.83) 

182.98 

1.7796 

(0.32174) 

2.9267 

(0.6335) 

0.7829 

0.00000 


*SW- Shaphiro wilk’s statistic for normality testing 


Table- 9: AUometric model fitted on the ABG-dbh simulated dataset 


Parameters 

Estimate 

Standard 

Confidence Interval 



Value 

Error 

Minimum 

Maximum 


a 

0.74 

0.255 

0.2291 

1.252 


b 

1.81 

0.098 

1.617 

2.011 

0.91 


Table-10: Summary Characteristics of growth attributes for the independent validation 
dataset 


Measure 

d 

variable 

Average 
(Min, Max.) 

Standar 

d 

Deviatio 

n 

Skewness 

(SEof 

Skewness) 

Kurtosis 

(SEof 

Kurtosis) 

SW 

Statisti 

c 

SWP- 

Value 

Dbh 




fcBIl 

0.9513 

0.03882 

■ 

165.96 

(39.30-275.92) 

68.61 

0.0039 

(0.512) 


0.966 

0.0671 


Table-11: Descriptive statistics of residuals for the ABG-dbh generalized regression equation. 


Average 
(Min, Max.) 

Standard 

Deviation 

Skewness 

(SEof 

Skewness) 

Kurtosis 

(SEof 

Kurtosis) 

SW 

Statistic 

SWP- 

Value 

-1.00 

(-12.60-9.76) 

7.196 

0.130 

(0.512) 

-1.359 

(0.9923) 

0.9333 

0.01792 
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Table-12r ^^dation of aUometric model on new independent dataset by fitting a linear 
equation between predicted and observed ABG values (predicted ABG=a+b*observed ABG). 


Parameter 

Estimate 

Asymptotic 

Standard 

Error 

Wald ConfideBce Interval 

95% 

R" 

Lower 

Upper 

fi 

3.17 

'4.388 ^ 

-6.04 

12.39 

0.98 

— 

0.98 

0.024 

"0.935 

1.03 


Fig. 5: Points generated from individual regressions for Eucalyptus and the resulting 
generalized regression line. 



Fig. 6: Graphs of the linear function fitted between the predicted and observed values 
(predicted biomass= a+b*observed biomass) for Eucalyptus species in India. (For a perfect 
model, ‘a’ should be zero and ‘b’ should be unity along with absolute one R^-value). 
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Fig* 7: Various plots of residual diagnostics on the aboveground biomass validation dataset 
for the allometnc function. 



Autocorrelation Plot 
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j^g. 8t Case wise plot of the predicted and observed values with raw residuals for ABG-dbh 
generalized equation. 



I 

OQ 

S 


Independent validation dataset 


Case Number 


-♦-obsabg 
predabg 
—Jr residual 



Ph.D. Thesis on “DocnmenUtioii, Analysis and Modeling of Eucalyptus.. Jit National Level” by NIGHAT JABEIN (2008) 


150 








Chapter-7: Bevehipmeat af Tree Growth Mo<ielliBg welisitc,-.. 


Development of Tree Growth Modelling Website in Indian Perspective 

and its launching on World Wide Web. 

7.1 Introduction 

The present era has seen an exponential growth and diversification in all forms of 
information, which is sometimes called, as information explosion. It has been made 
possible due to the impact of computer technology on the modem society. 
Computerized mformation systems have influenced nearly all types of organizations, 
whether small or large, public or private, national or multinational. Information 
systems viz (decision support system, expert system etc.) exist for almost all the 
fields. Information systems are concerned with data capture, storage, analysis and 
retrieval. In the context of agricultural in general and agroforestry in particular, they 
are vital to assist decision making in a short time fiame, potentially allowing 
decisions to be made and practices to be auctioned in real time. There are many tree 
species which have been recommended to be adopted in agroforestry systems, 
prominent among them are Eucalyptus, Poplar, Casurina equisitifolia^ Prosopis, etc. 
Eucalyptus species has been widely accepted by the farmers not only in pure 
plantations but also even with more enthusiasm along with crops i.e. in agroforestry 
because of its short rotation cycle, rapid growth and excellent fibre qualities. 
Eucalyptus species has also been recommended as among the six tree species to be 
used for agroforestry in India identified by the Planning Commission through its 
“Task Force on Greening India” that published its report in June, 2001. With the 
increased recognition of fast growing tree species as an alternative to wood 
production, many scientists and development workers have focused their interest on 
short rotation tree species. However, this rapidly rising interest is not matched by a 
corresponding level of available information, which is crucial to their proper 
assessment, evaluation and use in agroforestry research and development. 
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Consequently the task of systematic collection, storage, evaluation and dissemination 
of information on fast growing tree species calls for attention as a matt er of urgency. 

As one of the major problems, peculiar to these species, is the process of data 
collection and management. Attention is drawn to the need to facil i tate dissemination 
of this information. DBMS (Database Management System) is one such tool, which is 
invariably used in almost all displines of forestry/agroforestry and agriculture for 
dissemenation of information through world wide web (Ajit et al, 2004[5]). A 
database provides a conq)uter platform for tabulation, storage arid retrieval of huge 
information and with the rapid development of computers and Intemet, it is now 
possible to access the databases globally which has resulted into availabihty of 
information at fingertip (Ajit et al., 2005[7]). Some databases on forestry and related 
aspects has been developed abroad. Landis et al, (1999)[106] developed Global 
Forest Information Service (GFIS) under the auspices of the Intemational Union of 
Forest Research Organizations (lUFRO) to provide worldwide access to forest 
information. ‘MPTSys’ - a software on multipurpose tree species (Venkateswarlu, 
1994) [2 17], which is a user-fiiendly microcomputer software package has been 
designed to help scientists, administrators and extension workers to organize, manage 
and share research information on multipurpose trees. EUCALIST and TREDAT, 
developed at CSIRO Division of Forestry and Forest Products, Australia, contains 
more than 100000 records on the natural occurrence of Eucalyptus in Australia 
(Brown, 1989) [22]. lADSS-Information and Decision Support System: a user 
fiiendly software package was designed to help scientists for managing field research 
data and other information on multipurpose trees and their potential for producing 
fuelwood and other products (Cady et al, 1989)[26]. Carlowitz (1989)[28] at ICRAF 
led the team responsible for establishing the MPTS database. MPTS database of 
ICRAF was intended to be used as a decision support tool for selection of 
agroforestry tree species (Schroder and Jaenicke, 1994)[174]. It was further revised. 
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refined, modified and renamed as the Agroforestree database (AFX) and released as a 
stand alone application on CD-ROM (Salim et al, 1998)[169]. Rose (1989)[168] 
expressed the serious concern over the need and expected benefits from global 
databases and research networking, with particular reference to growth models of 
MPTS in the tropics. Moreover, databases experiences in tropical forestry are also 
outlined with respect to work done by CATIE (the MIRA database). 

Hann, 2004[69] developed ORGANON, an individual tree growth model for 
Douglas-fir and western hemlock that uses a hst of trees, each with exact 
measurements, as input data. The user can specify periods of growth, in five year 
increments, and management such as thinning, fertilizing, and pruning. The program 
produces stand statistics at each step, and yield tables after final harvest of the stand. 

Kitchen, 2008[92] highlighted exan^les of spatial information collection and 
processing to accomplish real- or near real-time management operations in precision 
agriculture. 

EUis et al., 2005[52] developed SEADSS (Southeastern Agroforestry Decision 
Support System), a web-based application to assist landowners and extension agents 
in the Southeast United States to evaluate potential sites and suitable tree and shrub 
species for agroforestry planning. 

Ellis and Schoeneberger, 2004[53] highlighted the appheation of different Computer- 
based Decision Support Tools (DST) like databases, geographical information 
system, models, knowledge-base or expert systems, and Xiybrid' decision support 
systems in agroforestry research and development. 
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Food and Agricultural organization (FAO, 2007)[55] Fas developed database on 
tropical forages that contains description of more than 600 grassland species and a 
linked picture gallery of photos. 

Infact, the work on development of databases on agroforestry in India was initiated 
by this team itself. Moreover, development of an exclusive site dedicated to tree 
Growth Modelling in India, was initiated as one of the objectives of this thesis work 
itself. This website is probably, the first of its kind in India. An attempt has been 
made in this study to provide maximum possible information on tree growth 
modeling in general and eucalyptus modelling in a hohstic manner. However, no such 
databases could be traced in India related to forestry and agroforestry. This website 
will not only help to avoid duplication of efforts on modelling but also provide a 
single platform, where all pubUshed information on Eucalyptus modelling in India 
would be available at stretch, with an in depth pros and cons. 

7.2 Materials and methods 

Secondary data/information on eucalyptus species has been collected, evaluated and 
synthesized freely from the pubUshed information available in the form of journals, 
books, bulletins, reports, newsletters, articles, leaflets etc. from the 38 centers of All 
India Coordinated Research Project on Agroforestry and ICAR institutes Uke 
NRCAF, IGFRI, CRIDA, CAZRI, CARI, CSSRI, etc and ICFRE institutes/ NGO’s 
etc. The information was compiled, collated and fed into the database. HTML and 
DHTML have been used for creating the dynamic pages of the database at the front 
end. Tree gowth modelling website is a compilation of a biomass and volume 
equations pertaining to eucalyptus published in Indian Context. This website has been 
hosted at the server of Indian Agricultural Statistics Research Institute (lASRI), New 
Delhi and these dynamic web pages are fully fimctional and one can visit at URL 
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htt p://iiiirror4asri.res.iii/net/tfflm/index.html. to view, browse and download 
information on the subject 

7.3 Results and Discussions 

This site entitled Tree Growth Modelling: Indian Experiences" has been designed 
with the basic objective of providing a platform where the user can have the basic 
concepts used in development of tree growth models at one end and the 
comprehensive collection of already developed model for a specific tree species at 
other end so as to use them directly for prediction purposes. There are a number of 
fast growing tree species like Eucalyptus, Poplus etc. that have been widely accepted 
by the Indian farmers in agroforestry/block plantations. The basic reason for adoption 
being the short rotation nature, which provides farmers substantial returns at early 
age. Keeping these considerations in mind, at the initial phase of development of this 
site, we have concentrated on Eucalyptus species modelling. This site is a unique 
collection of published biomass and volume equations, on Eucalyptus, off course, in 
Indian context only. At this stage, as per the objectives of this thesis, we have 
concentrated on eucalyptus only. However, other species adopted for agroforestry in 
India viz Poplar, Casurina, Prosopis, Albizzia etc. will be added in near future. 

The general portion of this site describes definition of modelling by various aufiiors 
and principles of modeling (Fig. 2); types of modelling describing about various 
types of models like descriptive models, prescriptive models, predictive models and 
simulation models (Fig. 3); model development process depicting various steps of 
model development process like model estimation and model validation (Fig. 4); tree 
growth modelling techniques describes essential steps of modelling tree growth like 
sampling of trees, measurement of height, diameter measurement, volume 
conq)utation, height measurement, estimation of biomass components etc. (Fig. 5); 
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models relevant to agroforestry provides the wo rking outlines of the models 
relevant to agroforestry worldwide (Fig. 6); and finally we have added tree growth 
modelling publication displaying a compilation of publication (of our team at 
KRCAF) on tree growth modeling along with full length papers available in pdf 
format for downloading (Fig. 7). 

At present stage we have concentrated only on eucalyptus species as per objective of 
this thesis, Accordingly, in this site under Eucal 5 ^tus introduction (Fig. 8); one can 
find general introduction of eucalyptus and its common uses and importance; Further 
more we have added published biomass equation in India (Fig. 9), like wise published 
volume equations in India along with abstracts (Fig. 10), and lastly a photo gallery 
displaying few snapshots of ongoing Eucalyptus trials at the Central Research Farm 
of NRCAF, Jhansi (Fig. 11). 

This “Tree Growdh Modelling website” has been designed to help for organizing, 
managing and sharing research information Eucalyptus modelling in India. HTML, 
DHTML have been used for creating the dynamic pages of the database at the firont 
end and the source codes written in HTML for designing the home page of Tree 
Growth Modeling Web site has been appended in Annexure 1. The software works 
on a HTTP Server and cater the requirement of client that may on any computer 
connected with Internet and having a graphic web browser. The unified database has 
been equipped witii interactive web pages to browse information on Eucalytus 
introduction, modelling introduction, types of modelling, model development 
process, detailed information on biomass and volume models in India related to 
Eucalyptus! species. 
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Fig. 1: Home page of Tree Growth Modelling Website 
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Annexures 


Annexure. 1: HTML source code of the homepage of Tree Growth 
Modelling Website. 
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BLOCK PLANTATIONS IN SEMI-ARID CONDITIONS OF INDIA 

AJIt, Ritu, Ntgiiat Jabeen, Rajender Singh, A.K.i^n(ta, Ram NewaJ and O.P.Chaturvedi 

N^)rK^ Research Cen^ for Agro forestry 
GwaBorRoa(l.Jhai^m003(UP),lndBa 

E.maH: uma^it^fria@ysdvM.(XJ.in; umaaj^3Uf^a@redffinail.com, aji^maS.nrcaf.emet.ki 

p experiment for growth and yield anaiysis/modeling of Eucalyptus based systems in Bundelkhand 
fiion was initiated at NRCAF, Jhansi with four clones of EucalyfOus terehcomis (C-3. C- 6 , C-7 and 
^0 obtained from ITC Bhadrachalam, Andhra Pradesh) in different systems and spacing, namely 
liisilviculture (5x4, 10x2, 10x5, 8x4, 5x5 m). compact block (3x3, 2.5x2.5 m) and boundaiy plan- 
ons (2.5m). The growth oteervations were rKMrded atter every three months tti ^e Myem 
I after six month in second year onwards. Observed parameters include total heig . ’ 

1 canopy of the each individual tree. Atotal of 1335 trees from the complete ^ 

tsured for the find year. However, from second year only 51 0 trees were ^ 

leriment for further growth observation. Eight trees (two from each w , 
npact block plantation of spa cing 2.5 x 2.5 m and other eight were harvested from agnsrlvicult ^ 

it SYMPOSIUM ON AGROFORESTRY UVHJHOOOS^I^. ^RON«^ PROTECTJON AN^ ^ 

. PROOIK:TION (IS-IS t^CEMBER. 2006) WRC-AGROFORESTRY, JHANSI (U.P.) 
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Height Estimatioii Model for Eucalyptus tereticomis Grown ynder 

Semiarid Conditions of India 

AJit, RStu Srivastava, Nighat Jabeen, A.K. Handa, Rajendra Singh and 

O.P. Chaturvedi 

National Research Centre for Agro forestry 
Gwalior Road, Jhansi-284003 India 


abstract 

An experiment on structural analysis and modelling of Eucalyptus based system was initiated at NRCAF, Jhansi in 2003 
with three systems (agrisilviculture, block plantation and boundary plantation) and four <*)nes obtained from ITC 
Bhadrachalam. To estimate the height of large trees, three types of models viz Richards, Schumacher, Allometric vrere 
attempted with dbh as predictor and the Revalues were almost comparable In ttie three systems. However, when these 
models were extrapolated to examine their behaviour outside the observed range of dbh, the allometric function resulted 
in reasonably acceptable predictions even quite outside the observed range also, whereas the othertwo namely Richards 
and Schumacher leads to merely constant estimation for the extrapolated region. Accordingly, the allometric models are 
proposed for estimating the height values by simply measuring the dbh of the standing Eucalyptus trees in these systems. 

Key words: Eucalyptus, height diameter models, allpmetric equations, valids^m 


1. INTRODUCT»ON 

In forestry, mathematical equations form the basis 
of many prediction studies. Typical example of 
such equation are height prediction model based 
on diameter, such models gain importance as 
modem scientific management relies heavily on 
p well defined growth and yield models that can be 
used to access the status of stands at any point of 
time. Eucalyptus tereticomis being exotic to India, 
its cultivation has spread to nearly the whole of the 
Indian subcontinent. Large-scale plantations of E. 
tereticomis were established between 1970 and 
1985. Nearly 2% of the cultivable land in 
northwestern India is planted with E. 
teref/corn/s.The average productivity of E. 
tereticomis plantations is 10 mVha/year on 
forestlands and 15-20 m®/ha/year on farmlands. 
The average productivity of commercial clones is 
about 20-25 m%a/year and many farmers have 
achieved up to 50 m®/ha/year. Infact, for the 
computation of productivity of such plantations, 
precise measurements of growth attributes (height 
and diameter) are essential. It has been noticed 
during the couree of experimentation that it 
becomes extremely tedious to accurately measure 
the height of large trees (more than 10 meter), 
whereas the dbh (diameter at 1 .37m height) can be 
easily measured and accordingly height-diameter 
models provide a handy tool for prediction of height 
4 on the basis of dbh values. Recently height- 
diameter curves have been frequency developed 
in forestry studies (Sharma and Parton, 2007; 
Lootens et ai., 2007; Casteo ef a/., 2006; Reed et 
ai., 2003; Mehtatalo, 2005; Inoue and Yoshida, 
2004; Sochacki et a/, 2007). Bachpai et ai , 2005 
developed multiple regression models for 


prediction of height using dbh and age for 
Bambusa tuida Ro)d3. plus clumps in regions of 
Assam and Meghalaya. Tiwari and Gadow, 2003 
developed a height-diameter relationship for 
Prosopis cineraria for the hot arid areas of Indra. 
Tewari et ai., 2002, have attempted height-age 
and diameter-age for irrigated plantation of 
Rajasthan; Soares and Tome, 2002, have 
developed height-diameter equation for 
Euca/ypfus in first rotation. 

Present study deals with the development of 
height-dbh model under different systems, utilizing 
the data of initial years of growth, for application to 
a much wider range (viz maturity and even 
senescence age). 

2. MATERIALS AND METHODS 

2.1 Study site and experimental details 

The experiment is being carried out at the 
Central Research Farm of National 
Research Centre for Agroforestry, Jhansi 
which is located at an elevation of 300 m 
above sea level and is situated between 
24”11' N latitude and 78‘“17' E longitude 
having tropical semi arid climate with mean 
annual rainfall of 900 mm. More than 75% of 
the rainfall is received during monsoon 
season (last week of June to first week of 
September). Mean maximum temperature 
ranges from 47.4° C (June) to 23.5° C 
(December) and mean minimum 
temperature from 27.2° C (June) to 4.1° C 
(December). The soil of the exfterimental 
area is sandy loam with low organic carbon, 
nitix»gen, phosphorous and medium to high 
in potassium. The study was initiated in 
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August, 2003 with the objective of studying 
the structure, biomass, productivity and 
modelling the growth of Eucalyptus based 
systems. The experiment was laid with three 
systems of Eucalyptus tereticornis namely 
agrisilviculture (AS), block plantation (CB) 
and boundary plantation (BP). Tree spacings 
considered in agrisilviculture are 5x4 m, 
10x2 m, 10x5 m, 8x4 m and 5x5 m; compact 
block are 3x3 m and 2.5x2.5 m and boundary 
plantation is 2.5 m. Wheat and Black gram 
were grown in the interspaces during Rabi 
and Kharif seasons respectively under 
boundary plantations and agrisilviculture 
systems. Four Eucalyptus clone namely C-3, 
C-6, C-7 and C-10 were obtained from ITC, 
Bhadrachalam (A.P.) and planted in field 
during August, 2003. In this article, the 
development and validation of height- 
diameter models have been reported. 
Growth observations (dbh-diameter at 
breast height and height) were recorded 
quarterly in the first year and twice a year in 
the subsequent period. 

2.2 Data and observations 

The available data includes the tree wise 
observations on height in meter taken from 


the ground level to tip of the shoot of the main 
bole and dbh (diameter at breast height) in 
centimeter taken at a height of 1 .37m at the 
main bole) for a period of 28 months. During 
the first year, observations were recorded 
four times in a year on 1 335 trees distributed 
over the three systems viz agrisilviculture, 
compact block and boundary plantation. 
Whereas from second year onwards 
observations were recorded twice a year on 
510 selected trees distributed proportionally 
over the three systems. 

2.3 Statistical analysis and candidate 
function 

Systat software (Wilkinson etaL, 1996) was 
used for computation of descriptive statistics 
(mean, SD, skewness, kurtosis etc) and 
fitting of non-linear equations (estimates of 
model parameters, asymptotic standard 
error of estimate, confidence interval, 
values) and plotting of various graphs 
pertaining to residual diagnostics and model 
validation (probability plot of residuals, auto 
correlation plots etc.) 

The following functions were attempted to fit the 

observed data 


Name of function 

Functional form 

Parameter interpretation 

Allometric 

Height = a*(DBH)'’ 

a= Scaling parameter 
b = Power coefficients 

Richards 

Height =a*[1-exp(-k*DBH)l 

a = Asymptote 

k = Rate of approach to the asymptote 
m = Shape parameter 

Schumacher 

Height=1 .37+a*exp(-b/DBH) 

a = Asymptote 
b= Shape parameter 


3. RESULTS AND DISCUSSION 

The observed values for height and dbh exhibited 
the approximate normality of variates (Table-1) 
under study. To get an idea of the shape of the 
function to be fitted on the data, a scatter plot of 
height vs dbh was initially drawn. It was clear from 
this scatter plot that the candidate functions 
usually adopted for modeling height-dbh curves 
viz Richards, Schumacher ar»d Allometric/power 
will fit well the observed pooled data 
set(AS+CB+BP). 

ThepaiaTieiaresimgtesatongwitholherr^stedslalisfi^ 
of the Richaixfe ftrrdion lied on pooled date has been 
conpied h Tstte 2, axl the lesiitent eoj^ion was 
Height=18.03*l1-exp {0.089*dbh)l with R" 

(obs. vs pred.) = 0.^. Similarly ttie parameter 
estimates of schumax^r and alkMTretric functions 
have been compiled in Tatfe 2, and the resultant 
equations were He^ =1 .37 + 1 8.95*exp (-6.79/dbh) 


with R^ (obs. vs pred.) = 0.94 for sdiumadier and 
Height = 1 .76*(dbh)‘’™ with R^ (obs. vs pred.) = 0.97 
for allometric fonctions resp^diveiy. The plots of 
these fitted functions along with the observed data 
points have been comfxlai in Fig.1 . It is dear from 
the above computaticwts and the grafrfrs of tire fitted 
functions that the R^ -values are comparafcrfe for ail 
the three functions but was obsereved to be 
maximum for allometric one. However, as has been 
pointed out by ottier researchers also in the recent 
artides (/^it et at., 2006; Prajneshu and Chandran, 
2005), that R^' value alone should not be used to 
judge the. best fitting function and it is equally 
important to consider the validation and more 
importantly the behavkx of the fitted fonction outside 
the observed range of the independent variate. 
Accordingly, to jut^ ttie pn^liction capabilities of 
ttiese functions in the extrapolated region, the 
function curves for the three-fitted function were 
drawn (Fig. 2). 
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Table 1 . Summary characteristics of growth attributes on the observed data set. 


Measured variable 

Average 
(Min, Max.) 

Standard 

Deviation 

Skewness 

SEof 

Skewness 

Kurtosis 

SEof 

Kurtosis 

Height 

6.09 

(1.08-14.90) 

3.615 

0.689 

0.045 

-0.55 

0.090 

dbh 

5.04 

(0.20-14.33) 

3.690 

0.530 

0.045 

-0.84 

0.090 


Table 2. Estimation of parameters of various functions fitted on pooled dataset (AS+CB+BP). 


SN 

Fitted Function 

Parameter 

Estimate value 

Standard Error 

Confidence interval 

Min Max 

Revalue 

1 

Wlometric 

A 

1.76 

0.015 

1.73 

1.79 

0.97 



B 

0.78 

0.004 

0.78 

0.79 


2 

Richards 

A 

18.03 

0240 

17.50 

18.50 

0.95 



K 

0.089 

0.001 

0.08 

0.09 




M 

0.011 

0.009 

0.00 

0.03 

I 

. i 

' ' ■' ! 

3 

Schumacher 

A 

18.95 

0.150 

18.60 

19.20 

'■ 0.94 I 



B 

6.79 

0.060 

6.66 

6.91 

i 

I 




Fig^. Extrapolated pfe#::tkm of haght ^owth with r^pect to DBH for diftorent functions fitted on pooled dataseti 
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It is very clear from these comparative extrapolated 
prediction graphs of the three fitted functions that 
the allometric function results in reasonably 
acceptable estimations even quite outside the 
observed range also, whereas the other two 
namely richards and schumacher leads to merely 
constant estimation of size for the extrapolated 
range. Therefore the allometric function which 
meets both, the criteria's of high R^-value and 
reasonably acceptable extrapolated predictions, 
was preferred over the other two. 

To judge the accuracy of prediction, the developed 
allometric equation was thoroughly evaluated 
through residual diagnostics. The error of 


prediction termed as residual is computed as the 
difference in the observed and predicted values. 
Theoretically, the residual should be independently 
and normally attributed with mean zero and 
constant variance. These assumptions were 
evaluated through pertinent graphs (Fig.3). The 
plot of residuals against their expected values 
clearly portrayed the normality of residuals. The 
plot of residual vs. independent variate confined 
that the residual are not continuously being 
over/under estimated and the plot of residual 
against estimate indicated that the residual have 
constant variance. Thus the proposed allometric 
equation fulfilled the regression requisites. 



Fig. 3. Different plots of residual diagnostics for the allometric function fitted to the pooled dataset. 


The allometric function for the pooled data 
(including all the three systems i.e AS+CB+BP) 
resulted in He^ht = 1.76*(dbh) with (obs. vs 
pred.)=0.97. To compare the three systems height- 
DBH models based on allometric Unctions were 
also fitted separately system wise (Fig.4) and the 
resultant equations (Tabie-3)are: 

Height = 1.87*(dbh) with R" =0.96 for 
agrisilviculture 

Height = 1 .61*(dbh) ““ with R" =0.95 for impact 
block 

Height = 1.65*(dbh)°‘“' with R^ =0.97 for boundary 
plantation respectively 


To compute the growth performance in totality, all 
the three fitted allometric curves (agrisilviculture, 
compact block and boundary plantation), were 
drawn in one single frame (Fig. 5), and it is clear 
from this plot that growth was better in 
agrisilviculture for the first year, whereas compact 
block and boundary plantation picked up the 
growth in the second year. The proposed models 
can be used for predicting the height of standing 
trees at any stage without destructive harvesting, 
by simply measuring the diameter atbreast height. 


Table 3. Estimation of parameters for the allometric function fitted on individual datasets under three 
systerTK. 


S.N. 

System 

Parameters 

Estimate 

Value 

Standard 

Error 

Confidence Interval 
Minimum Maximum 

R2 

1 

Agrisilvicuiture 

A. 

1.87 

0.020 



0.96 



B. 

0.75 

0.005 




2 

Compact Block 

A. 

1.61 

0.030 



0.95 



B. 

0.83 

0.009 




3 

Boundary Plantation 

A. 

1.65 

0.030 



0.97 



B. .'.i 

' 0.81 . 

0.010 
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Fig. 4. MIometric function fitted on individual datasete under the ttire© systann*. 



Fig 5. Coinpartson of the fitted sdoineMc ec^iaion in dHlsiinit sy^enns 


All used functions (allometric. richarcte and 
sdvumacher function) bseed on the crista of Ri* 
vs pred.) values, gave conparEtele rrautts. 
The allometric fum^ticm lead to the reasonably 
acceptable predictkxis when these functions were 
validated on the criteria of extrapolated estfinattons 
outside the observed range followed by 
Schumacher function and richards’s funr^on 
ranked at the last Since model selection usteg 
nonlinear regression is an inherently subjective 
process, and somewhat data dependent bas^ 
the final evaluation of the model on the va^ctefon 
procedure through pertinent graphs and 
extrapolated prediction seems to be a viatote 
alternative when the prime objective of nrodefing 
exercise is prediction. 
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